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<?xml version="1.0" encoding="UTF-8" ?> 


Ө <!DOCTYPE plist PUBLIC "-//Apple Computer//DTD PLIST 1.0/ 
JEN" "http://www 


apple.com/DTDs/PropertyList-1.0.dtd"» 


Ө <plist version="1.0"> 
Ө <dict> 
Ө <key>Major Version</key><integer>1</integer> 
«Кеу»Міпог Version</key><integer>1</integer> 
-- snip -- 
e <key>Tracks</key> 
<dict> 
<key>2438</key> 
<dict> 
<key>Track ID</key><integer>2438</integer> 
<key>Name</key><string>Yesterday</string> 
<key>Artist</key><string>The Beatles</string> 


<key>Composer</key> 
<string>Lennon [John], McCartney [Paul]</string> 


<key>Album</key><string>Help!</string> 
</dict> 
-- snip -- 
</dict> 
о <key>Playlists</key> 
<array> 
<dict> 
<key>Name</key><string>Now</string> 


<key>Playlist ID</key> 


<integer>21348</integer> 
-- snip -- 

<array> 

<dict> 


<key>Track ID</key> 
<integer>6382</integer> 


</dict> 
-- snip -- 
</array> 
</dict> 
</array> 
</dict> 


</plist> 
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def findDuplicates(fileName): 


print('Finding duplicate tracks in %s...' % fileName) 


# read in a playlist 


e plist = plistlib.readPlist(fileName) 
# get the tracks from the Tracks dictionary 
Ө tracks = plist['Tracks' | 
# create a track name dictionary 
e trackNames = {} 
# iterate through the tracks 
ο for trackId, track in tracks.items(): 
try: 
e name = track['Name'] 
duration - track['Total Time'] 
# look for existing entries 
о if name in trackNames: 
# if a name and duration match, increment 
the count 


cond 


# round the track length to the nearest se 


ο if duration//1000 == trackNames [name] 
[0]//1000: 


count = trackNames[name] [1] 


Ө trackNames[name] = (duration, count+1 
else: 


# add dictionary entry as tuple (duration, 
count) 


© trackNames[name] = (duration, 1) 
except: 
# ignore 


pass 
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Я store duplicates as (name, count) tuples 
e dups - [] 
for k, v in trackNames.items(): 
e if v[1] > 1: 
dups.append((v[1], k)) 
# save duplicates to a file 
e if len(dups) > 0: 
rint("Found %d duplicates. Track names saved to d 
up.txt" % len(dups)) 
else: 
print("No duplicate tracks found!") 
ο Т = open("dups.txt", "w") 
for val in dups: 
e f.write("[%d] %5\п" % (val[0], val[1])) 


f.close() 
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def findCommonTracks(fileNames) : 
# a list of sets of track names 
9 trackNameSets = [] 
for fileName in fileNames: 
# create a new set 
e trackNames - set() 
# read in playlist 
e plist = plistlib.readPlist(fileName) 
# get the tracks 
tracks = plist['Tracks' | 
# iterate through the tracks 
for trackId, track in tracks.items(): 
try: 
# add the track name to a set 
о trackNames.add(track['Name' 1) 
except: 


# ignore 


pass 
# add to list 
e trackNameSets.append(trackNames) 
# get the set of common tracks 
о commonTracks = set.intersection(*trackNameSets ) 
# write to file 
if len(commonTracks) > 0: 
ο f = open("common.txt", "w") 
for val in commonTracks: 
S = "%s\n" $5 val 
Ө f.write(s.encode("UTF-8")) 
f.close() 
print("%d common tracks found. " 
"Track names written to common.txt." % len(c 
ommonTracks)) 
else: 


print("No common tracks!") 
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def plotStats (fileName): 
# read in a playlist 
e plist = plistlib.readPlist(fileName) 
# get the tracks from the playlist 
tracks = plist['Tracks' | 
# create lists of song ratings and track durations 
e ratings = [] 
durations - [] 
4 iterate through the tracks 
for trackId, track in tracks.items(): 
try: 
e ratings.append(track['Album Rating']) 


durations.append(track['Total Time']) 


ехсерї: 
# ignore 


pass 


# ensure that valid data was collected 
о if ratings == [] or durations == |1: 
print("No valid Album Rating/Total Time data in %s 
." % fileName) 


return 
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4 scatter plot 

e х = np.array(durations, np.int32) 
# convert to minutes 

e x = x/60000.0 

e y = np.array(ratings, np.int32) 


о pyplot.subplot(2, 1, 1) 


pyplot.plot(x, у, 'о') 
pyplot.axis([0, 1.05*np.max(x), -1, 110]) 


pyplot.xlabel('Track duration' ) 


o © о © 


pyplot.ylabel('Track rating') 


# plot histogram 
pyplot.subplot(2, 1, 2) 

© pyplot.hist(x, bins=20) 
pyplot.xlabel('Track duration' ) 


pyplot.ylabel('Count') 


# show plot 


® pyplot.show() 
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def main(): 
# create parser 
descStr = """ 


This program analyzes playlist files (.xml) exported f 
rom iTunes. 


e parser = argparse.ArgumentParser(description=descStr) 
# add a mutually exclusive group of arguments 


о group = parser.add mutually exclusive group() 


# add expected arguments 


group.add argument('-- 
common', nargs='*', dest='plFiles', required=False) 


group.add argument('-- 
stats', dest-'plFile', required-False) 


Ө group.add argument('-- 
dup', dest='plFileD', required=False) 


# parse args 


о args = parser.parse args() 


if args.plFiles: 
# find common tracks 
findCommonTracks(args.plFiles) 
elif args.plFile: 
# plot stats 
plotStats(args.plFile) 
elif args.plFileD: 
# find duplicate tracks 
findDuplicates(args.plFileD) 
else: 


print("These are not the tracks you are looking f 
or.") 
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εγ/ρ!8γ!151/ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
import re, argparse 
import sys 
from matplotlib import руріої 
import plistlib 


import numpy as np 


def findCommonTracks(fileNames) : 
Find common tracks in given playlist files, 
and save them to common.txt. 
# a list of sets of track names 
trackNameSets = [] 
for fileName in fileNames: 
# create a new set 
trackNames = set() 
# read in playlist 
plist = plistlib. readPlist (fileName) 
# get the tracks 
tracks = plist['Tracks' | 
# iterate through the tracks 


for trackId, track in tracks.items(): 


try: 
# add the track name to a set 
trackNames.add(track['Name']) 
except: 


# ignore 


pass 


# add to list 
trackNameSets.append(trackNames ) 
# get the set of common tracks 
commonTracks = set.intersection(*trackNameSets) 
# write to file 
if len(commonTracks) » 0: 
f = open("common.txt", 'w') 
for val in commonTracks: 
S = "%s\n" % val 
f.write(s.encode("UTF-8")) 
f.close() 
print ("%d common tracks found. " 


"Track names written to common.txt." % le 
n(commonTracks)) 


else: 


print("No common tracks!") 


def plotStats(fileName): 


Plot some statistics by reading track information f 
rom playlist. 


# read in a playlist 


plist = plistlib.readPlist(fileName) 


ce 


o9 


# get the tracks from the playlist 
tracks = plist['Tracks' | 
# create lists of song ratings and track durations 
ratings = [] 
durations = [] 
# iterate through the tracks 
for trackId, track in tracks.items(): 
try: 
ratings.append(track['Album Rating']) 
durations.append(track['Total Time']) 
except: 
# ignore 
pass 
# ensure that valid data was collected 
if ratings -- [] or durations -- []: 


print("No valid Album Rating/Total Time data in 
fileName) 


return 


# scatter plot 
X= np.array(durations, np.int32) 
# convert to minutes 


х = x/60000.0 


йет 


у = np.array(ratings, np.int32) 
pyplot.subplot(2, 1, 1) 

pyplot.plot(x, у, 'о') 

pyplot.axis([0, 1.05*np.max(x), -1, 1101) 
pyplot.xlabel('Track duration') 


pyplot.ylabel('Track rating') 


# plot histogram 
pyplot.subplot(2, 1, 2) 
pyplot.hist(x, bins-20) 
pyplot.xlabel('Track duration') 
pyplot.ylabel('Count') 

# show plot 


pyplot.show() 


findDuplicates (fileName): 


Find duplicate tracks in given playlist. 


print('Finding duplicate tracks іп 55...' % fileNam 


# read in playlist 


plist = plistlib. readPlist (fileName) 


# get the tracks from the Tracks dictionary 
tracks = plist['Tracks' | 
# create a track name dictionary 
trackNames = {} 
# iterate through tracks 
for trackId, track in tracks.items(): 
try: 
name = track['Name'] 
duration = track['Total Time'] 
# look for existing entries 
if name in trackNames: 


# if a name and duration match, increme 
nt the count 


# round the track length to the nearest 


second 
if duration//1000 == trackNames [name] 
[0]//1000: 
count = trackNames[name] [1] 
trackNames[name] = (duration, count 
+1) 
else: 
# add dictionary entry as tuple (durati 
on, count) 


trackNames[name] = (duration, 1) 


except: 


# ignore 

pass 
# store duplicates as (name, count) tuples 
dups = [] 
for k, v in trackNames.items(): 

Lf v[1] > 1: 

dups.append((v[1], k)) 
# save duplicates to a file 
if len(dups) > 0: 


print("Found %d duplicates. Track names saved t 
o dup.txt" % len(dups)) 


else: 
print("No duplicate tracks found!") 
f = open("dups.txt", 'w') 
for val in dups: 
f.write("[%d] %sNn" % (val[0], val[1])) 


f.close() 


# gather our code in a main() function 
def main(): 

# create parser 

descStr = """ 


This program analyzes playlist files (.xml) exporte 


d from iTunes. 


parser = argparse.ArgumentParser(description=descSt 


# add a mutually exclusive group of arguments 


group = parser.add mutually exclusive group() 


# add expected arguments 


group.add argument ( ' - - 
common', nargs='*', dest='plFiles', required=False) 


group.add argument ( ' - - 
stats', dest='plFile', required=False) 


group.add argument ( ' - - 
dup', dest='plFileD', required=False) 


# parse args 


args = parser.parse args() 


if args.plFiles: 
# find common tracks 
findCommonTracks(args.plFiles) 
elif args.plFile: 
4 plot stats 
plotStats(args.plFile) 


elif args.plFileD: 


# find duplicate tracks 
findDuplicates(args.plFileD) 
else: 


print("These are not the tracks you are looking for 


# main method 


if name == ' main ': 


main() 


1.5 [I 
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$ python playlist.py --common test-data/maya.xml test- 
data/rating. xml 


000000 


5 common tracks found. Track names written to common.txt. 
$ cat common. txt 

God Shuffled His Feet 

Rubric 

Floe 

Stairway To Heaven 

Pi's Lullaby 


moksha:playlist mahesh$ 


ОООО0000000000000000 


$ python playlist.py --stats test-data/rating. xml 


Π1-1ΠΠΠΗΠΗΠΗΠΗΠ 
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• ПО000000 

• ПОД000000000 
• ДОДО)000000 
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ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΕωΥΤΙΕΓΙΠΠΠΠΠΠΠΠΠΗ 
LOEO euren nan p nnnm n a d pd d ECC 
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ΕΕ ΕΕ ΕΕΣ ΕΕ 7 ООООО000000000000000Ш0 X 0 
у ПОП 


ПИ ИП 
x = r cos( 0) 
y=rsin( Ө) 


000000000000000 Ө ΠΠΠΠΠΠΠΠΠΠΧΠΥΠΠΠΠΠΗ 
0000000000000Х *+ Y*= R #0000 e 0000 


πμ. x ПО y 0002-2000000 





x = rcos|6) 
y = rsin(6) 


02-2 ΠΠΠΠΠΠΠΠ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΙ( 8, 2)0 

0000000 xy 000000000000000000000 х = а 

+ rcos( 0 )[] у = b + rcos( Θ)ΠΠΠΠΠΠΠΠΠΠΠΗ 
000000 


2.1.1 ΠΠΠΠΠ 


02-ЗОДОДОО000000000000000000000000000 
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ΠΠΠΡΟΠΠΠΠΠ ΓΠΠΠΠΠΠΙΠ! = PC / ΓΠΠΠΠΠΠΠΠΠ 
ООО00000000000000000РО0000000000000 


1 = k 
r—R (а — k) со8( 0) + lkcos (--»)) 
y=R (а — k)sin(8) + 1k sin ( п г)) 
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ООО000000000000000000000000000000000 
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О000000000000006СорО000000000000000 
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О0000013000000000004А000000000000000 
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2.1.2 II 
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ОООО0000000000000000 


О0ро00000вигбтеррородоррордроп 
drawcircle.py[][]Python[]HLLL 


import math 


Ө import turtle 


à draw the circle using turtle 
def drawCircleTurtle(x, y, r): 


# move to the start of circle 


ө turtle.up() 
© turtle.setpos(x + г, у) 
о turtle.down() 


# draw the circle 
© for i in range(0, 365, 5): 
о a = math.radians(i) 
ο 


turtle.setpos(x + r*math.cos(a), у + r*math.sin(a 


)) 


@ drawCircleTurtle(100, 100, 50) 


© turtle.mainloop() 


ПӨПППППЕч Н ероооооооооО 
ДдгаумуСїгс!еТиге()ПППППӨПППирОПППП 
ΡγΥΓΠΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
О00000000000 


Le pp pa pp Dp D bx + г, у)000(х0 
У)О000000000000000000өП009омп()00ө 
О000гапое(о, 365, 5)0000005000000010000 
360ΠΠΠΙΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘΠΠΠΠΠΠΠΠ 
ОСО000000000000000000 


ОФПОДДОДо0000000000000000000000000 
ОДОО00000000000000000000000000МО0000 
ОО000000000300000000000000000 


ПӨПІ ППағгау/СігсіеТигйео ПП ППӨППІІ 


mainloop ОДООЄкіпеегобДОДОДООДОООООО 
[]|Tkinter(]Python[][]T]GUI[]E]L] 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
2.2 ПЦ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
‚ turtle ПП 


. FE ERIK H Tri; 
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ΠΠΠΠΠΠ51ΡΓΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
агам()0000000000000чрааѓе() 0000000000 


ΠΠΠΠΠΠΠΠΠ5ΡΙΓΟΠΠΠΠΠΠΠΠΠΠΠΠΗ 
SpiroAnimator{][] 


0000000000000000 2. 400 
2.3.1 Spiro 
П005ріго00000 


# а с1аѕ5 that draws a Spirograph 
class Spiro: 
# constructor 


def init (self, xc, yc, col, В, г, 1): 


# create the turtle object 
e self.t = turtle.Turtle() 

# set the cursor shape 
о self.t.shape('turtle') 

4 set the step in degrees 


e self.step = 5 


# set the drawing complete flag 


о self.drawingComplete = False 


# set the parameters 


e self.setparams(xc, yc, col, R, г, 1) 


# initialize the drawing 


о self.restart() 


ПӨП05ріго0000000000чеїероооооо000000 
ΠΠΗΠΗΠΠΘΙΠΠΠΗΠΠΠΗΠΠΠΗΠΗ 
https://docs.python.org/3.3/library/ 
turtle.htm! [T] ]r]turtten rin np mnn mm e L) 
ПОСОООООООО5000еОО00000000000000000 
ПОП 


ОӨПӨППО000000000000000 
2.3.2 DUDO 


μμ αλά гагп5()ПППППППЦПП5р!го[ППШП 


4 set the parameters 





def setparams(self, xc, yc, col, R, г, 1): 
# the Spirograph parameters 


e self.xc = xc 


e 


ο 


6 


self.yc = yc 
self.R = int(R) 
self.r = int(r) 
self.l = l 


self.col = col 


# reduce r/R to its smallest form by dividing with 
the GCD 


gcdVal = gcd(self.r, self.R) 
self.nRot = self.r//gcdVal 
# get ratio of radii 

self.k = r/float(R) 

# set the color 
self.t.color(*col) 

# store the current angle 


self.a = 0 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘΠΠΠΠΠΠΠΠΠΠ ΑΒ ΠΓΠΗ 
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2.3.3 restart()|_ |] 
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# restart the drawing 
def restart(self): 
4 set the flag 
e self.drawingComplete - False 
# show the turtle 
e self.t.showturtle() 


# go to the first point 


Ө self.t.up() 
ο А, k, 1 = self.R, self.k, self.l 

а = 0.0 

x = R*((1-k)*math.cos(a) + l*k*math.cos((1- 
k) *a/k) ) 

y = R*((1-k)*math.sin(a) - l*k*math.sin((1- 
k) *a/k) ) 
о self.t.setpos(self.xc + x, self.yc + у) 
ο self.t.down() 


ΠΠΗΠΠΠΗΠαΓαννίπο(οπιρ]εῖε[1ΠΠΠΠΠΠΠΠΠΠΠΗ 
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2.3.4 агам() || 


агам()00000000000000 


# draw the whole thing 
def draw(self): 
# draw the rest of the points 
R, К, 1 = self.R, self.k, self.l 
e for i in range(0, 360*self.nRot + 1, self.step): 
a = math.radians(i) 


x = R*((1-k)*math.cos(a) + l*k*math.cos((1- 


k)*a/k)) 
y = R*((1-k)*math.sin(a) - l*k*math.sin((1- 
k)*a/k)) 
self.t.setpos(self.xc + x, self.yc + y) 
# drawing is now done so hide the turtle cursor 
e self.t.hideturtle() 


ΠΘΙΙΓΠΙΠΠΠΠΠΙΠΠΗΗΠΠΠΗΠΠΠΠΗΞ60ΠΠηΒοϊ/1ῃ 
шт" XU У ЮДДФеДОДОДОДООО 


2.3.5 ΠΠΠΗ 
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4 update by one step 
def update(self): 


# skip the rest of the steps if done 


e if self.drawingComplete: 
return 
4 increment the angle 
6 self.a += self.step 
# draw a step 
R, k, L = self.R, self.k, self.l 
# set the angle 
e a - math.radians(self.a) 


x= self.R*((1-k)*math.cos(a) + l*k*math.cos((1- 


k)*a/k)) 

у = self.R*((1-k)*math.sin(a) - l*k*math.sin((1- 
k)*a/k)) 

self.t.setpos(self.xc + x, self.yc + y) 

# if drawing is complete, set the flag 
ο if self.a >= 360*self.nRot: 


self.drawingComplete = True 
# drawing is now done so hide the turtle curso 


self.t.hideturtle() 
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2.3.6 SpiroAnimator]|] 


зрігоАпіплагоГоДОООООО0000000000000000 
ООООО0000000000000000000000000000000 
ОООДО0000000000000000000000019000 


# a class for animating Spirographs 
class SpiroAnimator: 
# constructor 
def init (self, М): 
# set the timer value in milliseconds 
о self.deltaT = 10 
# get the window dimensions 
e self.width = turtle.window width() 
self.height = turtle.window height() 
# create the Spiro objects 
e self.spiros = [] 
for i in range(N): 
# generate random parameters 


о rparams = self.genRandomParams() 


4 set the spiro parameters 

Θ Spiro = Spiro(*rparams) 
self.spiros.append(spiro) 
# call timer 


о turtle.ontimer(self.update, self.deltaT) 
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орорроФеооодор0000000 
genRandomParams ПОДОО00000000 


2.3.7 депКапаотРагатѕ()[[ 


П00депАапаотрРагатѕ()О000000000000 
Spiro ПП ПП 


# generate random parameters 
def genRandomParams (self): 


width, height = self.width, self.height 


e R = random.randint(50, min(width, height)//2) 
e r = random.randint(10, 9*R//10) 
e l = random.uniform(0.1, 0.9) 
о хс = random.randint(-width//2, width//2) 
e yc - random.randint(-height//2, height//2) 
о col = (random.random(), 
random. гапдот(), 
random. random() ) 
ο return (xc, yc, col, R, г, 1) 
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# restart spiro drawing 
def restart(self): 


for spiro in self.spiros: 


# clear 

spiro.clear() 

# generate random parameters 
rparams - self.genRandomParams() 
4 set the spiro parameters 
spiro.setparams(*rparams) 

# restart drawing 


spiro.restart() 


О000005рігорообООО000000000000000000 
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2.3.9 update()[][] 
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def update(self): 
# update all spiros 
e nComplete - 0 
for spiro in self.spiros: 
# update 
e spiro.update() 
# count completed spiros 
e if spiro.drawingComplete: 


nComplete += 1 


4 restart if all spiros are complete 
о if nComplete == len(self.spiros): 
self.restart() 
# call the timer 


e turtle.ontimer(self.update, self.deltaT) 
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я toggle turtle cursor on and off 
def toggleTurtles(self): 
for spiro in self.spiros: 
if spiro.t.isvisible(): 
spiro.t.hideturtle() 
else: 


spiro.t.showturtle() 


2.3.11 0000 
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# save drawings as PNG files 
def saveDrawing(): 
# hide the turtle cursor 
e turtle.hideturtle() 
# generate unique filenames 


e dateStr = (datetime.now()).strftime( "9sd9sb9sY - 


fileName = 'spiro-' + dateStr 
print('saving drawing to %s.eps/png' % fileName) 
# get the tkinter canvas 
e canvas - turtle.getcanvas() 
# save the drawing as a postscipt image 
о canvas.postscript(file = fileName + '.eps') 
# use the Pillow module to convert the postscript 
image file to PNG 
e img = Image.open(fileName + '.eps') 
о img.save(fileName + '.png', 'png') 
# show the turtle cursor 


[7] turtle. showturtle() 
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e parser = argparse.ArgumentParser(description=descStr) 


# add expected arguments 
e parser.add argument ( ' - - 
sparams', nargs=3, dest='sparams', required=False, 


help="The three arguments in spara 
ms: R, r, 1.") 


# parse args 


e args = parser.parse args() 


Le DU bp aaadadaeu m --spearamsrp 
орреородоорооопорро 


ΠΠΠΠΠΠΠΠΠΠΠωΤΙΕΠΠΗ 


# set the width of the drawing window to 80 percent of 
the screen widt 


e turtle.setup(widthz0.8) 


4 set the cursor shape to turtle 


о turtle.shape('turtle') 


# set the title to Spirographs! 
e turtle.title("Spirographs!") 

# add the key handler to save our drawings 
ο turtle.onkey(saveDrawing, "s") 

# start listening 


о turtle. 115 Теп () 


Я hide the main turtle cursor 


о turtle.hideturtle() 


О®П005еёир()0800000000008000000000000 
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замеО гамі па 00050000000000Ф0000 
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# check for any arguments sent to -- 
Sparams and draw the Spirograph 


e if args.sparams: 
6 params = [float(x) for x іп args.sparams] 
# draw the Spirograph with the given parameters 
col = (0.0, 0.0, 0.0) 
spiro = Spiro(0, 0, col, *params) 
spiro.draw() 
else: 
# create the animator object 
Ө spiroAnim = SpiroAnimator(4) 
# add a key handler to toggle the turtle cursor 
о turtle.onkey(spiroAnim.toggleTurtles, "t") 
# add a key handler to restart the animation 


ο turtle.onkey(spiroAnim.restart, "space") 


4 start the turtle main Loop 


Ө turtle.mainloop() 
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ег/ spirograph/spiro.py ΠΠΠΠΠΠΠΠΗ 


import sys, random, argparse 





import numpy as np 

import math 

import turtle 

import random 

from PIL import Image 

from datetime import datetime 


from fractions import gcd 


# a class that draws a Spirograph 


class Spiro: 
# constructor 


def init (self, xc, yc, col, К, г, 1): 


# create the turtle object 
self.t = turtle.Turtle() 
# set the cursor shape 


self.t.shape('turtle') 


# set the step in degrees 
self.step = 5 
4 set the drawing complete flag 


self.drawingComplete - False 


4 set the parameters 


self.setparams(xc, yc, col, R, г, 1) 


# initialize the drawing 


self.restart() 


4 set the parameters 
def setparams(self, xc, yc, col, К, г, 1): 


# the Spirograph parameters 


self.xc 


| 
X 
о 


self.yc 


II 
< 
[e] 


self.R = int(R) 
self.r = int(r) 
self.l = l 

self.col = col 


# reduce r/R to its smallest form by dividing with 
the GCD 


gcdVal = gcd(self.r, self.R) 
self.nRot = self.r//gcdVal 
# get ratio of radii 

self.k = r/float(R) 

# set the color 
self.t.color(*col) 

# store the current angle 


self.a = 0 


# restart the drawing 

def restart(self): 
# set the flag 
self.drawingComplete = False 
# show the turtle 


self.t.showturtle() 


# go to the first point 


self.t.up() 

R, k, 1 = self.R, self.k, self.l 

а = 0.0 

x = R*((1-k)*math.cos(a) + l*k*math.cos((1-k)*a/k)) 
y - R*((1-k)*math.sin(a) - l*k*math.sin((1-k)*a/k)) 


self.t.setpos(self.xc + x, self.yc + у) 


self.t.down() 


# draw the whole thing 


def draw(self): 


k)*a/k)) 


k)*a/k)) 


# draw the rest of the points 


R, 


k, 1 = self.R, self.k, self.l 


for i in range(0, 360*self.nRot + 1, self.step): 


a 


math.radians(i) 


x R*((1-k)*math.cos(a) + l*k*math.cos((1- 


y = R*((1-k)*math.sin(a) - l*k*math.sin((1- 


self.t.setpos(self.xc + x, self.yc + y) 
# drawing is now done so hide the turtle cursor 


self.t.hideturtle() 


4 update by опе step 
def update(self): 
# skip the rest of the steps if done 
if self.drawingComplete: 
return 
# increment the angle 
self.a += self.step 
# draw a step 
R, k, 1 = self.R, self.k, self.l 


# set the angle 


a = math.radians(self.a) 

X = self.R*((1-k)*math.cos(a) + l*k*math.cos((1- 
k) *a/k) ) 

у = self.R*((1-k)*math.sin(a) - l*k*math.sin((1- 
K)*a/k)) 


self.t.setpos(self.xc + x, self.yc + y) 
# if drawing is complete, set the flag 
if self.a »- 360*self.nRot: 
self.drawingComplete - True 
# drawing is now done so hide the turtle cursor 


self.t.hideturtle() 


# clear everything 


def clear(self): 


self.t.clear() 


# а class for animating Spirographs 
class SpiroAnimator: 
# constructor 
def init (self, М): 
# set the timer value in milliseconds 
self.deltaT = 10 
# get the window dimensions 
self.width = turtle.window width() 
self.height = turtle.window height () 
# create the Spiro objects 
self.spiros = [] 
for i in range(N): 
# generate random parameters 
rparams = self.genRandomParams() 
# set the spiro parameters 
spiro = Spiro(*rparams) 
self.spiros.append(spiro) 
# call timer 


turtle.ontimer(self.update, self.deltaT) 


# restart spiro drawing 
def restart(self): 
for spiro in self.spiros: 
# clear 
spiro.clear() 
# generate random parameters 
rparams = self.genRandomParams() 
# set the spiro parameters 
spiro.setparams(*rparams) 
# restart drawing 


spiro.restart() 


# generate random parameters 
def genRandomParams (self): 
width, height = self.width, self.height 


R = random.randint(50, min(width, height)//2) 


г = random.randint(10, 9*R//10) 
l 2» random.uniform(0.1, 0.9) 
xc = random.randint(-width//2, width//2) 


yc - random.randint(-height//2, height//2) 
col = (random.random(), 
random.random(), 


random. random() ) 


return (xc, ус, col, В, г, 1) 


def update(self): 
# update all spiros 
nComplete = 0 
for spiro in self.spiros: 
# update 
spiro.update() 
# count completed spiros 
if spiro.drawingComplete: 
nComplete += 1 
4 restart if all spiros are complete 
if nComplete == len(self.spiros): 
self.restart() 
# call the timer 


turtle.ontimer(self.update, self.deltaT) 


# toggle turtle cursor on and off 
def toggleTurtles(self): 
for spiro in self.spiros: 
if spiro.t.isvisible(): 


spiro.t.hideturtle() 


else: 


spiro.t.showturtle() 


# save drawings as PNG files 
def saveDrawing(): 
4 hide the turtle cursor 
turtle.hideturtle() 
# generate unique filenames 
dateStr = (datetime.now()).strftime ( "%0%0%Ү -%SH%M%S" ) 
fileName = 'spiro-' + dateStr 
print('saving drawing to %s.eps/png' % fileName) 
# get the tkinter canvas 
canvas = turtle.getcanvas() 
# save the drawing as a postscipt image 
canvas.postscript(file = fileName + '.eps') 


# use the Pillow module to convert the poscript image f 
ile to PNG 


img = Image.open(fileName + '.eps') 
img.save(fileName + '.png', 'png') 
# show the turtle cursor 


turtle.showturtle() 


# main() function 


def main(): 
# use sys.argv if needed 
print('generating spirograph...') 
# create parser 


descStr = """This program draws Spirographs using the T 
urtle module. 


When run with no arguments, this program draws random S 
pirographs. 


Terminology: 


R: radius of outer circle 
r: radius of inner circle 


1: ratio of hole distance to г 


parser = argparse.ArgumentParser(description=descStr) 


# add expected arguments 


parser.add argument ( - - 
sparams', nargs=3, dest-'sparams', required=False, 


help="The three arguments in sparam 
© Rs Гу 429 


# parse args 


args = parser.parse args() 


# set the width of the drawing window to 80 percent of 
the screen width 


turtle.setup(widthz0.8) 


# set the cursor shape to turtle 


turtle.shape('turtle') 


# set the title to Spirographs! 
turtle.title("Spirographs!") 

# add the key handler to save our drawings 
turtle.onkey(saveDrawing, "s") 

# start listening 


turtle.listen() 


4 hide the main turtle cursor 


turtle.hideturtle() 


# check for any arguments sent to -- 
sparams and draw the Spirograph 


if args.sparams: 


params = [float(x) for x in args.sparams] 


# draw the Spirograph with the given parameters 
col = (0.0, 0.0, 0.0) 
Spiro = Spiro(0, 0, col, *params) 
spiro.draw() 
else: 
# create the animator object 
spiroAnim = SpiroAnimator (4) 
# add a key handler to toggle the turtle cursor 
turtle.onkey(spiroAnim.toggleTurtles, "t") 
# add a key handler to restart the animation 


turtle.onkey(spiroAnim.restart, "space") 


4 start the turtle main loop 


turtle.mainloop() 


# call main 
if name == ' main ': 


main() 
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$ python spiro.py 
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$ python spiro.py --sparams 300 100 0.9 
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# recursive Koch snowflake 

def kochSF(x1, yl, x2, y2, t): 
# compute intermediate points p2, p3 
if segment length > 10: 


# recursively generate child segments 


# flake #1 

kochSF(x1, yl, p1[0], ρ1[1], t) 

# flake #2 

kochSF(p1[0], ρ1[1], p2[0], p2[1], t) 
# flake #3 

kochSF(p2[0], ρ2[1], p3[0], p3[1], t) 
# flake #4 

kochSF(p3[0], p3[1], x2, y2, t) 


else: 


ΠΠΠΠΠΠΠΠΠΠΠΠΠ http://electronut.in/koch- 
snowflake-and-the-thue- morse-sequence/ 
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>>> import numpy as пр 
>>> import matplotlib.pyplot as plt 


>>> import matplotlib.animation as animation 
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Ө >>> 


x = np.array([[0, 0, 255], [255, 255, 0], І0, 255, 011) 
Ө >>> plt.imshow(x, interpolation='nearest' ) 


plt.show() 
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np.random.choice([0, 255], 4*4, р=[0.1, 0.9]).reshape(4, 4) 


000000 
array([[255, 255, 255, 2551, 
|255; 255, 255, 255], 
[255, 255, 255, 255], 


[255, 255, 255, 0]]) 
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def addGlider(i, 1, grid): 
"""adds a glider with top left cell at (i, j)""" 
Ө glider = np.array([[0, 0, 255], 
|255; 0, 255]; 
[0, 255, 255]]) 
Ө grid[i:i43, j:j+3] = glider 
© grid = np.zeros(N*N).reshape(N, N) 


Ө addGlider(1, 1, grid) 
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right = ΜΗ 
m = grid[i][j-*1] 
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>>> N = 16 
>>> 1l = 14 
>>> 12 = 15 


>>> (11141) ΕΝ 


>>> (1241) ΕΝ 
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Обрр0Рерооорооромородрородроорроооообо0 
ООО00000000000 


right = grid[i][(j-*1)9N] 
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# apply Conway's rules 
if grid[i, j] == ON: 
e if (total < 2) or (total > 3): 
newGrid[i, j] = OFF 
else: 
if total == 3: 


e newGrid[i, j] = ON 
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ПОМПО000000МП 
ООО00000000000 
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Я main() function 


def main(): 
# command line argumentss are in sys.argv[1], sys.argv[2], 


# sys.argv[0] is the script name and can be ignored 
@ # parse arguments 
parser = argparse.ArgumentParser(description="Runs Conway's 
Game of Life 
Simulation.") 
# add arguments 


Ө parser.add argument ('--grid- 
size', dest='N', required=False) 


Ө parser.add argument ('--mov- 
file', dest='movfile', required=False) 


Ө parser.add argument('-- 
interval', dest='interval', required=False) 


Ө parser.add argument ( ' - - 
glider', action='store true', required=False) 


args = parser.parse args() 
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# set grid size 
N = 100 
if args.N and int(args.N) > 8: 


N = int(args.N) 


# set animation update interval 
updateInterval = 50 
if args.interval: 


updateInterval = int(args.interval) 


4 declare grid 
Ө grid = np.array([]) 
# check if "glider" demo flag is specified 
if args.glider: 
grid = np.zeros(N*N).reshape(N, М) 


addGlider(1, 1, grid) 


else: 
# populate grid with random on/off - more off than on 


grid = randomGrid(N) 
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4 set up the animation 
Ө fig, ax = plt.subplots() 
img = ax.imshow(grid, interpolation='nearest' ) 


Ө ani = animation.FuncAnimation(fig, update, fargs= 
(img, grid, М, ), 


frames=10, 
interval=updatelnterval, 


save count=50) 


# number of frames? 
# set the output file 
if args.movfile: 


ani.save(args.movfile, fps=30, extra args-['- 
vcodec', 'libx264']) 


plt.show() 
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глазСег /сопмау/сопмау.РУППООООООО 
import sys, argparse 
import numpy as np 
import matplotlib.pyplot as plt 


import matplotlib.animation as animation 


ON = 255 
OFF = 0 


vals = [ON, OFF] 


def randomGrid(N): 
"""returns a grid of NxN random values""" 


return np.random.choice(vals, N*N, p- 
[0.2, 0.8]).reshape(N, N) 


def addGlider(i, j, grid): 


"""adds a glider with top-left cell at (i, j)""" 


glider = пр.аггау([[0, 0, 255], 
[255, 0, 255], 
[0, 255, 255]]) 


grid[i:i*3, j:j+3] = glider 


def update(frameNum, img, grid, N): 
# copy grid since we require 8 neighbors for calculation 
# and we go line by line 
newGrid = grid.copy() 
for i in range(N): 
for j in range(N): 


# compute 8- 
neghbor sum using toroidal boundary conditions 


# x and y wrap around so that the simulation 

# takes place on a toroidal surface 

total = int((grid[i, (j-1)%N] + grid[i, (j+1)%N] + 
grid[(i-1)%N, j] + grid[(i+1)%N, j] + 


grid[(i-1)9N, (j-1)%N] + grid[(i- 
1)%N, (j+1)%N] + 


grid[(i+1)%N, (j- 
1)%N] + grid[(i+1)%N, (j+1)%N])/255) 


# apply Conway's rules 
if grid[i, j] == ON: 


if (total < 2) or (total > 3): 


newGrid[i, j] = OFF 
else: 
if total == 3: 
newGrid[i, j] = ON 


# update data 
img.set data(newGrid) 
grid[:] = newGrid[:] 


return img, 


4 main() function 


def main(): 
# command line arguments are in sys.argv[1], sys.argv[2], 


# sys.argv[0] is the script name and can be ignored 
# parse arguments 
parser = argparse.ArgumentParser(description="Runs Conway's 
Game of Life 
simulation.") 
# add arguments 


parser.add argument('--grid- 
size', dest='N', required=False) 


parser.add argument ( ' - -mov- 
file', dest='movfile', required=False) 


parser.add argument('-- 
interval', dest='interval', required=False) 


parser.add argument ( ' - - 
glider', action='store_true', required=False) 


parser.add argument ( ' - - 
gosper', action='store true', required=False) 


args = parser.parse args() 


# set grid size 
N = 100 
if args.N and int(args.N) > 8: 


N = int(args.N) 


# set animation update interval 
updateInterval = 50 
if args.interval: 


updateInterval = int(args.interval) 


# declare grid 
grid = np.array([]) 
# check if "glider" demo flag is specified 
if args.glider: 
grid = np.zeros(N*N).reshape(N, N) 


addGlider(1, 1, grid) 


else: 
# populate grid with random on/off - more off than on 


grid = randomGrid(N) 


# set up the animation 
fig, ax = plt.subplots() 
img = ax.imshow(grid, interpolation='nearest' ) 


ani = animation.FuncAnimation(fig, update, fargs= 
(img, grid, N, ), 


frames-10, 
interval=updatelnterval, 
save count=50) 

# number of frames? 

# set the output file 

if args.movfile: 


ani.save(args.movfile, fps=30, extra args-['- 
vcodec', 'libx264']) 


plt.show() 
4 call main 
if name == ' main ': 


main() 
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A = sin (2n ft) 


ΠΕ : ΠΠΠΠΠΠΠ #0000 £ О0000000000000000 


А = sin (2n П/ R) 
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import numpy as np 


import wave, math 


sRate = 44100 
nSamples = sRate * 5 
Ө x = np.arange(nSamples)/float(sRate) 
Ө vals = np.sin(2.0*math.pi*220.0*x) 
Ө data = np.array(vals*32767, 'intl6').tostring() 


file = wave.open('sine220.wav', 'wb') 


Ө file.setparams((1, 2, sRate, nSamples, 'NONE', 'uncompres 
sed')) 


file.writeframes(data) 


file.close() 
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0000000000Кагріц5-5ігопо ПОДОДОПОД00000 
ОДООО0000000000000000000000004.400 
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ПППППППППКагрім5-5ігопд ПО0О00000000000 
ОДОДОООПРУЄпоправФацерроооб"авФск "00 
Pythonllcollections] ОДООО000000000000000 
О0000000000000000009еачеПо0000000000 
ΠΠΠΠΠΠΠΠΠΠΑ͂-5ΠΠΠΠΠΠΠΠΠΠΠΠΟ(1})Π΄ΠΠΠΗΗ 
О00000000аеачероо00000000000000000 
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04-5 ΠΠαθαιεΠΠΠΠΠΗ 
ППЦПППШППШПППРУ һопЦПЦПдедие[ 


>>> from collections import deque 
Ө >>> d = deque(range(10)) 

>>> print(d) 

deque([9,- 1, 2, 3, 4, 5, 6, 7, 8, 9]) 
Ө >>> d.append(-1) 


>>> print(d) 


deque([0, 1, 2, 3, 4, 5, 6, 7, 8, 9, -1]) 
© >>> d.popleft() 

0 

>>> print d 


deque([1, 2, 3, 4, 5, 6, 7, 8, 9, -1]) 
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О000аеачеПВ00000000000Кагріиѕ-5ёгопоП 
ΗΗΠΗΠΠΗ 


я generate note of given frequency 
def generateNote( freq): 
nSamples = 44100 
sampleRate = 44100 
М = int(sampleRate/freq) 
4 initialize ring buffer 
€ buf = deque([random.random() - 0.5 for i in range(N)]) 
4 initialize samples buffer 
Ө samples = np.array([0]*nSamples, 'float32') 
for i in range(nSamples): 


e samples[i] = buf[0] 


о ауд = 0.996*0.5*(buf[0] + buf[1]) 
buf .append (avg) 


buf.popleft() 


# convert samples to 16-bit values and then to a string 
# the maximum value is 32767 for 16-bit 
e samples = np.array(samples*32767, 'іпі16') 


о return samples.tostring() 
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def writeWAVE(fname, data): 
# open file 
Ө file = wave.open(fname, 'wb') 


# WAV file parameters 


nChannels = 1 
sampleWidth = 2 
frameRate = 44100 
nFrames = 44100 


# set parameters 


file.setparams((nChannels, sampleWidth, frameRate, nFrames, 
'NONE', 'noncompressed')) 
Ө file.writeframes(data) 


file.close() 


Le DO wAVDBBTBeUmn p b d p © 
орробрродоооореФвроророо 
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DOLI Python pygametn LILIWAYVII II 
рудаплербПоОДООДООПРуєпопОООООД0000000 
ο πο μμ.” 


О0000000000003оѓеРІауег0000000 


def init (self): 


e pygame.mixer.pre init(44100, -16, 1, 2048) 
pygame.init() 
4 dictionary of notes 
ө self.notes = {} 
# add a note 
def add(self, fileName): 
© self.notes[fileName] = pygame.mixer.Sound(fileName) 
# play a note 
def play(self, fileName): 
try: 
о self.notes[fileName].play() 
except: 
print(fileName + ' not found!') 
def playRandom(self): 
"""play a random note""" 
e index = random.randint(0, len(self.notes) -1) 
о note = list(self.notes.values()) [index] 


note.play() 
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О00000000 


4.3.5 main()[ 
ПО00таіп()О000000000000000000000000 


parser = argparse.ArgumentParser(description="Generating so 
unds wi 


Karplus String Algorithm") 
# add arguments 


Ө parser.add argument('-- 
display', action='store true’, required=False) 


parser.add _argument ( ' - 
play', action='store | п” required=False) 


parser.add argument ( ' - 
piano', ο core πα required=False) 


args = parser.parse args() 


# show plot if flag set 
if args.display: 
gShowPlot = True 


plt.ion() 


# create note player 


nplayer = МоїеР1ауег() 


print('creating notes...') 
for name, freq in list(pmNotes.items()): 
fileName = name + '.wav' 
if not os.path.exists(fileName) or args.display: 
data = generateNote(freq) 
print('creating ' + fileName + '...') 
writeWAVE(fileName, data) 
else: 


print('fileName already created. skipping...') 


# add note to player 


nplayer.add(name + '.wav') 


4 play note if display flag set 
if args.display: 
nplayer.play(name + '.wav') 


time.sleep(0.5) 


# play a random tune 


if args.play: 


while True: 
try: 
e nplayer.playRandom() 
# rest - 1 to 8 beats 
о rest = np.random.choice([1, 2, 4, 8], 1, 
р=[0.15, 0.7, 0.1, 0.051) 
time.sleep(0.25*rest[0]) 
except KeyboardInterrupt: 


exit() 
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ОООДОДО0000000000000000 
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рр/біоб/плазіег/кагріи5/к5.рУПО 
ітрогі 5у5, 05 
import time, random 
import wave, argparse, pygame 
import numpy as np 
from collections import deque 


from matplotlib import pyplot as plt 


# show plot of algorithm in action? 


gShowPlot = False 


# notes of a Pentatonic Minor scale 
Я piano C4-E(b)-F-G-B(b)-C5 


pmNotes = {'C4': 262, 'Eb': 311, 'F': 349, 'G':391, 'Bb':46 
6} 


# write out WAV file 


def writeWAVE(fname, data): 


# open file 

file = wave.open(fname, 'wb') 
# WAV file parameters 
nChannels - 1 

ѕатр1еміа+һ = 2 

frameRate = 44100 

nFrames = 44100 


# set parameters 


file.setparams((nChannels, sampleWidth, frameRate, nFrames, 
'NONE', 'noncompressed')) 
file.writeframes(data) 


file.close() 


4 generate note of given frequency 
def generateNote(freq): 
nSamples = 44100 
sampleRate - 44100 
М = int(sampleRate/freq) 
8 initialize ring buffer 
buf = deque([random.random() - 0.5 for i in range(N)]) 
# plot of flag set 


if gShowPlot: 


axline, = plt.plot(buf) 
# initialize samples buffer 
samples = np.array([0]*nSamples, 'float32') 
for i in range(nSamples): 

samples[i] = buf[0] 

avg = 0.995*0.5*(buf[0] + buf[1]) 

buf.append(avg) 

buf.popleft() 

# plot of flag set 

if gShowPlot: 

if i % 1000 == 0: 
axline.set ydata(buf) 


plt.draw() 


# convert samples to 16-bit values and then to a string 
# the maximum value is 32767 for 16-bit 
samples = np.array(samples*32767, 'intl6') 
return samples.tostring() 
# play a WAV file 
class NotePlayer: 
# constructor 
def init (self): 


pygame.mixer.pre init(44100, -16, 1, 2048) 


pygame.init() 
4 dictionary of notes 
self.notes = {} 
# add a note 
def add(self, fileName): 
self.notes[fileName] = pygame.mixer.Sound (fileName) 
# play a note 
def play(self, fileName): 
try: 
self.notes[ fileName] .play() 
except: 
print(fileName + ' not found!') 
def playRandom(self): 
"""play a random note""" 
index = random.randint(0, len(self.notes) -1) 
note = list(self.notes.values()) [index] 


note.play() 


4 main() function 
def main(): 
4 declare global var 


global gShowPlot 


parser = argparse.ArgumentParser(description="Generating so 
unds with 


Karplus String Algorithm") 
# add arguments 


parser.add argument('-- 
display', action-'store true', required-False) 


parser.add argument('-- 
play', action-'store true', required-False) 


parser.add argument('-- 
piano', action-'store true', required-False) 


args = parser.parse агд5() 


# show plot if flag set 
if args.display: 
gShowPlot - True 


plt.ion() 


# create note player 


nplayer = NotePlayer() 


print('creating notes...') 
for name, freq in list(pmNotes.items()): 


fileName = name + '.wav' 


if not os.path.exists(fileName) or args.display: 


data = generateNote(freq) 


print('creating ' + fileName + '...') 
writeWAVE(fileName, data) 


else: 


print('fileName already created. skipping 


# add note to player 


nplayer.add(name + '.wav') 


4 play note if display flag set 
if args.display: 
nplayer.play(name + '.wav') 


time.sleep(0.5) 


# play a random tune 
if args.play: 
while True: 


try: 


nplayer.playRandom() 


# rest - 1 to 8 beats 


rest = np.random.choice([1, 2, 4, 8], 1, 


p=[0.15, 0.7, 0.1, 0.05]) 


time.sleep(0.25*rest[0]) 
except KeyboardInterrupt: 


exit() 


# random piano mode 
if args.piano: 
while True: 
for event in pygame.event.get(): 
if (event.type == pygame.KEYUP) : 
print("key pressed") 
nplayer.playRandom() 


time.sleep(0.5) 


# call main 
if name == ' main ': 


main() 


4.5 ΠΠ 
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$ python3 ks.py —display 
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5 python ks.py -ріау 
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ООО000000000000 
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ООО0О0000000000000000000000000000000 
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г ПОООООООО 
г ПОООООООО 
г ПОДОД0000 

• ПОО0000000000 

• ДО00000000 


ОООДО000000000000000000000000000 


5.2 [000 
ОБОДООООПОРУЄАЛоГОДО 


: пипоруД0000000000000000 
- ΠΙΔΙΡΙΟΙΙΙΟΠΠΠΠΠΠΠΠΠΠΠΗ 
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5.3.1 QQUUUUUUUUUU 


О0000000п um руДбО0ОО0ОО0О00000000000000 
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Ө import math 


O import numpy as np 
© width, height = 640, 480 


@ pos = ος 0, height/2.0] + 10*np.random.rand(2*N). 
shape(N, 


Ө angles = 2*math.pi*np. random. rand(N) 


Ө vel = np.array(list(zip(np.sin(angles), np.cos(angles)))) 
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>>> zip([®, 1, 21, ІЗ, 4, 51) 


[(0, 3), (1, 4), (2, 5)] 
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def applyBC(self): 
"""apply boundary conditions""" 
deltaR = 2.0 
for coord in self.pos: 
e if coord[0] > width + deltaR: 
coord[0] = - deltaR 
if coord[0] < - deltaR: 
coord[0] = width + deltaR 
1f coord[1] > height + deltaR: 
coord[1] = - deltaR 
if coord[1] < - deltaR: 


coord[1] = height + deltaR 


ОӨПОООхОО00000000000000000000000000 
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fig = plt.figure() 


ax = plt.axes(xlim=(0, width), ylim=(0, height) ) 


9 
pts, = ax.plot([], [], markersize=10, c='k', marker='o', 1 
s='None' ) 


beak, = ax.plot([], [], markersize=4, c='r', marker='o', 1 
s-'None') 


anim - animation.FuncAnimation(fig, tick, fargs- 
(pts, beak, boids), 


interval=50) 


Lee n iba d ets mam beakr in 
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000000 


Ө vec = self.pos + 10*self.vel/self.maxVel 


Ө beak.set data(vec. reshape(2*self.N) 
[::2], vec.reshape(2*self.N)[1::2]) 


педородобрродорробометоооробо01080000000 
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[]reshape[]matplotlib[][][|set даѓаП01::2 1000 
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5.3.4 ПОООООО 


ООДООРУСпоПОДОООООЗООО0000"питруро 
ОР ОДОДОД00000000000000пепору000 


import питру as np 


from scipy.spatial.distance import squareform, pdist, cdist 


def test2(pos, radius): 
# get distance matrix 
e distMatrix = squareform(pdist(pos)) 
# apply threshold 
e D = distMatrix < radius 
# compute velocity 
e vel = pos*D.sum(axis=1).reshape(N, 1) - D.dot(pos) 


return vel 


Πθ[[]ή5αιµαγεί[οΓπη()Πραίξι()ΠΠΠΠδείρυ[Π[Πῃ 
000000000000000000000000000000000000 
ο ο μμ - 


>>> import numpy as np 
>>> from scipy.spatial.distance import squareform, pdist 
>>> x = np.array([[0.0, 0.0], [1.0, 1.0], [2.0, 2.611) 


>>> squareform(pdist(x) ) 


аггау([[ 0. , 1.41421356, 2.82842712], 
[ 1.41421356, 0. , 1.41421356], 


[ 2.82842712, 1.41421356, 0. ]]) 


5дчагеїогтп (0000003 x 3000000 ij ООООР iL 
P УДООООООО0000Ф6000000000000000000300 
ΕΕ ΕΕ 


>>> squareform(pdist(x)) < 1.4 
array([[ True, False, False], 
[False, True, False], 


[False, False, True]], dtype=bool) 


< III pp DO DEO 8 0000000000О 
ОООО0000000000000000000000000000 
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теі ООО000000000000000000РуһопП00 
DOLIO testi ПсеѕЕ2П000000000еѕ.рур000 


>>> [from timeit import timeit 


>>> 
timeit('testl(pos, 100)", 'from test import testl, М, pos, 


width, height', 

number=100) 

7.880876064300537 

>> 

timeit(' μι 100)', 'Тгот test import test2, М, pos, 
width, height' 

number=100) 


0.0369696 
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def applyRules(self): 
# apply rule #1: Separation 
D = distMatrix < 25.0 
e 
УТИ ИЕ 1) - D.dot(sel 


e self.limit(vel, self.maxRuleVel) 


4 distance threshold for alignment (different from separati 


D = distMatrix < 50.0 


# apply rule #2: Alignment 
© vel2 = D.dot(self.vel) 
self. limit(vel2, self.maxRuleVel) 


vel += vel2; 


# apply rule #3: Cohesion 
о vel3 = D.dot(self.pos) - self.pos 
self. limit(vel3, self.maxRuleVel) 


vel += vel3 


return vel 
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ООО000000000000000000000000000000000 
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# add a "button press" event handler 


cid = fig.canvas.mpl connect('button press event', buttonPr 
ess 


перпотрі ςσοππθςϊ(}ΠΠΗΠΠΙ8ϊρ|ο{1ΠΘΠΠΠΠΠΠΠ 
που Рге55 0000000 


ОООО0000000000000000000000000 


def buttonPress(self, event): 
"""event handler for matplotlib button presses""" 
# left-click to add a boid 

Ө if event.button is 1: 


о self.pos = np.concatenate((self.pos, 


np.array([[event.xdata, event.ydata]])), 


axis=0) 


4 generate a random velocity 


© angles = 2*math.pi*np. random. гапа (1) 


v = np.array(list(zip(np.sin(angles), np.cos(angles)))) 
self.vel = np.concatenate((self.vel, v), axis=0) 


self.N += 1 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘΠΠΠ(Ενεπῖ.χα818Π 
емепі.удаѓа)П0000000000000000000ө0000 
πο μου 
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# right-click to scatter boids 
e elif event.button is 3: 
# add scattering velocity 


self.vel += 0.13 
(self.pos - np.array([[event.xdata, event.ydata]])) 


ПФПППИПИПИППИПИПИПИППӨППППИПИПЯПГГІГ 
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ОО000000000030000000000000000000 
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ООО000000000000000 


9 
parser = argparse.ArgumentParser(description="Implementing 
Craig 

Reynolds's Boids...") 


# add arguments 


parser.add argument ( ' - -num- 
boids', dest='N', required-False) 


args - parser.parse args() 


Я set the initial number of boids 
N = 100 
if args.N: 


М = int(args.N) 


4 create boids 


boids = Boids(N) 


таіп()00000®0000000000000000000 
argparse[|[|[] 


5.3.8 Boids|] 


ΠΠΠΠΠΒο!45ΠΠΠΠΠΠΠΠΠ 


class Boids: 
"""Class that represents Boids simulation""" 
def init (self, М): 
"""initialize the Boid simulation""" 
# initial position and velocities 
о 
self.pos = [width/2.0, height/2.0] + 10*np.random.rand(2*N) 
.reshape(N, 2) 
# normalized random velocities 
angles - 2*math.pi*np.random.rand(N) 
LE - np.array(list(zip(np.sin(angles), np.cos(angles) 
self.N = М 
# minimum distance of approach 


self.minDist = 25.0 


# maximum magnitude of velocities calculated by "rules" 
self.maxRuleVel = 0.03 
# maximum magnitude of the final velocity 


self.maxVel = 2.0 


BoidD pp ppp mp p d eri mm pr D 
ΕΠΗ 


роіаѕ.їск()008000000000000000000000 


def tick(frameNum, pts, beak, boids): 
#print frameNum 
"""Урдаїе function for animation""" 
boids.tick(frameNum, pts, beak) 


return pts, beak 


ООО000000000000000000000000000000000 
00000000 


def limitVec(self, vec, maxVal): 
"""limit the magnitude of the 2D vector""" 
mag = norm(vec) 


if mag » maxVal: 


vec[0], vec[1] = vec[0]*maxVal/mag, vec[1]*maxVal/mag 


e def limit(self, X, maxVal): 


"""limit the magnitude of 2D vectors in array X to 
maxValue" пп 


тог мес 1п Х: 


self.limitVec(vec, maxVal) 


ΠΘΓΊΠΠΠΙΙΠΙΙΕ)ΠΠΠΠΠΠΠΠΠΠΠΠΗΠΠΠΠΠΠΠΠΠΠΗ 


5.4 [| 
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https://github.com/electronut/pp/blob/ 


master/boids/boids.pyQ000000 


import 
import 
import 
import 


import 


from scipy.spatial.distance import squareform, pdist, cdist 


Sys, argparse 

math 

numpy as np 
matplotlib.pyplot as plt 


matplotlib.animation as animation 


from numpy.linalg import norm 


width, 


height = 640, 480 


class Boids: 


"""Class that represents Boids simulation""" 


def init (self, N): 


self.pos = [width/2.0, height/2.0] + 10*np.random.rand(2*N) 


"""initialize the Boid simulation""" 


4 initial position and velocities 


.reshape(N, 2) 


4 normalized random velocities 


angles 2*math.pi*np. random. гапа (№) 


self.vel = np.array(list(zip(np.sin(angles), np.cos(angles) 
))) 

self.N = М 

# minimum distance of approach 


self.minDist = 25.0 


# maximum magnitude of velocities calculated by "rules" 
self.maxRuleVel = 0.03 
# maximum maginitude of the final velocity 


self.maxVel = 2.0 


def tick(self, frameNum, pts, beak): 
"""Update the simulation by one time step.""" 
# get pairwise distances 
self.distMatrix = squareform(pdist(self.pos)) 
# apply rules: 
self.vel += self.applyRules() 
self.limit(self.vel, self.maxVel) 
self.pos += self.vel 
self.applyBC() 
# update data 


pts.set data(self.pos.reshape(2*self.N)[::2], 


self.pos.reshape(2*self.N)[1::2]) 
vec = self.pos + 10*self.vel/self.maxVel 
beak.set data(vec.reshape(2*self.N)[::2], 


мес. reshape(2*self.N)[1::2]) 


def limitVec(self, vec, maxVal): 
"""limit the magnitide of the 2D vector""" 
mag = norm(vec) 


if mag » maxVal: 


vec[0], vec[1] = vec[0]*maxVal/mag, vec[1]*maxVal/mag 


def limit(self, X, maxVal): 


"""limit the magnitide of 2D vectors in array X to ma 
xValue" пп 


тог мес in Х: 


self.limitVec(vec, maxVal) 


def applyBC(self): 
"""apply boundary conditions""" 
deltaR = 2.0 
for coord in self.pos: 


if coord[0] > width + deltaR: 


coord[0] = - deltaR 

if coord[0] < - deltaR: 
coord[0] = width + deltaR 

if coord[1] > height + deltaR: 
coord[1] = - deltaR 

if coord[1] < - deltaR: 


coord[1] = height + deltaR 


def applyRules(self): 


# apply rule #1: Separation 


D = self.distMatrix < 25.0 


vel = self.pos*D.sum(axis=1).reshape(self.N, 1) - D.dot(sel 
f.pos) 


self.limit(vel, self.maxRuleVel) 


# distance threshold for alignment (different from separati 
on) 


D = self.distMatrix < 50.0 


# apply rule #2: Alignment 
vel2 = D.dot(self.vel) 
self.limit(vel2, self.maxRuleVel) 


vel += vel2; 


4 apply rule #3: Cohesion 
vel3 = D.dot(self.pos) - self.pos 


self. limit(vel3, self.maxRuleVel) 


vel += vel3 


return vel 


def buttonPress(self, event): 


"""avent handler for matplotlib button presses""" 
# left-click to add a boid 


if event.button is 1: 


self.pos = np.concatenate((self.pos, 


np.array([[event.xdata, event.ydata]])), 
axis=0) 
# generate a random velocity 


angles = 2*math.pi*np. random. гапа (1) 


у = np.array(list(zip(np.sin(angles), np.cos(angles) ))) 


self.vel = np.concatenate((self.vel, v), axis=0) 


self.N += 1 


# right-click to scatter boids 


elif event.button is 3: 
# add scattering velocity 


self.vel += 0.1* 
(self.pos - np.array([[event.xdata, event.ydata]])) 


def tick(frameNum, pts, beak, boids): 
#print frameNum 
"""update function for animation""" 
boids.tick(frameNum, pts, beak) 


return pts, beak 


4 main() function 
def main(): 

# use sys.argv if needed 

print('starting boids...') 
parser = argparse.ArgumentParser(description-"Implementing 
Craig 

Reynold's Boids...") 
# add arguments 


parser.add argument('--num- 
boids', dest='N', required=False) 


args = parser.parse args() 


# set the initial number of boids 


М = 100 
if args.N: 


N = int(args.N) 


# create boids 


boids = Boids(N) 


# set up plot 
fig = plt.figure() 


ax = plt.axes(xlim=(0, width), ylim=(0, height) ) 


pts, = ax.plot([], ІІ, markersize=10, c='k', marker='o', ls 
-'None') 
реак, = ax.plot([], ГІ, markersize-4, c='r', marker='o', ls 
-'None') 


anim - animation.FuncAnimation(fig, tick, fargs- 
(pts, beak, boids), 


interval=50) 


# add a "button press" event handler 


cid = fig.canvas.mpl connect('button press event', boids.bu 
ttonPress) 


plt.show() 


# call main 
if _ name == ' main ': 


main() 


5.5 ΗΠ 
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$ python3 boids.py 
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06-1 ΑΞΟΙΙΠΠΗΠΠΠΜΠΗ 


6.1 ШШ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 
00000006-200000000000000000000000000 
и 


А5ЗСПОДОООООООООО0О00О0000000А5СПО00000 
ОДДОДАСВвОропорбобрб000000000000000 
О000А5СИПООООО" 00” О00000000000000000 





06-2 ΠΠΠΠΠΠΠΠ 


ОО000000000000800000000000000000000 
[0,255 JO8 000000000008 0000000000000 
255 П0000000000000000 


О000000000М x МОПОДОДООООДОМОМОАЗСІО 
ОООООООООООООооОооооооооооооооо00000 
ΠΠΠΠΠΠΠΠΑΘΕΟΙΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΙΟ,2551ΠΠ 
ППИПИППППАБСИЦППИПИПИППИПИПИППГИГПИП 


ΠΠΠΑ5ΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ(ουτίετ 
ООООО0000000000000000000000000000000 
ООО0000000000000000 


О0000“000“О0000000000000000000000000 
ОО00000000000000000000000000А5С1000 
ОО0000000000000000А5СИО0000000000000 
ОООООООООООООооОооооооооооо000000000 
ОООО0000000000000000 


О0000000000000000000000Соигіег000000 
ОООДО0000000000000000000000 


О000000000А5С00000000 

100000000000 

21000007 x МОООЦ 
ЗОДОМОДОДОООО00000000000 
А0ООО0О00000000000000000000000А5СП000 
ЗООООПА5СФСПОДОДООООООДОД00000000 


6.2 (ТІП 


ΠΠΠΠΠΠΠΡΙ!ΠοννΠΡγ{ΠπΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ООО00000000000000000пипруд 0000000 


6.3 || 


О0000000000000А5СІО0000000000000000 
ОО000000000000000000000000А5СІП0000 
ОДОО000000000000000000000000000000000 
О0000000 


О00000000006.400 


6.3.1 ΠΠΠΠΠΠΠΠΠ 


ОООО000000000000000000000000000000 
ASCIIQUU 


# 70 levels of gray 


Ө gscalel = 0. анна, 
()1{}[]?-_+~<>i!lI;:,\"^ 


# 10 levels of gray 


Ө gscale2 = "@%#*+=-:. " 


еппрозсат'е10700000000ФПравсаїтегрДрпП 
1О0ДО000000000000000000000000000000 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΡϑΗΙ Bourkel 
[Character Representation of Grey Scale 
Ітадеѕ0000 
http://paulbourke.net/dataformats/asciiart/ 


ПП 
ООООО00000000000000000000000000000 


# open the image and convert to grayscale 
Ө image = Image.open( fileName) .convert("L") 
# store the image dimensions 
Ө W, Н = image.size[0], image.size[l] 
# compute the tile width 
Ө w = W/cols 
# compute the tile height based on the aspect ratio and sca 
le of the font 
Ө h = w/scale 
# compute the number of rows to use in the final grid 


Ө rows = int(H/h) 


[Je [ηηιασς.ορεη()/ΠΠΠΠΠΠΠΠΗ 
Image.convert( 000000000000“ e00 
luminance[|||[|[ |[ IILI IDD] 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘΠΠΠΠΠΠΠΠΠΠΠ 
Осої5 0 ОДОО00000000000000000000000000 
ш кши 


ΠΠΠΠΠΠΠΠΠΠΘΠΠΠΠΠΠΠΠΠΠΠΠ5ς8!1ΕΠΠΠΠΠΠΠΠ 
ΠΠΠΠΘΠΠΠΠΠΠΠΠΠΠΠΗΠΠΠ 


ООО000000000000000000000000000000000 
ΠΠΠΠΠΠ5ς8!ΕΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΟ.43ΠΠΠ 


Соигіег 000000000 

6.3.2 ПШПШ 
О00000000000000000000009еёАхегадец()0 
ΕΠΗ 


Ө def getAverageL (image): 

# get the image as a numpy array 
O im - np.array(image) 

# get the dimensions 
© w,h = im.shape 

# get the average 


Ө return np.average(im. reshape (w*h) ) 


ПӨППППППППР!. тадеППППППӨППП!тадеП 
ΠΠΠΠπαπΙΡΥ/ΠΠΓΠΠΠΙΠΙΠΗΠΠΠΗΠΠΗΠΠΗΠΠΗΠΠΠΠΗ 
Le rp m bab p OO Dd e mr 
питру.амегаде( ОДОДОДОДОДОДОДО00 
питру.геѕћаре()058000000(%0һ)000000000 
ОДО000000000000мР-8)000 
питру.амегаде( ОДОДОДОДОДОДОДОДО 


6.3.3 ΠΠΠΠΠΑΞΕΟΗΠΠ 
О0000000000000А5С1000 


# an ASCII image is а list of character strings 


Ө aimg = [] 
# generate the List of tile dimensions 
for j in range(rows): 
e yl - int(j*h) 
y2 = int((j+1)*h) 
# correct the last tile 
if j == rows-1: 
y2 =H 
# append an empty string 
© aimg.append("") 
for i in range(cols): 
# crop the image to fit the tile 
о х1 = int(i*w) 
х2 = int((i-1)*w) 
# correct the last tile 
e if i == С015-1: 
x2 = W 
# crop the image to extract the tile into another Image obj 
ect 
о img = ітаде.сгор((х1, yl, x2, y2)) 
# get the average luminance 


© avg = int(getAverageL (img) ) 


# look up the ASCII character for grayscale value (avg) 
if moreLevels: 

Ө gsval = gscalel[int((avg*69)/255)] 
else: 

© gsval = gscale2[int((avg*9)/255)] 
# append the ASCII character to the string 


® aimg[j] += gsval 


ППППИПИППАБСИЦППИПИППИПИПИПИППППФГІІ 
О00000000000000000000000өП000000000 
ООООО0000000УПОЗООО00000000000000000 
ΕΕ ΕΕ ΕΕ 


ООО000000000000000000000000000000000 
ОООДОО0О0000000000000У0 00У0000000000 
ООО000000000000000000 


Le inni Sc nip pd m p d dp Op Od DULL 
ОООО0000000000000000000000000 


прерроросододродосохордоевроророробох 
О0000000У0000000®П000ітаде.сгор()00000 
00000000009еѓАуегадец)00®П0100006.3.2П 
О000000000000006000000000100255 1000 
Г0091000100000000000000009ѕсаіе200000 
ОО000000000000000А5СИПОӨПО0О00000000 


00000000007 000000000000000®0000000 
О00000А5СПоѕма00000000000000 


6.3.4 ППП 
О000000000000000000000000Оагарат5еп) 


parser = argparse.ArgumentParser(description="descStr" ) 
# add expected arguments 


Ө parser.add argument ( ' - 
file', dest-'imgFile', ο nus) 


O  parser. add 2. 
scale', dest='scale' гест аа- False) 


Ө  parser.add argument ( - 
out', dest-'outFile', геи ος 


Ө  parser.add argument(' -- 
cols', dest='cols', required=False) 


Ө parser.add argument (' -- 
morelevels', dest='moreLevels', action='store true') 


перроровододбродоворбоорододовсоооо00 
орреодробороррерроАасенпдорро0оророесо 
[L]--morelevelsT ППОДООООД0000000000 


6.3.5 [ΠΙΑΞΕΟΙΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ОДОДОДООДА5СФНООДООО0000000000 


4 open a new text file 


Ө f = open(outFile, 'w') 


# write each string in the list to the new file 
Ө for row in aimg: 
f.write(row + '\n') 
# clean up 


Ө f.close() 


ΠΘΠΠΠΠΠΠΠοΡ8π{}ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ο ыш ШШ 


6.4 [II 


ΠΠΗΠΠΠΑΣΟΙΙΠΗΗΠΗΠΗΠΠΠΗ 
https://github.com/electronut/pp/ 


blob/master/ascii/ascii.py[]L]LILILIL]LILILI 
import sys, random, argparse 
import numpy as np 
import math 


from PIL import Image 


4 grayscale level values from: 


# http://paulbourke.net/dataformats/asciiart/ 


# 70 levels of gray 


gscalel = "$@B<8&WM#* oahkbdpqwmZO00QLCIUYXzcvunxrj ft/N | ()1{} 
[]?-_+~<>i!1lI;:,\"% 


# 10 levels of gray 


gscale2 = '@%#*+=-:. 


def getAverageL (image): 


Given PIL Image, return average value of grayscale value 
# get image as numpy array 

im - np.array(image) 

# get the dimensions 

w,h = im.shape 

# get the average 


return np.average(im.reshape(w*h)) 


def covertImageToAscii(fileName, cols, scale, moreLevels): 


Given Image and dimensions (rows, cols), returns an m*n lis 
t of Images 


4 declare globals 


global gscalel, gscale2 


# open image and convert to grayscale 

image = Image.open(fileName).convert('L') 

# store the image dimensions 

W, H = image.size[0], image.size[l] 

print("input image dims: %d x πα" % (W, Н)) 

4 compute tile width 

w = W/cols 
# compute tile height based on the aspect ratio and scale o 
f the font 

h - w/scale 

# compute number of rows to use in the final grid 


rows = int(H/h) 


print("cols: %d, rows: %d" % (cols, rows)) 


print("tile dims: %d x %d" % (м, ћ)) 


# check if image size is too small 
if cols > W or rows > Н: 
print("Image too small for specified cols!") 


ехії (0) 


Я an ASCII image is a list of character strings 


aimg = [] 


# generate the List of tile dimensions 
for j in range(rows): 
yl = int(j*h) 
y2 = int((j+1)*h) 
# correct the last tile 
if j == rows-1: 
y2 =H 
# append an empty string 
aimg.append("") 
for i in range(cols): 
# crop the image to fit the tile 
х1 = int(i*w) 
х2 = int((i-1)*w) 
# correct the last tile 
if i == со15-1: 
x2 = W 
# crop the image to extract the tile into another Image obj 
ect 
img = image.crop((xl, yl, x2, y2)) 
# get the average luminance 


avg = int(getAverageL (img) ) 


# look up the ASCII character for grayscale value (avg) 


if moreLevels: 

gsval = gscalel[int((avg*69)/255)] 
else: 

gsval = gscale2[int((avg*9)/255)] 
# append the ASCII character to the string 


aimg[j] += gsval 


# return text image 


return aimg 


# main() function 
def main(): 


# create parser 


descStr = "This program converts an image into ASCII art." 
parser = argparse.ArgumentParser(description=descStr) 
# add expected arguments 


parser.add argument('-- 
file', dest-'imgFile', required-True) 


parser.add argument('-- 
scale', dest-'scale', required-False) 


parser.add argument('-- 
out', dest-'outFile', required-False) 


parser.add argument('-- 
cols', dest='cols', required-False) 


parser.add argument('-- 
morelevels', dest='moreLevels', action='store true’) 


# parse arguments 


args = parser.parse args() 


imgFile = args.imgFile 
# set output file 
outFile = 'out.txt' 
if args.outFile: 
outFile = args.outFile 
# set scale default as 0.43, which suits a Courier font 
scale = 0.43 
if args.scale: 
scale = float(args.scale) 
# set cols 
cols = 80 
if args.cols: 
cols = int(args.cols) 
print('generating ASCII art...') 
# convert image to ASCII text 


aimg = covertImageToAscii(imgFile, cols, scale, args.moreLe 
vels) 


4 open a new text file 
f = open(outFile, 'w') 
# write each string in the list to the new file 
for row in aimg: 
f.write(row + '\n') 
# clean up 
f.close() 


print("ASCII art written to 555" % outFile) 


# call main 
if name == ' main ': 


main() 


6.5 [IASC 


ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΗΠΠΠΠΠΠΠΠΠαΦἴα/Γοῦοῖ.Ιρα[] 
ООО00000000000000 


$ python ascii.py --file data/robot.jpg --cols 100 


06-ЗПОДА5СЄЛООДОООДОДОгобогї.)Р9)000 





06-3 аѕсіі.руП0000 
ПООООО00000А5СИОО0000 


6.6 ПП 


ОО00000000000000000000000А5С0000000 
ООООО0000000000000000000000000000000 
ОДООДООД00000А5СПООДОД00000000 


6.7 ΠΠ 
ОО00000000000000А5СП0000 


1ΠΠΠΠΠΠΠ--5ς8!61.ΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠ5ς8!ἑΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 


ОООО000000000000000 


200000000000--іпмегДОДА5СПОДООООД00000 
ОД00000000000000000255П0000000000 


ЗОООДОД00000000000000000000000000000 
О00000000000000А5СИО000000000 


python3 аѕсіі.ру --map "@$%^`." 


О0000000000600000000000А5С0000 
Η΄ ©"ППО00000“. "0000002550 


070 00000 





О00000000000007-20000000000000000000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


ОО00000000000000000000000000000000*0 
О“ОО00000000000000000000000000000000 
ο. μμ. 


О000000000000000007-2000000000000000 
ООО0О0000000000000000000000000000000 
ООО000000000000000000000000000000000 
ОО0000" 00000000000 


О000000000000РуєһопрбО000000000000000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 000000 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕ ΕΕ 


ОООДО0000000000000000 


. ΠΡΥ΄ΠΟΠΠΠΠΠΡΙΓΠΠΠΠΠΠ 
- ПОДОООДАЄВОП 


• 00000 
• ПОДОДО000000000000000 
-ΠΗΠΗΠΠΠΠΗΠΗΒΘΒΗΠ 





07-1 ДОО0000 


7.1 OUUU 


ООО0О0000000000000000000000000000000 
ОДООО000000000000000000000000000М x М 
ОМОДООМОДООООбОДОДО000000000000000 


1П00000000000000000000000 
2 ПО000000000000М х МОДОДОД 


ЗПОООО0000000000000000000 
4000000000000М x МООООООДОДОДООООП 
7.1.1 000000 


О007-2000000000000000М x МО000 
ШЕГШЫ 


(^w, МУ) 


і сову пп у) 


4— (iw, (ih) 
W 





07-2 ПО0000 


07-200000000000000000000000х00000000 
YOUU 


О0000000000000000000000010)0000000000 
О ему, 1)ОО0О0000001+ у, ()+1)*һ)ППП\МУПП 
ОО0000000000РІН1О000000000000000 


7.1.2 ΠΠΠΠΠ 


О00000000000000000000000000000000008 
О000000000000800000010,255 1000000000 М 
DODODERGBDUDUL 


"+ та осо тм gib gab gu Di b DO bb 
r.g.b "Ἢ ЛНР . 


ΠΠΠΠΠΠΒΘΒΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠΠΠΗΠΗΠΠΠΠΠΠΠΠΗΠΒΟΒΠΠΠΠΠΗΠΗΠΠΠΠΗΠΗ 


7.1.3 (ПП! 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ. 
ОДОДОООД00000000000000000А08 00000000 
ОДДОДОДАСЄВОДОГО 


ОДОДОДОД00000000000000А6 ВОБООО000000 
ОООО000000000000000000000000000000 


012 == V ( Ру = )2 "E (41 — go) и [ by =- b» )2 


ΠΠΠΠΠΠίΓι,σι, 03) 0 672, 92; 0 2) 0000 
ОДОДОООД000000898ОДОДОО0000000АСВОП 
ОООДОО000000000000000000000019000 





7.2 [ПП 


ООДООД0РИТомУОООДОООДООО00000000000000 
пипруд0000000 


7.3 || 


ООООО000000000000000000000000000000 
ВС ВООДОДОД0О00000000000000000000000 
MI μμ. 


7.3.1 ΠΠΠΠΠΠ 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕ 


def getImages(imageDir) : 
given a directory of images, return a list of Images 


Ө files = os.listdir(imageDir) 
images = [] 


for file in files: 


Ө 
filePath = os.path.abspath(os.path.join(imageDir, file)) 


try: 


4 explicit load so we don't run into resource crunch 
© fp = open(filePath, "rb") 

im = Image.open(fp) 

images. аррепа (іт) 

# force loading the image data from file 
о im. Тоад () 

# close the file 
6 fp.close() 

except: 
# skip 
print("Invalid image: 55" % (filePath, ) ) 


return images 


[je [rrnos.tistdir0TlimageDir n np d 
ОО000000000000000000000РІЫ Ιπη80 ΕΠ 


[Je[jos.path.abspath()[]os.path.join (LLL 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗΡΥΙΠΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 

ΠΠΠΠΠΠΠΓΝΟΟΙΡΘΓΠΠΠΠΠΠΠΠΠΠΠΠΗ 

ОсЛРосХо ам 00000000000000000000000000 

[Windows V ILinux[/[][] 


ОДОООООРИОІ ппадерородорорр0000 
Ітаде.ореп()0000000000000000000000000 
О0000000000000000000Руһпоһпо00000000 


Πιπιαςς.ορεΠ()ΠΗΠΠΠΙΡΠΠΡΙΙΟΠΗΠΠΗΠΠΠΠΠΗ 
ΕΕ ΕΕ ΕΕ 


Πθ[]ΠΠορεπ/()ΠΠΗΠΠΠΠΗΠΠΠΠΗΠΠΠΗΗΠΠΗΗ 
іпаде.ореп ОДДОО0000іп0000000 


ОДорепо00000000000Ф9000іпаде.іоад O00 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΡΕ ΕΕ 


Le rp m pp ap mp III! 


7.3.2 000000000000 


ООО000000000000000000000000000000000 
П000009еѓАуегадевев()00000000 


def getAverageRGB(image): 


return the average color value as (г, 9, b) for each input 
image 


# get each tile image as a numpy array 
Ө im = np.array(image) 

# get the shape of each input image 
Ө w,h,d = im.shape 


# get the average RGB value 


Ө return tuple(np.average(im.reshape(w*h, d), axis=0) ) 


Ljerjinumpyriinimagen ип bd 
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def splitImage(image, size): 


given the image and dimensions (rows, cols), return an m*n 
115% of images 


Ө у, H = image.size[0], image.size[l] 
Ө т, п = size 
€ у, Б = int(W/n), int(H/m) 

# image list 

imgs = [] 

# generate a list of dimensions 


for j in range(m): 


for i in range(n): 
# append cropped image 
ο imgs.append(image.crop((i*w, j*h, (131) Ἐν, (j+1)*h))) 


return imgs 
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def getBestMatchIndex(input_avg, avgs): 


return index of the best image match based on average RGB v 
alue distance 


# input image average 


avg = input_avg 


# get the closest RGB value to input, based on RGB distance 


index = 0 
Ө min_index = 0 
Ө min dist = float("inf") 
Ө for val in avgs: 
о dist = ((val[0] - avg[0])*(val[0] - avg[0]) + 
(val[1] - avg[1])*(val[1] - avg[1]) + 


(val[2] - avg[2])*(val[2] - avg[2])) 
e if dist « min dist: 
min dist - dist 
min index = index 


index += 1 


return min index 
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def createImageGrid(images, dims): 


given a list of images and a grid size (m, n), create a gri 
d of images 


€ m, n - dims 


# sanity check 


assert m*n == len(images) 


# get the maximum height and width of the images 
# don't assume they're all equal 
Ө width = max([img.size[0] for img in images] ) 


height = max([img.size[1] for img in images] ) 


# create the target image 


Ө grid img = Image.new('RGB', (n*width, m*height) ) 


4 paste the tile images into the image grid 
for index in range(len(images) ): 
о гом = int(index/n) 


о col = index - n*row 


о 
grid img.paste(images[index], (col*width, row*height) ) 


return grid img 


бепрородррооордазвентоо0000 
сгеаїѓеітадебгіа()00000000000000аѕѕегі 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕΝ 
броброророровкеврородрдбо00000000000000 
ПООДОО0ОД000000000000000000000000000 
ПООДД0000000000000 


ОӨПОПО0О000000000000000000000000000 
ООО0О0000000000000000000000000000000 
ОООО000000000000000000000000 


Lern bimag eni nib m m p ab 
О0000000000000000өП0000000000000 
Ітаде.раѕёе()0000000000000000 
Ітаде.раѕїе()00000000000тааеПо000000 
ООО00000Д0000000000000000000000000000 
ОО000000000000000000000000000003%гоу 


+ сої) ОДОМОДОДОО00( ому, СО!)ПППППШП 
ΠΠΠΘΠΠΠΠΠΗΠΠΠΠΠΠΗΘΠΠΠΠΠΠ 


7.3.6 [| III 


τως 
[| 


def createPhotomosaic(target image, input images, grid size 
, reuse images-True): 


creates a photomosaic given target and input images 


print('splitting input image...') 
# split the target image into tiles 


Ө target images = splitImage(target image, grid size) 


print('finding image matches...') 

# for each tile, pick one matching input image 
output images - [] 

# for user feedback 

count = 0 


O batch size = int(len(target images)/10) 


# calculate the average of the input image 
avgs = [] 
for img in input images: 


© avgs.append(getAverageRGB(img) ) 


for img in target images: 
# compute the average RGB value of the image 
ο avg = getAverageRGB(img) 
# find the matching index of closest RGB value 
# from a list of average RGB values 
match index = getBestMatchIndex(avg, avgs) 
output images.append(input images [match index] ) 
# user feedback 
ο 
S count > 0 and batch size > 10 and count % batch size is 


print('processed %d of %d...' 56 
(count, len(target images) ) ) 


count += 1 
# remove the selected image from input if flag set 
о if not reuse images: 


input images. remove (match) 


print('creating mosaic...') 


4 create photomosaic image from tiles 


© mosaic image = createImageGrid(output images, grid size) 


# display the mosaic 


return mosaic image 
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# parse arguments 
parser = argparse.ArgumentParser(description='Creates a pho 
tomosaic from 
input images’ ) 
# add arguments 


parser.add argument('--target- 
image', dest='target image', required=True) 


parser.add argument('--input- 
folder', dest='input folder', required=True) 


parser.add argument('--grid- 
size', nargs-2, dest-'grid size', required-True) 


parser.add argument('--output 
file', dest='outfile', required-False) 
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print('resizing images...') 
# for given grid size, compute the maximum width and height 
of tiles 
e dims = (int(target image.size[0]/grid size[1]), 
int(target image.size[1]/grid 517е[0])) 
print("max tile dims: %5" % (dims,)) 
# resize 
for img in input images: 


e img.thumbnail(dims) 
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import sys, os, random, argparse 

from PIL import Image 

import imghdr 


import numpy as np 


def getAverageRGB(image): 


return the average color value as (г, g, b) for each input 
image 


# get each tile image as a numpy array 
im = np.array(image) 

# get the shape of each input image 
w,h,d = im.shape 

# get the average RGB value 


return tuple(np.average(im.reshape(w*h, d), axis=0) ) 


def splitImage(image, size): 


given the image and dimensions (rows, cols), returns an m*n 
list of images 


= 
= 
II 


image.size[0], image.size[1] 
m, n = size 

w, h = int(W/n), int(H/m) 

# image list 


1πος = [] 


# generate а list of dimensions 
for j in range(m): 
for i in range(n): 


# append cropped image 


imgs.append(image.crop((i*w, j*h, (i+1)*w, (j+1)*h))) 


return imgs 


def getImages(imageDir) : 


given a directory of images, return a list of Images 
files = os.listdir(imageDir) 
images = [] 


for file in files: 


filePath = os.path.abspath(os.path.join(imageDir, file) ) 


try: 


# explicit load so we don't run into a resource crunch 
fp = open(filePath, "rb") 
im = Image. ореп (Тр) 
images . аррепа (іт) 


# force loading image data from file 


im. load() 
# close the file 
fp.close() 
except: 
# skip 
print("Invalid image: ο" 5 (filePath, ) ) 


return images 


def getImageFilenames(imageDir) : 


given a directory of images, return a list of image filenam 
es 


files = os.listdir(imageDir) 
filenames = [] 


for file in files: 


filePath = os.path.abspath(os.path.join(imageDir, file) ) 
try: 
imgType = imghdr.what(filePath) 
if imgType: 
filenames.append(filePath) 


except: 


# skip 
print("Invalid image: 55" % (filePath, ) ) 


return filenames 


def getBestMatchIndex(input avg, avgs): 


return index of the best image match based on average RGB v 
alue distance 


# input image average 


avg = input avg 


# get the closest RGB value to input, based on RGB distance 
index = 0 
min index = 0 
min dist = float("inf") 
for val in avgs: 
dist = ((val[0] - avg[0])*(val[0] - avg[0]) + 
(val[1] - avg[1])*(val[1] - avg[1]) + 
(val[2] - avg[2])*(val[2] - avg[2])) 
if dist < min dist: 


min dist = dist 


min index = index 


index += 1 


return min_index 


def createImageGrid(images, dims): 


given a list of images and a grid size (m, n), create a gri 
d of images 


m, n = dims 


# sanity check 


assert m*n == len(images) 


# get the maximum height and width of the images 
# don't assume they're all equal 
width = max([img.size[0] for img in images] ) 


height = max([img.size[1] for img in images] ) 


# create the target image 


grid img = Image.new('RGB', (n*width, m*height) ) 


4 paste the tile images into the image grid 
for index in range(len(images) ): 
row = int(index/n) 


col = index - n*row 


grid img.paste(images[index], (col*width, row*height) ) 


return grid img 


def createPhotomosaic(target image, input images, grid size 
, reuse images-True): 


creates photomosaic given target and input images 


print('splitting input image...') 
# split the target image into tiles 


target images - splitlImage(target image, grid size) 


print('finding image matches...') 


# for each tile, pick one matching input image 


output images = [] 
# for user feedback 
count = 0 


batch size = int(len(target images) /10) 


# calculate the average of the input image 
avgs = [] 
for img in input images: 


avgs.append (getAverageRGB (img) ) 


for img in target images: 
# compute the average RGB value of the image 
avg = getAverageRGB(img) 
# find the matching index of closest RGB value 
# from a list of average RGB values 
match index = getBestMatchIndex(avg, avgs) 
output images.append(input images[match index] ) 
# user feedback 
ч count > 0 and batch size > 10 and count % batch size is 


print('processed %d of %d...' % 
(count, len(target images))) 


count += 1 


# remove the selected image from input if flag set 
if not reuse images: 


input images. remove (match ) 


print('creating mosaic...') 
# create photomosaic image from tiles 


mosaic image = createImageGrid(output images, grid size) 


# display the mosaic 


return mosaic image 


# gather our code in a main() function 


def main(): 


# command line arguments are in sys.argv[1], sys.argv[2], 


# sys.argv[0] is the script name itself and can be ignored 
# parse arguments 
parser = argparse.ArgumentParser(description='Creates a pho 
tomosaic from 
input images' ) 
# add arguments 


parser.add argument('--target- 
image', dest-'target image', required=True) 


parser.add argument('--input- 
folder', dest='input folder', required=True) 


parser.add argument('--grid- 
size', nargs-2, dest-'grid size', required-True) 


parser.add argument('--output- 
file', dest='outfile', required-False) 


args - parser.parse args() 


###### INPUTS ###### 


# target image 


target image = Image.open(args.target image) 


# input images 
print('reading input folder...') 


input images = getImages(args.input folder) 


# check if any valid input images found 

if input images == []: 
print('No input images found in %s. Exiting.' % (args.input 
_folder, )) 


exit () 


# shuffle list to get a more varied output? 


random.shuffle(input images) 


# size of the grid 


grid size = (int(args.grid size[0]), int(args.grid size[1]) 
) 


# output 
output filename = 'mosaic.png' 
if args.outfile: 


output filename - args.outfile 


# reuse any image in input 


reuse images - True 


# resize the input to fit the original image size? 


resize input - True 


##### END INPUTS ##### 


print('starting photomosaic creation...') 


# if images can't be reused, ensure m*n <= num of images 
if not reuse images: 
if grid size[0]*grid size[1] > len(input images): 
print('grid size less than number of images') 
exit() 
# resizing input 
if resize input: 
print('resizing images...') 
# for given grid size, compute the maximum width and height 
of tiles 
dims = (int(target image.size[0]/grid size[1]), 
int(target image.size[1]/grid size[0])) 
print("max tile dims: %s" % (dims,)) 
# resize 
for img in input images: 


img.thumbnail(dims) 


4 create photomosaic 
mosaic image - createPhotomosaic(target image, input images 
, grid size, 


reuse images) 


4 write out mosaic 


mosaic image.save(output filename, 'PNG') 


print("saved output to %s" % (output filename,)) 


print('done.') 


4 standard boilerplate to call the main() function 
# to begin the program 
if name == ' main ': 


main() 


7.5 ПОО00000000 
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$ python photomosaic.py --target-image test- 
data/cherai.jpg --input-folder 


test-data/set6/ --grid-size 128 128 
reading input folder... 

starting photomosaic creation... 
resizing images... 

max tile dims: (23, 15) 

splitting input image... 

finding image matches... 

processed 1638 of 16384 ... 


processed 3276 of 16384 ... 


processed 4914 of 16384 ... 
Creating mosaic... 
saved output to mosaic.png 


done. 
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4 tile a graphics file to create an intermediate image of а 
set size 


def createTiledImage(tile, dims): 
# create the new image 

Ө img = Image.new('RGB', dims) 
W, Н = dims 
w, h = tile.size 


# calculate the number of tiles needed 


© 


cols = int(W/w) + 1 


© 


rows = int(H/h) + 1 
# paste the tiles into the image 
for i in range(rows): 


for j in range(cols): 


о img.paste(tile, (j*w, 1*һ)) 
# output the image 


return img 
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я create an image tile filled with random circles 
def createRandomTile(dims) : 
# create image 
Ө img = Image.new('RGB', dims) 
O draw = ImageDraw.Draw(img) 
# set the radius of a random circle to 1% of 
# width or height, whichever is smaller 


© r = int(min(*dims)/100) 


4 number of circles 
Ө n= 1000 
# draw random circles 


for i in range(n): 


# - 
r makes sure that the circles stay inside and aren't cut of 
f 


# at the edges of the image so that they'll look better whe 
n tiled 


Ө X, у = random.randint(0, dims[0]- 

г), random.randint(0, dims[1]-r) 

о 

fill = (random.randint(0, 255), random.randint(0, 255), 
random. гапаіпї (0, 255)) 

ο draw.ellipse((x-r, y-r, х+г, ytr), fill) 


return img 
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>>> import autos 
>>> img = autos.createRandomTile( (256, 256) ) 
>>> img.save('out.png') 


>>> exit() 
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8.3.3 ΠΠΠΠΠΠΠ 
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ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕ 


def createAutostereogram(dmap, tile): 
# convert the depth map to a single channel if needed 
Ө if dmap.mode is not 'L': 
dmap = dmap.convert('L') 
# if no image is specified for a tile, create a random circ 
les tile 
Ө if not tile: 
tile = createRandomTile((100, 100) ) 
# create an image by tiling 
© img = createTiledImage(tile, dmap.size) 
# create a shifted image using depth map values 
Ө sImg = img.copy() 
# get access to image pixels by loading the Image object fi 
rs 
Ө pixD = dmap.load() 
pixS = sImg.load() 


# shift pixels horizontally based on depth map 


Ө cols, rows = sImg.size 
for j in range(rows): 
for i in range(cols): 
xshift = pixD[i, 11/10 
xpos = i - tile.size[0] + xshift 


if xpos » 0 and xpos « cols: 


e o о o 


pixS[i, j] = pixS[xpos, 1] 
# display the shifted image 


return sImg 
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# create a parser 
parser = argparse.ArgumentParser(description="Autosterogram 
Ses) 


# add expected arguments 


Ө parser.add argument(' -- 
depth', dest-'dmFile', required-True) 


parser.add argument('-- 
tile', dest-'tileFile', required-False) 


parser.add argument('-- 
out', dest-'outFile', required-False) 


# parse args 
args = parser.parse args() 
# set the output file 
outFile = 'as.png' 
if args.outFile: 

outFile = args.outFile 
# set tile 
tileFile = False 
if args.tileFile: 


tileFile = Image.open(args.tileFile) 
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8.4 ПП 


ΕΕ ΕΕ ΕΕ ΕΕΒΕ ΕΒΕ 
https://github.com/electronut/pp/blob/ 
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import sys, random, argparse 


from PIL import Image, ImageDraw 


# create spacing/depth example 


def createSpacingDepthExample(): 


tiles = [Image.open('test/a.png'), Image.open('test/b.png') 


Image.open('test/c.png')] 
img - Image.new('RGB', (600, 400), (0, 0, 0)) 
spacing = [10, 20, 40] 
for j, tile in enumerate(tiles): 
for i in range(8): 
img.paste(tile, (10 + i*(100 + j*10), 10 + j*100)) 


img.save('sdepth.png') 


# create an image filled with random circles 
def createRandomTile(dims): 

# create image 

img - Image.new('RGB', dims) 

draw - ImageDraw.Draw(img) 

# set the radius of a random circle to 1% of 

# width or height, whichever is smaller 

г = int(min(*dims)/100) 

# number of circles 

n = 1000 

# draw random circles 


for i in range(n): 


# - 
r makes sure that the circles stay inside and aren't cut of 
f 


# at the edges of the image so that they'll look better whe 
n tiled 
X, у = random.randint(0, dims[0]- 
r), random.randint(0, dims[1]-r) 
fill = (random.randint(0, 255), random.randint(0, 255), 
random. гапалп (0, 255)) 
draw.ellipse((x-r, y-r, x+r, y+r), fill) 
# return image 


return img 


# tile a graphics file to create an intermediate image of a 
set size 


def createTiledImage(tile, dims): 
# create the new image 
img - Image.new('RGB', dims) 
W, Н = dims 
w, h = tile.size 


# calculate the number of tiles needed 


cols int(W/w) + 1 
rows = int(H/h) + 1 


# paste the tiles into the image 


for i in range(rows): 
for j in range(cols): 
img.paste(tile, (j*w, i*h)) 
# output the image 


return img 


# create a depth map for testing 

def createDepthMap (dims): 
dmap = Image.new('L', dims) 
dmap.paste(10, (200, 25, 300, 125)) 
dmap.paste(30, (200, 150, 300, 250) ) 
dmap.paste(20, (200, 275, 300, 375)) 


return dmap 


# given a depth map image and an input image, 
# create a new image with pixels shifted according to depth 
def createDepthShiftedImage(dmap, img): 

# size check 


assert dmap.size == img.size 


# create shifted image 


sImg = img.copy() 


# get pixel access 
pixD = παρ. load() 
pixS = sImg.load() 
# shift pixels output based on depth map 
cols, rows = sImg.size 
for j in range(rows) : 
for i in range(cols): 
xshift = pixD[i, j]/10 
xpos = i - 140 + xshift 
if xpos > 0 and xpos < cols: 
pixS[i, j] = pixS[xpos, j] 
# return shifted image 


return sImg 


# given a depth map (image) and an input image, 
# create a new image with pixels shifted according to depth 
def createAutostereogram(dmap, tile): 

# convert the depth map to a single channel if needed 

if dmap.mode is not 'L': 

dmap = dmap.convert('L') 

# if no image is specified for a tile, create a random circ 
les tile 


if not tile: 


tile = createRandomTile((100, 100)) 
# create an image by tiling 
img = createTiledImage(tile, dmap.size) 
# create a shifted image using depth map values 
sImg = img.copy() 


# get access to image pixels by loading the Image object fi 
rst 


pixD = dmap.load() 
pixS = sImg.load() 
# shift pixels horizontally based on depth map 
cols, rows = sImg.size 
for j in range(rows): 
for i in range(cols): 
xshift = pixD[i, j]/10 
xpos = i - tile.size[0] + xshift 
if xpos > 0 and xpos « cols: 
pixS[i, j] = pixS[xpos, j] 
# return shifted image 


return sImg 


# main() function 


def main(): 


# use sys.argv if needed 

print('creating autostereogram...') 

# create parser 
parser = argparse.ArgumentParser(description="Autosterogram 
Suc) 

# add expected arguments 


parser.add argument('-- 
depth', dest-'dmFile', required-True) 


parser.add argument('-- 
tile', dest='tileFile', required-False) 


parser.add argument('-- 
out', dest='outFile', required-False) 


# parse args 
args = parser.parse args() 
# set the output file 
outFile = 'as.png' 
if args.outFile: 

outFile = args.outFile 
# set tile 
tileFile = False 
if args.tileFile: 

tileFile = Image.open(args.tileFile) 
# open depth map 


dmImg = Image. open(args.dmFile) 


# create stereogram 
asImg = createAutostereogram(dmImg, tileFile) 
# write output 


asImg.save(outFile) 


# call main 
if name == ' main ': 


main() 


8.5 ΠΠΠΠΠΗΠΗΠΠΗ 
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$ python3 autos.py --depth data/stool-depth.png 
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$ python3 autos.py --depth data/stool-depth.png – 
tile data/escher-tile. jpg 
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geometry.hubpages.com/hub/Free-M-C- 
Escher-Tessellation-Background-Patterns- 
Tiling-Lizard-Background/ 
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import sys 
from OpenGL.GLUT import * 


from OpenGL.GL import * 


def display(): 
glClear (GL COLOR BUFFER BIT|GL DEPTH BUFFER BIT) 
glColor3f (1.0, 1.0, 0.0) 


glBegin(GL_ QUADS) 


glVertex3f (-0.5, -0.5, 0.0) 
glVertex3f (0.5, -0.5, 0.0) 
glVertex3f (0.5, 0.5, 0.0) 
glVertex3f (-0.5, 0.5, 0.0) 
glEnd() 


glFlush(); 


glutInit(sys.argv) 
glutInitDisplayMode(GLUT SINGLE|GLUT RGB) 
glutInitWindowSize(400, 400) 
glutCreateWindow("oldgl") 
glutDisplayFunc(display) 


glutMainLoop() 
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Ө #version 330 соге 
Ө in vec3 aVert; 


© uniform mat4 uMVMatrix; 


Ө uniform mat4 uPMatrix; 
Ө out vec4 vCol; 


void main() { 
// apply transformations 
@ gl Position = uPMatrix * uMVMatrix * vec4(aVert, 1.0); 
// set color 
@ vCol = vec4(1.0, 0.0, 0.0, 1.0); 


} 
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Ө #version 330 core 
Ө in vec4 vCol; 
© out vec4 fragColor; 
void main() { 
// use vertex color 


Ө fragColor = vCol; 
} 
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class RenderWindow: 
"""GLFW Rendering window class""" 


def init (self): 


# save current working directory 


сма = os.getcwd() 


# initialize glfw 


e glfw.glfwlnit() 


# restore cwd 


os.chdir(cwd) 


# version hints 
e 
glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MAJOR, 3) 


glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MINOR, 3) 


glfw.glfwwindowHint(glfw.GLFW OPENGL FORWARD COMPAT, GL TRU 
E) 


gl fw.gl fwWindowHint (glfw.GLFW OPENGL PROFILE, 


gl fw.GLFW OPENGL CORE PROFILE) 


# make a window 
self.width, self.height = 640, 480 
self.aspect = self.width/float(self.height) 


e 
self.win = glfw.glfwCreateWindow(self.width, self.height, 


b'simpleglfw') 


4 make the context current 


о gl fw.glfwMakeContextCurrent(self.win) 


HE QUUUGLFW I I e@ 00000 pen GET TETETETL] 
OpenGL 3.30000000Ф0000000640 x 480000 
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OpenGL] 


ООО00000000000 
# initialize GL 
ө glViewport(0, 0, self.width, self.height) 
e glEnable(GL DEPTH TEST) 


e glClearColor(0.5, 0.5, 0.5, 1.0) 
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О00006РМПОООООООООООВООООООООО0000 
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4 set window callbacks 


glfw.glfwSetMouseButtonCallback(self.win, self.onMouseButto 


glfw.glfwSetKeyCallback(self.win, self.onKeyboard) 


glfw.glfwSetWindowSizeCallback(self.win, self.onSize) 


ООО000000000000000000000000000000000 
ООО000000000000 


0000 
0000000 


def onKeyboard(self, win, key, scancode, action, mods): 
#print 'keyboard: ', win, key, scancode, action, mods 
e if action -- glfw.GLFW PRESS: 
4 ESC to quit 
if key == glfw.GLFW KEY ESCAPE: 
e self.exitNow = True 
else: 


# toggle cut 


self.scene.showCircle = not self.scene.showCircle 


О000000000опКеуроага()000000000000000 
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def onSize(self, win, width, height): 
#print 'onsize: ', win, width, height 
self.width = width 
self.height - height 
self.aspect - width/float(height) 


e glViewport(0, 0, self.width, self.height) 


LDOBHOBUOBODUS I ViewportO ООДОО0000000000 
ороророФр0р000000Оміаєо height 000000 
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def run(self): 
# initializer timer 
e glfw.glfwSetTime(0) 
t = 0.0 
e 
while not glfw.glfwWindowShouldClose(self.win) and not self 
.exitNow: 


4 update every x seconds 


e currT = glfw.glfwGetTime() 


if currT - t > 0.1: 
4 update time 
t = currT 
# clear 


ο 
glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


# build projection matrix 


e 
pMatrix = glutils.perspective(45.0, self.aspect, 0.1, 100.0 
) 


о 
mvMatrix = glutils.lookAt([0.0, 0.0, -2.01, [0.0, 0.0, 0.0) 


, 


[0.0, 1.0, 0.0]) 


4 render 

ο self.scene.render(pMatrix, mvMatrix) 
# step 

Ө self.scene.step() 

о glfw.glfwSwapBuffers(self.win) 


# poll for and process events 
Φ glfw.glfwPollEvents() 


# end 


gl fw.glfwlerminate() 
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ПОСОООООООООООООООООоОООоОоооооооооооооо0000000000 
ОДООД0000000000000000000000000000000000000000000 


О000000000000000000005П0000000000000000000000 
σι ΓΝΝΗΠΠΠΗ 


9.4.3 Scene|] 
П00005сепероо000000000000000 


class Scene: 
"" OpenGL 3D scene class""" 
# initialization 
def init (self): 
# create shader 


e self.program - glutils.loadShaders(strVS, strFS) 


6 glUseProgram(self.program) 


О5сеперроодроОоПО000000000000000000000 
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self.pMatrixUniform = glGetUniformLocation(self.program, b' 
uPMatrix' ) 


self .mvMatrixUniform = glGetUniformLocation(self.program, b 


'uMVMatrix') 


4 texture 


self.tex2D = glGetUniformLocation(self.program, b'tex2D') 


OUUUUUgIGEetUniformLocation )ПОООООО 
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я define triangle strip vertices 
e vertexData = numpy.array ( 
[-0.5, -0.5, 0.0, 
0.5, -0.5, 0.0, 
-0.5, 0.5, 0.0, 


0.5, 0.5, 0.0], numpy.float32) 


# set up vertex array object (\/АО) 

e self.vao = glGenVertexArrays(1) 
glBindVertexArray(self.vao) 
# vertices 

e self.vertexBuffer = glGenBuffers(1) 
glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


# set buffer data 


ο 
glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexData 


GL_STATIC DRAW) 
# enable vertex array 
e glEnableVertexAttribArray(0) 


4 set buffer data pointer 


о 
glVertexAttribPointer(0, 3, GL FLOAT, GL FALSE, 0, None) 
à unbind МАО 


ο glBindVertexArray(0) 
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# texture 


self.texId = glutils.loadTexture('star.png' ) 
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# step 
def step(self): 
# increment angle 
о self.t = (self.t + 1) % 360 


# set shader angle in radians 


e 
glUniformlf(glGetUniformLocation(self.program, 'uTheta'), 


math.radians(self.t)) 
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# render 
def render(self, pMatrix, mvMatrix): 


# use shader 


e glUseProgram(self.program) 


# set projection matrix 
e 


glUniformMatrix4fv(self.pMatrixUniform, 1, GL FALSE, pMatri 
x) 


# set modelview matrix 


glUniformMatrix4fv(self.mvMatrixUniform, 1, GL FALSE, mvM 
atrix) 


# show circle? 
e 
glUniformli(glGetUniformLocation(self.program, b'showCircle 
'), 


self.showCircle) 


# enable texture 
glActiveTexture(GL TEXTUREO) 
glBindTexture(GL TEXTURE 2D, self.texId) 


glUniformli(self.tex2D, 0) 


# bind VAO 


ο glBindVertexArray(self.vao) 


# draw 
Ө glDrawArrays(GL TRIANGLE STRIP, 0, 4) 
à unbind МАО 


© glBindVertexArray (0) 
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#version 330 core 
€ layout(location = 0) in vec3 aVert; 


Ө uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 


uniform float uTheta; 
© out vec2 vTexCoord; 


void main() { 
// rotational transform 

Ө mat4 rot = mat4( 
vec4(cos(uTheta), sin(uTheta), 0.0, 0.0), 
vec4(-sin(uTheta), cos(uTheta), 0.0, 0.0), 
vec4(0.0, 0.0, 1.0, 0.0), 
vec4(0.0, 0.0, 0.0, 1.0) 
); 


// transform vertex 
gl Position = uPMatrix * uMVMatrix * rot * vec4(aVert, 1.0) 


// set texture coordinate 
Ө vTexCoord = aVert.xy + vec2(0.5, 0.5); 


} 
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#version 330 соге 


Ө in vec4 vCol; 


in vec2 vTexCoord; 


Ө uniform sampler2D tex2D; 


© uniform bool showCircle; 
Ө out vec4 fragColor; 


void main() { 
if (showCircle) { 
// discard fragment outside circle 
e if (distance(vTexCoord, vec2(0.5, 0.5)) > 0.5) { 


discard; 


е1<е { 
о fragColor = texture(tex2D, vTexCoord); 


} 


else { 

© fragColor = texture(tex2D, vTexCoord) ; 
} 

} 
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import OpenGL 


from OpenGL.GL import * 


import numpy, math, sys, os 


import glutils 
import glfw 


st rVS = nn 


#version 330 core 
layout(location = 0) in vec3 aVert; 


uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 


uniform float uTheta; 
out vec2 vTexCoord; 


void main() { 

// rotational transform 

mat4 rot = mat4( 
vec4(cos(uTheta), sin(uTheta), 0.0, 0.0), 
vec4(-sin(uTheta), cos(uTheta), 0.0, 0.0), 
vec4(0.0, 0.0, 1.0, 0.0), 
vec4(0.0, 0.0, 0.0, 1.0) 
); 


// transform vertex 


gl Position = uPMatrix * uMVMatrix * rot * vec4(aVert, 1.0) 


, 


// set texture coordinate 


vTexCoord = aVert.xy + vec2(0.5, 0.5); 


} 


strFS = nuu 


Zversion 330 core 


in vec2 vTexCoord; 


uniform sampler2D tex2D; 


uniform bool showCircle; 


out vec4 fragColor; 


void main() { 
if (showCircle) { 
// discard fragment outside circle 


if (distance(vTexCoord, vec2(0.5, 0.5)) > 0.5) { 


discard; 
} 
else { 
fragColor = texture(tex2D, vTexCoord) ; 
} 
} 
else { 
fragColor = texture(tex2D, vTexCoord) ; 
} 


} 


class Scene: 
""" OpenGL 3D scene class""" 
Я initialization 
def init (self): 
# create shader 


self.program = glutils.loadShaders(strVS, strFS) 


glUseProgram(self.program) 


self.pMatrixUniform = glGetUniformLocation(self.program, b' 
uPMatrix' ) 


self .mvMatrixUniform = glGetUniformLocation(self.program, b 
'uMVMatrix') 


4 texture 


self.tex2D = glGetUniformLocation(self.program, b'tex2D') 


# define triange strip vertices 
vertexData = numpy.array ( 
[-0.5, -0.5, 0.0, 
0.5, -θ.5, 0.0, 
-0.5, 0.5, 0.0, 


0.5, 0.5, 0.0], numpy.float32) 


4 set up vertex array object (\/АО) 

self.vao = glGenVertexArrays(1) 
glBindVertexArray(self.vao) 

# vertices 

self.vertexBuffer = glGenBuffers(1) 
giBindBuffer(GL ARRAY BUFFER, self.vertexBuffer) 


4 set buffer data 


glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexData 


, 


GL STATIC DRAW) 
4 enable vertex array 
glEnableVertexAttribArray(0) 


# set buffer data pointer 


glVertexAttribPointer(0, 3, GL FLOAT, GL FALSE, 0, None) 
à unbind МАО 


glBindVertexArray(0) 


# time 
self.t = 0 
# texture 


self.texId = glutils.loadTexture('star.png' ) 


# show circle? 


self.showCircle = False 


# step 
def step(self): 
4 increment angle 
self.t = (self.t + 1) % 360 


4 set shader angle in radians 
glUniformlf(glGetUniformLocation(self.program, 'uTheta'), 


math.radians(self.t)) 


# render 
def render(self, pMatrix, mvMatrix): 
# use shader 


glUseProgram(self.program) 


# set projection matrix 


glUniformMatrix4fv(self.pMatrixUniform, 1, GL FALSE, pMatri 


x) 


# set modelview matrix 


glUniformMatrix4fv(self.mvMatrixUniform, 1, GL FALSE, mvMat 
rix) 


# show circle? 
glUniformli(glGetUniformLocation(self.program, b'showCircle 
"), 


self.showCircle) 


# enable texture 
glActiveTexture(GL TEXTUREO) 
glBindTexture(GL TEXTURE 2D, self.texId) 


glUniformli(self.tex2D, 0) 


3 bind VAO 
glBindVertexArray(self.vao) 

# draw 

glDrawArrays(GL TRIANGLE STRIP, 0, 4) 
à unbind МАО 


glBindVertexArray (0) 


class RenderWindow: 
"""GLFW Rendering window class""" 


def init (self): 


# save current working directory 


cwd = os.getcwd() 


# initialize glfw - this changes cwd 


glfw.glfwInit() 


4 restore cwd 


os.chdir(cwd) 


# version hints 
glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MAJOR, 3) 
glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MINOR, 3) 


glfw.glfwwindowHint(glfw.GLFW OPENGL FORWARD COMPAT, GL TRU 
E) 
gl fw.gl fwWindowHint (glfw.GLFW OPENGL PROFILE, 
glfw.GLFW OPENGL CORE PROFILE) 
# make a window 


self.width, self.height = 640, 480 


self.aspect = self.width/float(self.height) 


self.win = glfw.glfwCreateWindow(self.width, self.height, 


b'simplegl fw' ) 
# make context current 


gl fw.glfwMakeContextCurrent(self.win) 


# initialize GL 
glViewport(0, 0, self.width, self.height) 
glEnable(GL DEPTH TEST) 


glClearColor(0.5, 0.5, 0.5, 1.0) 


4 set window callbacks 


glfw.glfwSetMouseButtonCallback(self.win, self.onMouseButto 
n) 


glfw.glfwSetKeyCallback(self.win, self.onKeyboard) 


glfw.glfwSetWindowSizeCallback(self.win, self.onSize) 


# create 3D 


self.scene = Scene() 


# exit flag 


self.exitNow = False 


def onMouseButton(self, win, button, action, mods): 
#print ‘mouse button: ', win, button, action, mods 


pass 


def onKeyboard(self, win, key, scancode, action, mods): 
#print 'keyboard: ', win, key, scancode, action, mods 
if action -- glfw.GLFW PRESS: 
# ESC to quit 
if key == glfw.GLFW KEY ESCAPE: 
self.exitNow - True 
else: 


# toggle cut 


self.scene.showCircle = not self.scene.showCircle 


def onSize(self, win, width, height): 
#print 'onsize: ', win, width, height 
self.width = width 


self.height height 


self.aspect - width/float(height) 


glViewport(0, 0, self.width, self.height) 


def run(self): 


# initializer timer 
gl fw.gl fwSetTime(0) 
t = 0.0 
while not glfw.glfwWindowShouldClose(self.win) and not self 
.exitNow: 
# update every x seconds 
currT = glfw.glfwGetTime() 
if currT - t 20.1: 
# update time 
t = currT 


4 clear 


glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


4 build projection matrix 


pMatrix = glutils.perspective(45.0, self.aspect, 0.1, 100.0 
) 


mvMatrix = glutils.lookAt([0.0, 0.0, -2.0], [0.0, 0.0, 0.0] 


, 


[0.0, 1.0, 0.0]) 


# render 


self.scene.render(pMatrix, mvMatrix) 


# step 


self.scene.step() 


о Ћи. а fwSwapBuf fers (self .win) 
# poll for and process events 


glfw.glfwPollEvents() 


# end 


glfw.glfwTerminate() 


def step(self): 
4 clear 


glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


4 build projection matrix 


pMatrix = glutils.perspective(45.0, self.aspect, 0.1, 100.0 
) 


mvMatrix = glutils.lookAt([0.0, 0.0, -2.0], [0.0, 0.0, 0.0] 


, 


[0.0, 1.0, 0.0]) 
# render 


self.scene.render(pMatrix, mvMatrix) 


# step 


self.scene.step() 


glfw.SwapBuffers(self.win) 
# poll for and process events 


glfw.PollEvents() 


# main() function 
def main(): 
print("Starting simpleglfw. " 
"Press any key to toggle cut. Press ESC to quit.") 
rw = RenderWindow( ) 
rw.run() 
# call main 
if name == ' main ': 


main() 


9.6 [OpenGL 
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$ python simpleglfw.py 
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def loadTexture(filename): 

"""load OpenGL 2D texture from given image file""" 
img = Image.open(filename) 
imgData - numpy.array(list(img.getdata()), np.int8) 
texture = glGenTextures(1) 
glBindTexture(GL TEXTURE 2D, texture) 


glPixelStorei(GL UNPACK ALIGNMENT, 1) 


о о о о © © 


glTexParameterf (GL TEXTURE 2D, GL TEXTURE WRAP S, GL CLAMP 
Το EDGE) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE WRAP T, GL CLAMP 
Το EDGE) 
© 


glTexParameterf (GL TEXTURE 2D, GL TEXTURE MAG FILTER, GL L 
INEAR) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE MIN FILTER, GL_L 
INEAR) 

ө 

glTexlImage2D(GL TEXTURE 2D, ©, GL АСВА, img.size[0], ітд. 5 
ize[1], 


0, GL RGBA, GL UNSIGNED BYTE, imgData) 


return texture 


ЦӨППЦ!оааТехїиге()ЦЦЦПРу  оПЦЦЦПЦПР!ЦПЦ 
Ітадео00000000000өӨП000тадеђо000000 
ӘППпитруППӨППППППОрепое pp ET 
Орепс00000000000000ө000000000000 
ПОсехёигеПоО0000000000000000000000ө 
ПОО000000000010000000000000000100008 
ΠΠΠΠΠΠΠΦ/ΠΠΠΠΠΟΡρεπαΙ ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
ООДОО0000000000000000000000000005010 
ООбО000х0Уу000ФОДОДОО0О0О0О0О0000000000000 
0000000000000000000“ 0000“ оде оооооооо 
ПОДОО000000000000000000000000 


9.7 [I 


О0000000РуёћопОрепсі 000000000000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


9.8 [IL 
ООО000000000000 


100000000000000200000010000000000000 
О10100000000000000000000000000000000 
О000РуєһопО00000000000000чпіғогтр000 
ΕΕ ΕΕ ΕΕ ΕΕ 


2 П0000000000000000000000000000000000 
ООООО0000000000000000000000000000000 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕΝ 
ООООО000000000000004 х40000000000000 
1L.000000000g!TexParameterf()00 

GL TEXTURE WRAP ӘЛППІШІСІ. ΒΕΡΕΑΤΠΠ 


ЗООДОДО000000000000009-/0000000061511 
біл 0000 





09-7  0000000000000 


П101 [DUE 





ООО0О0000000000000000000000000000000 
ο. ο μ.μ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΣ 
ООО0О0000000000000000000000000000000 
ООО000000000000000000 


ООО00000Д0000000000000000000000000000 
ООО000000000000000000000000000000000 
ООО000000000000000000000000000000000 
О000000000000000000000001роагаітоП00 
О0000000000000000010-2000 


ООООО0000000000000000000000000000000 
По РОДООДОООДООООДОДО0О0О000000000000 
ООО000000000000000000000000000000000 
ОДОДОО0О0000000000000О0репс ППТ 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 





[10-1 ΠΠΠΠΠΠΠΠΠΠ 


ОООО0О0000000000000000000000000000000 
ОДООО00000000000000000000000000000000 
ОООО0О0000000000000000000000000000000 
ΕΕ ΕΕ ΕΕ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


• ПО0Д000000000 
- ПО РООООООПО 
- ПОДОДО0Д00000000000 


«ΏΠΟρεησ! ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΑΡΠΘΠΠΠΠΠΠ 
ООООО 
- ДОДОД000000000 


10.1 [I 


ООООО0000000000000000000000000000000 
О00000000000000000010-20000000000000 


010-2 50000000000000000 
ООО000000000000 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 

• П0000000000900020000000000000 

• ПООООО0000000000000000000000000 
• ПООДО000000000000000 

ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


10.1.1 [| III 
ΠΠΠΠΠΠΠΠΝΠΠΠΠΠ {ΠΠΠΠΠΠΠΠΠΠΠ 


I а 
Pi = Po + Wt аё 


000 P OOOO 20000 Р о ОДОДО0000 V о ОООООО 
Г ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


О000000000000000000000000000000000 
тшш z ППОООООООООООООООООНО 


10.1.2 ΠΠΠΠΠΠ 


ПИШШИШИШШИШИШИШИШШИШИП z О0000000000000 
ООО000000000000000000000000000000000 
ο. тле ОПП 





010-3 ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 


0010-300000 o 000000 х 00000000 6 00000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΡΕ 


ОДОО000000000000000000000000000000000 
ООО00000Д0000000000000000000000000000 
ОООО00000000000000000000000 


000 а О0900000000000000000200000000000 


ОД0000002.00000000000000 У ООО0000000 
О00000 


У = (cos(8) sin( Ф). sin e(8) sin( Ф). cos( Ф)) 


ΠΙΟΠΖΟΙΠΠΠΗΠΠΠΠΠΠΠΠΠΠΠΙΟΠ36ΟΙΠΠΠΠΠΠΠΗ 
ΕΠΗ 


9 = random! (0, 20] ).Ф = random! 0, 360] ) 


О00000000020-ЗОО0000000000000000000 
ОООО000000000000000 


ООО0000000000000000000010.1000000000 
ο ο. шшш ΓΠΗΠΗ 


Ба = 0.051 


О000000000000000000000000000.0500000 
О00200000000000000000000000000000000 


они 


10.1.3 (000 


О0000000000000000000000репец 
GL ΡΟΙΝΤΘΗΠΗΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΉΠΠΠ 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 0000000 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕ 
0000" 000" 000000000000000000000000000 
ООО0О0000000000000000000000000000000 
ООО000000000000000 


10.1.4 ПОрепсі 00000000 


П000000000000000репсО6000епаіпо00 
000000000000000000000000000000000000 
ПАІрһаППП 


ОООДО0000000000000000000000000000 
атрпаДОбОДОД000000000000000Оаїрвадоо( 
ООО000000000000000000000000000000000 
ΠΠΠΠΠΠΠΠΠΠΠ8!ΡΠ8ΠΠΠΠΠΠΠΠΠΠΠΒΗΘΒΑΠΠΠ 
32ΠΠΒΟΒΑΠΠΠΠΠΘΙΡΠΘ8ΠΠΠΠΠΙΘ,255ΙΠΟΠΠΠΠΗ 
ΠΠ255ΠΠΠΠΠΠΠ8!ΡΠ8ΠΠΠΠΠΠΠΠΠΠΠΠΠΠ8!ΡΠ8 
ОДДОООД00АЄвОрО0О000000000 


OpenGL p b Op E Od OU C CCIHIDUL 
10-400000000000000000000000000000 


ППОрепс:1000000аірпароб0000000000000 
О0000000000866Вво000000000000000арһа 
О00000000000000000000000000000000000 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ8!ΡΠΘΠΠΠΠΠΠΠ 
0000010. 3.50000000 


~ RGB = (0, 0, 0) 
At RGB = (255, 255, 255) 


Π10-4 η 





ΠΠΠΠΠΠΠΠΠ4!ρΠ8ΠΠΠΠΠΠΠΠΠΠΠΠΠ4!ΡΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
Θ!ΡΠ8ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ8!ΡΠ8ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘ!Ρ Π8ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


10.1.5 [JII 


ООО000000000000000000000000000000000 
ОООО00000000000000000000000000000000 


ООО000000000000000000000000000000000 
ОООООД0000000000000000000000000000000 
О0РОДО"0О000000 


О00000000000000000000000000000010-5 
ОООДО0000000000000000 


и u 


GENE LEE) 





(ZARR Міла 75 PIER) 


[10-5 DON 


ООО000000000000000000000000000000000 
ООООО0000000000000000000000000000000 
U п ОООООООД00000000000000000000000000 
ООООО0000000000000000000000000000000 
v 000 п П0О000000000000000000000000000 
ОООО000000000000000000000 


ООО000000200 = - м 0000000000 u 0000000 
О000000000000 и 000000100002000700700 
ОДООО000Д000000000000000000000000000 
ООО000000000000 rn 000 7 = м x n 0000000 
ООО00000000000 > Or ООО00000000000000 
ООООО00000000000000000 и =n x r= n x 
(их Π)ΗΠΗΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΞΠΗΠΠΠΠΠΠΠΗ 
ООО000000000000000000000000000000000 
О00000000000000000000000000000030000 
ОО0000000000000000000000000000000 A L 
ОО0000000000000 70 u On 00000000 


Те Wr ng 0 
R= Гу wy ny O 

по 43. n. Ἡ 

0 0 0 1 


евини 


10.1.6 ΠΠΠΠΠΠ 


ОДОДОДОД000000006 СЕРМ/СООООООО00О0000000 
ΕΕ ΕΕ ΕΕ 


10.2 [I 


ППППРУОреп 5 ΠΠΠΠΠΠΡγ΄ποη Орепо 000 
UUUUnumpy Dp EOD BEIDE TE TCI] 


10.3 [ICD 
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О00000000000000Рагісіеѕуѕёет0000000 
ж 


10.3.1 000000000 


ΠΠΠΠΠΠΗΠΠΠΠΠΠΗΠΝΒΟΓΠΗΠΠΠΠΠΗΠΗΠΗΠΠΠΗ 
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4 create Vertex Array Object (МАО) 

self.vao = glGenVertexArrays(1) 

# bind VAO 


glBindVertexArray(self.vao) 
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= 0.2 
θ quadV = [ 
-S, S, 0.0, 
-S, -5, 0.0, 
s, s, 0.0, 
s, -S, 0.0, 
s, s, 0.0, 
-S, -S, 0.0 
] 
@ vertexData = numpy.array(numP*quadV, numpy.float32) 
Ө self.vertexBuffer = glGenBuffers(1) 
Ө  glBindBuffer(GL ARRAY BUFFER, self.vertexBuffer) 
e 
glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexDat 
a, 
GL_STATIC_DRAW) 
# texture coordinates 
о quadT = | 
0.0, 1.0, 
0.0, 0.0, 


# vertices 


1.0, 1.0, 


1.0, 0.0, 

1.0, 1.0, 

0.0, 0.0 

] 
tcData = numpy.array(numP*quadT, numpy.float32) 
self.tcBuffer = glGenBuffers(1) 
glBindBuffer(GL ARRAY BUFFER, self.tcBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(tcData), tcData, GL STA 
TIC DRAW) 
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4 time lags 
e 
timeData = numpy.repeat(0.005*numpy.arange(numP, dtype=nump 
y.float32), 
4) 
self.timeBuffer = glGenBuffers(1) 


glBindBuffer(GL ARRAY BUFFER, self.timeBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(timeData), timeData, 


GL STATIC DRAW) 
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я velocites 
velocities = [] 
# cone angle 
e coneAngle = math.radians(20.0) 
# set up particle velocities 


for i in range(numP): 


# inclination 
e angleRatio = random.random() 
a = angleRatio*coneAngle 
# azimuth 
e t = random.random()*(2.0*math.pi) 
# get velocity on sphere 
о ух = math.sin(a) *math.cos(t) 
vy = math.sin(a)*math.sin(t) 
νζ = math.cos(a) 
# speed decreases with angle 
e speed = 15.0*(1.0 - angleRatio*angleRatio) 
# add a set of calculated velocities 
о velocities += 6*[speed*vx, speed*vy, speed*vz] 
# set up velocity vertex buffer 
self.velBuffer - glGenBuffers(1) 
glBindBuf fer (GL ARRAY BUFFER, self.velBuffer) 
ο velData = numpy.array(velocities, numpy.float32) 


glBufferData(GL ARRAY BUFFER, 4*len(velData), velData, GL S 
TATIC DRAW) 
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#version 330 core 


in vec3 aVel; 
in vec3 aVert; 
in float aTime0; 


in vec2 aTexCoord; 


uniform mat4 uMVMatrix; 
uniform mat4 uPMatrix; 


uniform mat4 bMatrix; 


uniform float uTime; 
uniform float uLifeTime; 
uniform vec4 uColor; 


uniform vec3 uPos; 


out vec4 vCol; 


out vec2 vTexCoord; 


ООО00000Д0000000000000000000000000000 
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void main() { 
// set position 
Ө float dt = uTime - aTime0; 
O float alpha = clamp(1.0 - 2.0*dt/uLifeTime, 0.0, 1.0); 
© if(dt < 0.0 || dt > uLifeTime || alpha < 0.01) { 
// out of sight! 
gl Position = vec4(0.0, 0.0, -1000.0, 1.0); 
} 
else { 
// calculate new position 


о vec3 accel = vec3(0.0, 0.0, -9.8); 


// apply a twist 
float PI = 3.14159265358979323846264; 


e 
float theta = mod(100.0*length(aVel)*dt, 360.0)*PI/180.0; 


о 
mat4 rot = mat4(vec4(cos(theta), sin(theta), 0.0, 0.0), 


vec4(-sin(theta), cos(theta), 0.0, 0.0), 
vec4(0.0, 0.0, 1.0, 0.0), 
vec4(0.0, 0.0, 0.0, 1.0)); 
// apply billboard matrix 
ο vec4 pos2 = bMatrix*rot*vec4(aVert, 1.0); 
// calculate position 


Ө 
vec3 newPos = pos2.xyz + uPos + aVel*dt + 0.5*accel*dt*dt; 


// apply transformations 


© 
gl Position = uPMatrix * uMVMatrix * vec4(newPos, 1.0); 


} 
// set color 

€ vCol = vec4(uColor.rgb, alpha); 
// set texture coordinates 


vTexCoord = aTexCoord; 
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е (5100сіапарОООДД00000100 
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cos(0) αἰπίθ) 0.0 0.0 

—sin(0) cos{@) 0.0 0.0 
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#version 330 core 


uniform sampler2D uSampler; 
in vec4 vCol; 
in vec2 vTexCoord; 


out vec4 fragColor; 


void main() { 


// get the texture color 


vec4 texCol = texture2D(uSampler, vec2(vTexCoord.s, vTexCoo 
rd.t)); 


// multiply texture color by set vertex color; use the v 
ertex color alpha 


O fragColor = vec4(texCol.rgb*vCol.rgb, vCol.a); 
} 
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100000/00000 

2000000000000 
ЗОДОДОД0000000000000000000 
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М = camera.eye - camera.center 
e N /= numpy.linalg.norm(N) 


U = camera.up 


U /= numpy.linalg.norm(U) 
R = numpy.cross(U, N) 
U2 = numpy.cross(N, R) 
о bMatrix = numpy.array([R[0], 02191, МІ91, 0.0, 
R[1], U2[1], N[1], 0.0, 


R[2], U2[2], N[2], 9.6, 
0.0, 0.0, 0.0, 1.0], numpy.float32) 


e 
glUniformMatrix4fv(self.bMatrixU, 1, GL TRUE, bMatrix) 
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О00000000000000000000000000000000ө 
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4 enable texture 
e glActiveTexture(GL TEXTUREO) 
о glBindTexture(GL_TEXTURE 20, self.texid) 


e glUniformli(self.samplerU, 0) 


Ө 


4 turn depth mask off 
if self.disableDepthMask: 


glDepthMask(GL_FALSE) 


# enable blending 
if self.enableBlend: 
glBlendFunc(GL SRC ALPHA, GL ONE) 


glEnable(GL BLEND) 


# bind VAO 
glBindVertexArray(self.vao) 
# draw 


glDrawArrays(GL TRIANGLES, 0, 6*self.numP) 
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glutils.loadTexture()[ 00000000 
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# a simple camera class 
class Camera: 
"""helper class for viewing""" 
Ө def init (self, eye, center, up): 
self.r - 10.0 
self.theta - 0 
self.eye = numpy.array(eye, numpy.float32) 
self.center - numpy.array(center, numpy.float32) 


self.up = numpy.array(up, numpy.float32) 


def rotate(self): 


"""rotate eye by one step""" 


e self.theta = (self.theta + 1) % 360 
4 recalculate eye 
e self.eye = self.center + numpy.array([ 
self.r*math.cos(math.radians(self.theta)), 
self.r*math.sin(math.radians(self.theta)), 


0.0], numpy.float32) 
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ПООООООООООООООО 
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er/ particle-system/ps.py[(]LLIL] 
import sys, random, math 


import OpenGL 


from OpenGL.GL import * 


import numpy 


import glutils 


strVS = 


#version 330 core 


in vec3 


in vec3 


aVel; 


aVert; 


in float aTime0; 


in vec2 


uniform 
uniform 
uniform 
uniform 
uniform 
uniform 


uniform 


aTexCoord; 


mat4 uMVMatrix; 
mat4 uPMatrix; 
mat4 bMatrix; 
float uTime; 
float uLifeTime; 
vec4 uColor; 


vec3 uPos; 


out vec4 vCol; 


out vec2 vTexCoord; 


void main() { 
// set position 
float dt = uTime - aTime0; 
float alpha = clamp(1.0 - 2.0*dt/uLifeTime, 0.0, 1.0); 
if(dt < 0.0 || dt > uLifeTime || alpha < 0.01) { 
// out of sight! 
gl Position = vec4(0.0, 0.0, -1000.0, 1.0); 
} 
else { 
// calculate new position 
vec3 accel = vec3(0.0, 0.0, -9.8); 
// apply a twist 


float PI = 3.14159265358979323846264; 
float theta = mod(100.0*length(aVel)*dt, 360.0)*PI/180.0; 


mat4 rot = mat4(vec4(cos(theta), sin(theta), 0.0, 0.0), 
vec4(-sin(theta), cos(theta), 0.0, 0.0), 
vec4(0.0, 0.0, 1.0, 0.0), 
vec4(0.0, 0.0, 0.0, 1.0)); 
// apply billboard matrix 
vec4 pos2 - bMatrix*rot*vec4(aVert, 1.0); 


// calculate position 


vec3 newPos = pos2.xyz + uPos + aVel*dt + 0.5*accel*dt*dt; 


// apply transformations 


gl Position = uPMatrix * uMVMatrix * vec4(newPos, 1.0); 
} 
// set color 
vCol = vec4(uColor.rgb, alpha); 
// set tex coords 
vTexCoord = aTexCoord; 


} 


strFS = nuu 


Zversion 330 core 


uniform sampler2D uSampler; 
in vec4 vCol; 
in vec2 vTexCoord; 


out vec4 fragColor; 


void main() { 


// get texture color 


vec4 texCol = texture(uSampler, vec2(vTexCoord.s, vTexCoord 
.t)); 


// multiply by set vertex color; use the vertex color al 
pha 


fragColor = vec4(texCol.rgb*vCol.rgb, vCol.a); 


} 


# a simple camera class 
class Camera: 
"""helper class for viewing""" 
def init (self, eye, center, up): 
self.r = 10.0 
self.theta = 0 
self.eye = numpy.array(eye, numpy.float32) 
self.center - numpy.array(center, numpy.float32) 


self.up = numpy.array(up, numpy.float32) 


def rotate(self): 
"""rotate eye by one step""" 
self.theta = (self.theta + 1) % 360 
4 recalculate eye 
self.eye = self.center + numpy.array([ 


self.r*math.cos(math.radians(self.theta)), 


self.r*math.sin(math.radians(self.theta)), 


0.0], numpy.float32) 


4 particle system class 


class ParticleSystem: 


# initialization 

def init (self, numP): 
4 number of particles 
self.numP - numP 
# time variable 
self.t - 0.0 


self.lifeTime = 5.0 


self.startPos - numpy.array([0.0, 0.0, 0.5]) 
# load texture 

self.texid = glutils.loadTexture('star.png') 
# create shader 


self.program = glutils.loadShaders(strVS, strFS) 


glUseProgram(self.program) 


# set sampler 


texLoc = glGetUniformLocation(self.program, b"uTex") 


glUniformli(texLoc, 0) 


# uniforms 


self.timeU = glGetUniformLocation(self.program, b"uTime") 


self.lifeTimeU = glGetUniformLocation(self.program, b"uLife 
Time") 


self.pMatrixUniform - glGetUniformLocation(self.program, b' 
uPMatrix') 


self.mvMatrixUniform = glGetUniformLocation(self.program, b 
"uMVMatrix") 


self.bMatrixU = glGetUniformLocation(self.program, b"bMatri 
х") 


self.colorU = glGetUniformLocation(self.program, b"uColor") 


self.samplerU = glGetUniformLocation(self.program, b"uSampl 
er") 


self.posU = glGetUniformLocation(self.program, b"uPos") 


# attributes 


self.vertIndex = glGetAttribLocation(self.program, b"aVert" 


self.texIndex = glGetAttribLocation(self.program, b"aTexCoo 


га") 


self.timeOIndex = glGetAttribLocation(self.program, b"aTime 
0") 


self.velIndex = glGetAttribLocation(self.program, b"aVel") 


# render flags 
self.enableBillboard = True 
self.disableDepthMask = True 


self.enableBlend = True 


# which texture to use 
self.useStarTexture = True 
# restart - first time 


self.restart (numP) 


# step 
def step(self): 
# increment time 


self.t += 0.01 


# restart particle system 


def restart(self, numP): 


4 set number of particles 


self.numP = numP 


# time variables 
self.t = 0.0 


self.lifeTime = 5.0 


# color 
self.colO = numpy.array([random.random(), random. random(), 


random.random(), 1.0]) 


# create Vertex Arrays Object (VAO) 
self.vao = glGenVertexArrays(1) 
# bind VAO 


glBindVertexArray(self.vao) 


# create attribute arrays and vertex buffers: 


# vertices 


vertexData = numpy.array(numP*quadV, numpy.float32) 
self.vertexBuffer - glGenBuffers(1) 


glBindBuf fer (GL ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexData 


, 


GL STATIC DRAW) 


# texture coordinates 
quadT = [ 

0.0, 1.0, 

0.0, 0.0, 

1.0, 1.0, 

1.0, 0.0, 

1.0, 1.0, 

0.0, 0.0 

] 


tcData = numpy.array(numP*quadT, numpy.float32) 


self.tcBuffer = glGenBuffers(1) 
glBindBuffer(GL_ARRAY BUFFER, self.tcBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(tcData), tcData, GL STA 
TIC DRAW) 


# time lags 
timeData = numpy.repeat(0.005*numpy.arange(numP, dtype=nump 
y.float32), 
4) 
self.timeBuffer = glGenBuffers(1) 


glBindBuffer(GL_ ARRAY BUFFER, self.timeBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(timeData), timeData, 
GL STATIC DRAW) 
# velocites 
velocities - [] 
я cone angle 
coneAngle = math.radians(20.0) 
# set up particle velocities 
for i in range(numP) : 
# inclination 
angleRatio = random. random ( ) 


a = angleRatio*coneAngle 


# azimuth 

t = random. random()*(2.0*math.pi) 

# get veocity on sphere 

ух = math.sin(a) *math.cos(t) 

vy = math.sin(a)*math.sin(t) 

vz = math.cos(a) 

# speed decreases with angle 

speed = 15.0*(1.0 - angleRatio*angleRatio) 

# add a set of calculated velocities 

velocities += 6*[speed*vx, speed*vy, speed*vz] 
4 set up velocity vertex buffer 
self.velBuffer - glGenBuffers(1) 
glBindBuf fer (GL ARRAY BUFFER, self.velBuffer) 
velData = numpy.array(velocities, numpy.float32) 


glBufferData(GL ARRAY BUFFER, 4*len(velData), velData, GL 5 
TATIC DRAW) 


4 enable arrays 
glEnableVertexAttribArray(self.vertIndex) 
glEnableVertexAttribArray(self.texIndex) 
glEnableVertexAttribArray(self.timeOIndex) 


glEnableVertexAttribArray(self.vellndex) 


4 set buffers 


glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


glVertexAttribPointer(self.vertIndex, 3, GL FLOAT, GL FALSE 
, 0, None) 


glBindBuf fer (GL ARRAY BUFFER, self.tcBuffer) 


glVertexAttribPointer(self.texIndex, 2, GL FLOAT, GL FALSE, 
0, None) 


glBindBuffer(GL ARRAY BUFFER, self.velBuffer) 


glVertexAttribPointer(self.velIndex, 3, GL FLOAT, GL FALSE, 
0, None) 


glBindBuffer(GL ARRAY BUFFER, self.timeBuffer) 


glVertexAttribPointer(self.timeOIndex, 1, GL FLOAT, GL FALS 
E, 0, None) 


4 unbind VAO 


glBindVertexArray (0) 


# render the particle system 


def render(self, pMatrix, mvMatrix, camera): 


# use shader 


glUseProgram(self.program) 


# set projection matrix 


glUniformMatrix4fv(self.pMatrixUniform, 1, GL FALSE, pMatri 
x) 


# set modelview matrix 


glUniformMatrix4fv(self.mvMatrixUniform, 1, GL FALSE, mvMat 
rix) 


# set up a billboard matrix to keep quad aligned to view di 
rection 
if self.enableBillboard: 

М = camera.eye - camera.center 

N /= numpy.linalg.norm(N) 

U = camera.up 

U /= numpy.linalg.norm(U) 

R = numpy.cross(U, N) 

02 = numpy.cross(N, В) 

bMatrix = numpy.array([R[0], U2[0], N[0], 0.0, 

R[1], U2[1], N[1], 0.0, 


R[2], U2[2], N[2], 0.0, 


0.0, 0.0, 0.0, 1.0], numpy.float32) 


glUniformMatrix4fv(self.bMatrixU, 1, GL TRUE, bMatrix) 
else: 
4 identity matrix 
bMatrix = numpy.array([1.0, 0.0, 0.0, 0.0, 
0.0, 1.0, 0.0, 0.6, 


0.0, 0.0, 1.0, 0.0, 


0.0, 0.0, 0.0, 1.0], numpy.float32) 


glUniformMatrix4fv(self.bMatrixU, 1, GL FALSE, bMatrix) 
# set start position 
glUniform3fv(self.posU, 1, self.startPos) 
# set time 
glUniformlf(self.timeU, self.t) 
#set lifetime 
glUniformlf(self.lifeTimeU, self.lifeTime) 
# set color 


glUniform4fv(self.colorU, 1, self.col0) 


4 enable texture 
glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 2D, self.texid) 


glUniformli(self.samplerU, 0) 


# turn depth mask off 
if self.disableDepthMask: 


glDepthMask(GL FALSE) 


# enable blending 
if self.enableBlend: 
glBlendFunc(GL SRC ALPHA, GL ONE) 


glEnable(GL BLEND) 


# bind VAO 

glBindVertexArray(self.vao) 

# draw 

glDrawArrays(GL TRIANGLES, 0, 6*self.numP) 
# unbind VAO 


glBindVertexArray (0) 


# disable blend 
if self.enableBlend: 


glDisable(GL BLEND) 


# turn depth mask on 
if self.disableDepthMask: 


glDepthMask(GL_TRUE) 


# disable texture 


glBindTexture(GL TEXTURE 2D, 0) 


TO و‎ 


10.5 [I 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


import sys, random, math 
import OpenGL 
from OpenGL.GL import * 
import numpy 


import glutils 


st rVS = nuu 


Zversion 330 core 


in vec3 aVert; 
uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 


out vec4 vCol; 


void main() { 
// apply transformations 
gl Position = uPMatrix * uMVMatrix * vec4(aVert, 
// set color 
vCol = vec4(0.8, 0.0, 0.0, 1.0); 


} 


strFS = пп 


#version 330 core 


in vec4 vCol; 


out vec4 fragColor; 


void main() { 
// use vertex color 
fragColor = vCol; 

} 


1.0); 


class Box: 
def init (self, side): 


self.side = side 


# load shaders 
self.program = glutils.loadShaders(strVS, strFS) 


glUseProgram(self.program) 


S = side/2.0 


vertices - [ 


5, -5,-5, 
5, 5, 5, 
5, 5, -5, 
5, 5, 5, 
5, -5,-5 


4 set up vertex array object (МАО) 

self.vao = glGenVertexArrays(1) 
glBindVertexArray(self.vao) 

# set up VBOs 

vertexData = numpy.array(vertices, numpy.float32) 
self.vertexBuffer = glGenBuffers(1) 


glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexData 


, 


GL STATIC DRAW) 


#enable arrays 


self.vertIndex = glGetAttribLocation(self.program, "aVert") 


glEnableVertexAttribArray(self.vertIndex) 


4 set buffers 


glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


glVertexAttribPointer(self.vertIndex, 3, GL FLOAT, GL FALSE 


, 0, None) 
# unbind VAO 


glBindVertexArray(0) 


def render(self, pMatrix, mvMatrix): 


# use shader 


glUseProgram(self.program) 


# set projection matrix 
glUniformMatrix4fv(glGetUniformLocation(self.program, 
{г1х'), 

1, GL FALSE, pMatrix) 

# set modelview matrix 
glUniformMatrix4fv(glGetUniformLocation(self.program, 
atrix'), 

1, GL FALSE, mvMatrix) 

# bind VAO 

glBindVertexArray(self.vao) 


# draw 


' uPMa 


‘uMVM 


glDrawArrays(GL_TRIANGLES, 0, 36) 
# unbind VAO 


glBindVertexArray(0) 
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class PSMaker: 
"""GLFW Rendering window class for Particle System""" 
def init (self): 
6 self.camera - Camera([15.0, 0.0, 2.5], 
[0.0, 0.0, 2.5], 
[0.0, 0.0, 1.0]) 
self.aspect - 1.0 
self.numP - 300 
self.t- 0 
# flag to rotate camera view 


self.rotate - True 


# save current working directory 


сыа = os.getcwd() 


# initialize glfw; this changes cwd 


e glfw.glfwlnit() 


# restore cwd 


os.chdir(cwd) 


# version hints 
gl fw.gl fwWindowHint(glfw.GLFW CONTEXT VERSION MAJOR, 3) 
gl fw.gl fwWindowHint(glfw.GLFW CONTEXT VERSION MINOR, 3) 


gl fw. gl fwWindowHint(glfw.GLFW_OPENGL FORWARD COMPAT, GL TRU 
E) 


gl fw.gl fwWindowHint (glfw.GLFW OPENGL PROFILE, 


gl fw.GLFW OPENGL CORE PROFILE) 


# make a window 
self.width, self.height = 640, 480 


self.aspect = self.width/float(self.height) 


self.win = glfw.glfwCreateWindow(self.width, self.height, 


b"Particle System") 


# make context current 


gl fw.glfwMakeContextCurrent(self.win) 


# initialize GL 
glViewport(0, 0, self.width, self.height) 
glEnable(GL DEPTH TEST) 


glClearColor(0.2, 0.2, 0.2,1.0) 

# set window callbacks 
glfw.glfwSetMouseButtonCallback(self.win, self.onMouseButto 
n) 

glfw.glfwSetKeyCallback(self.win, self.onKeyboard) 


glfw.glfwSetWindowSizeCallback(self.win, self.onSize) 


8 create 3D 
self.psys - ParticleSystem(self.numP) 


self.box = Box(1.0) 


# exit flag 


о self.exitNow = False 
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ѕер()0000000000000000 
def step(self): 

# increment time 

e self.t += 10 

о self.psys.step() 
# rotate eye 
if self.rotate: 

© self.camera.rotate() 
4 restart every 5 seconds 
if not int(self.t) % 5000: 


о self.psys.restart(self.numP) 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΘΠΠΠΠ 
ParticleSystem[step()000000000000000000 


ни 


10.6.2 [II 
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€ def onKeyboard(self, win, key, scancode, action, mods): 
#print 'keyboard: ', win, key, scancode, action, mods 
if action == glfw.GLFW PRESS: 
# ESC to quit 
if key == glfw.GLFW KEY ESCAPE: 
self.exitNow - True 
elif key -- glfw.GLFW KEY R: 
self.rotate - not self.rotate 
elif key == glfw.GLFW KEY B: 


# toggle billboarding 


self.psys.enableBillboard - not self.psys.enableBillboard 
elif key == glfw.GLFW KEY D: 


# toggle depth mask 


self.psys.disableDepthMask = not self.psys.disableDepthMask 
elif key == glfw.GLFW KEY T: 


# toggle transparency 


self.psys.enableBlend = not self.psys.enableBlend 


еПрдородрродрооороророоодррдрророороодр 
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10.6.3 00000 
ILIGLFWIIIIIILI IILI IILI III IILI IILI II 


def run(self): 
# initializer timer 
glfw.SetTime(0) 
t = 0.0 
о 
while not glfw.glfwWindowShouldClose(self.win) and not sel 
f.exitNow: 
# update every x seconds 
6 currT = glfw.glfwGetTime() 
if currT - t » 0.01: 
4 update time 


t = currT 


4 clear 


glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


4 render 


pMatrix = glutils.perspective(100.0, self.aspect, 0.1, 100. 
0) 


# modelview matrix 
mvMatrix = glutils.lookAt(self.camera.eye, self.camera.cent 
er, 


self.camera.up) 


# draw nontransparent object first 


© self.box.render(pMatrix, mvMatrix) 


# render 


о 
self.psys.render(pMatrix, mvMatrix, self.camera) 


# step 


e self.step() 


glfw.glfwSwapBuffers(self.win) 
# poll for and process events 
glfw.glfwPollEvents() 

# end 


glfw.glfwTerminate() 
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import sys, os, math, numpy 
import OpenGL 
from OpenGL.GL import * 
import numpy 
from ps import ParticleSystem, Camera 
from box import Box 
import glutils 


import glfw 


class PSMaker: 
"""GLFW Rendering window class for Particle System""" 


def init (self): 


self.camera = Сатега([15.0, 0.0, 2.5], 
[0.0, 0.0, 2.5], 
[0.0, 0.0, 1.0]) 


1.0 


self.aspect 
self.numP - 300 

self.t = 0 

# flag to rotate camera view 


self.rotate = True 


# save current working directory 


сма = os.getcwd() 


# initialize glfw; this changes сма 


glfw.glfwInit() 


4 restore cwd 


os.chdir(cwd) 


# version hints 


glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MAJOR, 3) 


glfw.glfwwindowHint(glfw.GLFW CONTEXT VERSION MINOR, 3) 


gl fw.gl fwWindowHint(glfw.GLFW OPENGL FORWARD COMPAT, GL TRU 
E) 


gl fw.gl fwWindowHint (glfw.GLFW OPENGL PROFILE, 


gl fw.GLFW OPENGL CORE PROFILE) 


# make a window 
self.width, self.height = 640, 480 


self.aspect = self.width/float(self.height) 


self.win = glfw.glfwCreateWindow(self.width, self.height, 
b"Particle System") 
4 make context current 


glfw.glfwMakeContextCurrent(self.win) 


# initialize GL 
glViewport(0, 0, self.width, self.height) 
glEnable(GL DEPTH TEST) 


glClearColor(0.2, 0.2, 0.2,1.0) 


4 set window callbacks 


glfw.glfwSetMouseButtonCallback(self.win, self.onMouseButto 
n) 


glfw.glfwSetKeyCallback(self.win, self.onKeyboard) 


glfw.glfwSetWindowSizeCallback(self.win, self.onSize) 


# create 3D 


self.psys = ParticleSystem(self.numP) 


self.box = Box(1.0) 


# exit flag 


self.exitNow = False 


def onMouseButton(self, win, button, action, mods): 


#print 'mouse button: ', win, button, action, mods 


pass 


def onKeyboard(self, win, key, scancode, action, mods): 
#print 'keyboard: ', win, key, scancode, action, mods 
if action == glfw.GLFW PRESS: 
# ESC to quit 
if key == glfw.GLFW KEY ESCAPE: 
self.exitNow - True 
elif key == glfw.GLFW KEY R: 
self.rotate - not self.rotate 


elif key == glfw.GLFW KEY B: 


# toggle billboarding 


self.psys.enableBillboard = not self.psys.enableBillboard 
elif key == glfw.GLFW KEY D: 


# toggle depth mask 


self.psys.disableDepthMask = not self.psys.disableDepthMask 
elif key == glfw.GLFW KEY T: 


# toggle transparency 


self.psys.enableBlend = not self.psys.enableBlend 


def onSize(self, win, width, height): 
#print 'onsize: ', win, width, height 
self.width = width 


self.height height 


self.aspect = width/float (height) 


glViewport(0, 0, self.width, self.height) 


def step(self): 
# increment time 
self.t += 10 
self.psys.step() 


# rotate eye 


if self.rotate: 
self.camera. rotate() 

# restart every 5 seconds 

if not int(self.t) % 5000: 


self.psys.restart(self.numP) 


def run(self): 
# initializer timer 
glfw.glfwSetTime(0) 
t = 0.0 
while not glfw.glfwWindowShouldClose(self.win) and not self 
.exitNow: 
# update every x seconds 
currT = glfw.glfwGetTime() 
if currT - t > 0.01: 
# update time 


t = currT 


4 clear 


glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


# render 


pMatrix = glutils.perspective(100.0, self.aspect, 0.1, 100. 
0) 


# modelview matrix 
mvMatrix = glutils.lookAt(self.camera.eye, self.camera.cent 
er, 


self.camera.up) 


# draw nontransparent object first 


self.box.render(pMatrix, mvMatrix) 


# render 


self.psys.render(pMatrix, mvMatrix, self.camera) 


# step 


self.step() 


gl fw.glfwSwapBuf fers (self .win) 
# poll for and process events 
glfw.glfwPollEvents() 

# end 


glfw.glfwTerminate() 


4 main() function 

def main(): 
# use sys.argv if needed 
print('starting particle system...') 
prog = PSMaker() 


prog. гип() 


# call main 
if name == ' main ': 


main() 
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5 python3 psmain.py 
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https://github.com/electronut/pp/tree/mast 
er/common/(/[]LL) 


11.4 00000 


ОООД0000000000000000000000009000 
volreader. py 000000000001 1.500 


def loadVolume(dirName): 


"""read volume from directory as a 3D texture""" 


# list images in directory 

Ө files = sorted(os.listdir(dirName) ) 
print('loading images from: %s' % dirName) 
imgDataList = [] 
count = 0 
width, height = 0, 0 
for file in files: 


6 
file path = os.path.abspath(os.path.join(dirName, file)) 


try: 
# read image 
© img = Image.open(file path) 


imgData = np.array(img.getdata(), np.uint8) 


# check if all images are of the same size 
о if count is 0: 
width, height = img.size[0], img.size[1] 
imgDataList.append(imgData) 


else: 


e 
if (width, height) == (img.size[0], img.size[1]): 
imgDataList.append(imgData) 
else: 


print('mismatch') 


raise RunTimeError("image size mismatch") 
count += 1 
#print img.size 
except: 
# skip 


print('Invalid image: %5' % file path) 


# load image data into single array 
depth = count 
Ө data = np.concatenate(imgDataList) 


print('volume data dims: %d %d %d' % (width, height, depth) 
) 


# load data into 3D texture 

@ texture = glGenTextures(1) 
glPixelStorei(GL UNPACK ALIGNMENT, 1) 
glBindTexture(GL TEXTURE 3D, texture) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP S, GL CLAMP . 
TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP T, GL CLAMP ` 
TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP В, GL CLAMP ` 


TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE MAG FILTER, GL LI 
NEAR) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE MIN FILTER, GL LI 
NEAR) 


O glTexlImage3D(GL TEXTURE 3D, 0, GL RED, 
width, height, depth, 0, 
GL RED, GL UNSIGNED BYTE, data) 
# return texture 


O return (texture, width, height, depth) 


loadVolume()DDmnmiesrttistdie 000000000 
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ООООД0000000000000000500000000 
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11.5 ΠΠΠΠΠΠΠΠΠ 


00000000000000 
https://github.com/electronut/pp/tree/mast 
er/volrender/[][]volreader.py[]L]L] 

import os 

import numpy as np 


from PIL import Image 


import OpenGL 


from OpenGL.GL import * 
from scipy import misc 


def loadVolume(dirName): 


"""read volume from directory as a 3D texture""" 


# list images in directory 

files = sorted(os.listdir(dirName) ) 
print('loading images from: %' % dirName) 
imgDataList = [] 

count = 0 

width, height = 0, 0 


for file in files: 


file path = os.path.abspath(os.path.join(dirName, file) ) 
try: 
# read image 
img = Image.open(file path) 


imgData = np.array(img.getdata(), np.uint8) 


# check if all are of the same size 

if count is 0: 
width, height = img.size[0], img.size[1] 
imgDataList.append(imgData) 


else: 


if (width, height) == (img.size[0], img.size[1]): 
imgDataList.append(imgData) 


else: 


print('mismatch') 
raise RunTimeError("image size mismatch") 
count += 1 
#print img.size 
except: 
4 skip 


print('Invalid image: %s' % file path) 


# load image data into single array 
depth = count 
data = np.concatenate(imgDataLlist) 


print('volume data dims: %d %d %d' % (width, height, depth) 
) 


# load data into 3D texture 

texture = glGenTextures(1) 
glPixelStorei(GL UNPACK ALIGNMENT, 1) 
glBindTexture(GL TEXTURE 3D, texture) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP S, GL CLAMP . 
TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP T, GL CLAMP ` 
TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE WRAP В, GL CLAMP ` 
TO EDGE) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE MAG FILTER, GL LI 
NEAR) 


glTexParameterf(GL TEXTURE 3D, GL TEXTURE MIN FILTER, GL LI 
NEAR) 
glTexlImage3D(GL TEXTURE 3D, 0, GL RED, 
width, height, depth, 0, 
GL RED, GL UNSIGNED BYTE, data) 
#return texture 


return (texture, width, height, depth) 


# load texture 

def loadTexture(filename): 
img = Image.open(filename) 
img data = np.array(list(img.getdata()), 'В') 
texture = glGenTextures(1) 
glPixelStorei(GL UNPACK ALIGNMENT, 1) 
glBindTexture(GL TEXTURE 2D, texture) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE WRAP S, GL CLAMP ` 
TO EDGE) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE WRAP T, GL CLAMP ` 


TO EDGE) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE MAG FILTER, GL LI 
NEAR) 


glTexParameterf(GL TEXTURE 2D, GL TEXTURE MIN FILTER, GL LI 
NEAR) 


glTexlImage2D(GL TEXTURE 2D, ©, GL АСВА, img.size[0], img.si 
ze[1], 
0, GL АСВА, GL UNSIGNED BYTE, img data) 


return texture 
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Ор00000000Каусиберрроророророрбододо 
ОДОООДОДООДОООРВОДОДООООООООДОД0ОООО 
00000000гаусибе. ру 000000011.700 


ООО000000000000 


#version 330 соге 


layout(location = 1) in vec3 cubePos; 


layout(location = 2) in vec3 cubeCol; 


uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 


out vec4 vColor; 


void main() 


{ 


// set back-face color 


vColor = vec4(cubeCol.rgb, 1.0); 


// transformed position 


vec4 newPos = vec4(cubePos.xyz, 1.0); 


// set position 


gl Position = uPMatrix * uMVMatrix * newPos; 


} 


Ө strFS = """ 


#version 330 соге 


in vec4 vColor; 


out vec4 fragColor; 


void main() 


fragColor = vColor; 


ΠΘΠΠΠΠΠΒΘν(ΘαὈΕΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
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# cube vertices 
e vertices - numpy.array([ 
0.0, 0.0, 0.0, 
1.0, 0.0, 0.0, 
1.0, 1.0, 0.0, 
0.0, 1.0, 0.0, 
0.0, 0.0, 1.0, 
1.0, 0.0, 1.0, 
1.0, 1.0, 1.0, 


0.0, 1.0, 1.0 


], numpy.float32) 


# cube colors 
e colors = numpy.array([ 
0.0, 0.0, 0.0, 
1.0, 0.0, 0.0, 


1.0, 1.0, 0.0, 


1.0, 1.0, 1.0, 
0.0, 1.0, 1.0 


], numpy.float32) 


# individual triangles 


e indices = numpy.array([ 
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ООО00000000000000 


# set up vertex array object (VAO) 


self.vao = glGenVertexArrays(1) 


glBindVertexArray(self.vao) 


# vertex buffer 
self.vertexBuffer = glGenBuffers(1) 
glBindBuffer(GL_ ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(vertices), vertices, GL 
_ STATIC DRAW) 


# vertex buffer — cube vertex colors 
self.colorBuffer = glGenBuffers(1) 
glBindBuffer(GL_ARRAY BUFFER, self.colorBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(colors), colors, GL STA 
TIC DRAW) 


# index buffer 

self.indexBuffer = glGenBuffers(1) 
e 
glBindBuffer(GL ELEMENT ARRAY BUFFER, self.indexBuffer); 
glBufferData(GL ELEMENT ARRAY BUFFER, 2*len(indices), indic 
es, 


GL STATIC DRAW) 


ОО000000000000000000000000000УАОПО00 


ПО000000000000000000өП0іпаісеѕ00000 
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def initFBO(self): 
# create frame buffer object 
self.fboHandle = glGenFramebuf fers (1) 
# create texture 
self.texHandle = glGenTextures(1) 
# create depth buffer 


self.depthHandle = glGenRenderbuf fers (1) 


# bind 


glBindFramebuffer(GL_FRAMEBUFFER, self. fboHandle) 


glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 20, self.texHandle) 


# set parameters to draw the image at different sizes 


e 
glTexParameteri(GL TEXTURE 2D, GL TEXTURE MIN FILTER, GL LI 


NEAR) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE MAG FILTER, GL LI 
NEAR) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE WRAP S, GL CLAMP ` 
TO EDGE) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE WRAP T, GL CLAMP _ 
TO EDGE) 

# set up texture 
glTexlImage2D(GL TEXTURE 2D, ©, GL ВОВА, self.width, self.he 
ight, 


0, GL RGBA, GL UNSIGNED BYTE, None) 


4 bind texture to FBO 


o 
glFramebufferTexture2D(GL FRAMEBUFFER, GL COLOR ATTACHMENTO 


, 


GL TEXTURE 2D, self.texHandle, 0) 


# bind 


e glBindRenderbuffer(GL_RENDERBUFFER, self.depthHandle) 


glRenderbufferStorage(GL RENDERBUFFER, GL DEPTH COMPONENT24 


, 


self.width, self.height) 


# bind depth buffer to FBO 


glFramebuf ferRenderbuffer(GL FRAMEBUFFER, GL DEPTH ATTACHME 
NT, 


GL RENDERBUFFER, self.depthHandle) 
# check status 
о status = glCheckFramebufferStatus (GL FRAMEBUFFER) 
if status == GL FRAMEBUFFER COMPLETE: 
pass 
#print "Тро %d complete" % self. fboHandle 
elif status == GL FRAMEBUFFER UNSUPPORTED: 
print "fbo %d unsupported" % self. fboHandle 
else: 


print "fbo за Error" 5 self.fboHandle 


О00000000000000000000000000000ө0000 

оррробррерррородоводбдоороФеррооооооро000 
О000000024000000000000000000006өП000 

ООО0О0000000000000000000000000000000 
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11.6.3 ΠΠΠΠΠΠΠΠ 
ООО00000000000000 


def renderBackFace(self, pMatrix, mvMatrix): 


"""renders back-face of гау- 
cube to a texture and returns ії""" 


# render to FBO 
e glBindFramebuffer(GL FRAMEBUFFER, self.fboHandle) 
# set active texture 
glActiveTexture(GL TEXTUREO) 
# bind to FBO texture 


glBindTexture(GL_TEXTURE 20, self.texHandle) 


# render cube with face culling enabled 


e 
self.renderCube(pMatrix, mvMatrix, self.program, True) 


# unbind texture 
e glBindTexture(GL TEXTURE 2D, 0) 
glBindFramebuffer(GL FRAMEBUFFER, 0) 


glBindRenderbuffer(GL_RENDERBUFFER, 0) 


# return texture ID 


ο return self.texHandle 
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11.6.4 ПШПШ 


ΕΕ ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 
T ο μμ. 
alse 


def renderFrontFace(self, pMatrix, mvMatrix, program): 
"""render front-face of гау-сире""" 
# no face culling 


self.renderCube(pMatrix, mvMatrix, program, False) 


11.6.5 [HII 
О0000гепаегсире()000000000000000000 


def renderCube(self, pMatrix, mvMatrix, program, cullFace): 


"""renderCube uses face culling if flag set""" 
glClear(GL_COLOR BUFFER BIT | GL DEPTH BUFFER ВІТ) 


# set shader program 


glUseProgram(program) 


# set projection matrix 


glUniformMatrix4fv(glGetUniformLocation(program, b'uPMatrix 
ET 


1, GL FALSE, pMatrix) 


4 set modelview matrix 
glUniformMatrix4fv(glGetUniformLocation(program, b'uMVMatri 
x'), 


1, GL FALSE, mvMatrix) 


# enable face culling 


glDisable(GL CULL FACE) 


e if cullFace: 
glFrontFace(GL CCW) 
glCullFace(GL FRONT) 


glEnable(GL CULL FACE) 


# bind VAO 


glBindVertexArray(self.vao) 


# animated slice 


e 
glDrawElements(GL TRIANGLES, self.nIndices, GL UNSIGNED SHO 


АТ, Мопе) 


# unbind VAO 


glBindVertexArray(0) 


# reset cull face 
if cullFace: 
# disable face culling 


glDisable(GL CULL FACE) 
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ООРВОДДОООООООООД00О00О0000000000000000 
o RayCubeDLDmg nr 


€ def reshape(self, width, height): 
self.width - width 
self.height - height 
self.aspect = width/float(height) 


4 re-create FBO 


self.clearFBO() 


self.initFBO() 


ΠΠΟΡΕΠΟΙ ΠΠΠΠΠΠΠΠΠΠΓ65Π8ρ6ε0ΠΠΘῃΠ 


11.7 (ПП 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΒΕ 
https://github.com/electronut/pp/tree/mast 


er/volrender/[][]raycube.py[][][] 
import OpenGL 
from OpenGL.GL import * 


from OpenGL.GL.shaders import * 


import numpy, math, sys 


import volreader, glutils 


st rVS = nuu 


#version 330 core 


layout(location = 1) in vec3 cubePos; 
layout(location = 2) in vec3 cubeCol; 
uniform mat4 uMVMatrix; 
uniform mat4 uPMatrix; 


out vec4 vColor; 


void та1п() 


{ 


} 


// set back face color 


vColor = vec4(cubeCol.rgb, 1.0); 


// transformed position 


vec4 newPos = vec4(cubePos.xyz, 1.0); 


// set position 


gl Position = uPMatrix * uMVMatrix * newPos; 


strFS = nuu 


Zversion 330 core 


in vec4 vColor; 


out vec4 fragColor; 


void main() 


{ 


fragColor = vColor; 


} 


class RayCube: 


"""class used to generate rays used in ray casting""" 


def init (self, width, height): 


"""RayCube constructor""" 


4 set dims 


self.width, self.height - width, height 


4 create shader 


self.program - glutils.loadShaders(strVS, strFS) 


# cube vertices 


vertices - numpy.array([ 


0.0, 0.0, 0.0, 
1.0, 0.0, 0.0, 
1.0, 1.0, 0.9, 


0.0, 1.0, 0.0, 


0.0, 0.0, 1.0, 


# сире 


colors 


1.0, 0.0, 1.0, 
1.0, 1.0, 1.0, 
0.0, 1.0, 1.0 


], numpy.float32) 


colors 
= numpy.array([ 
0.0, 0.0, 0.0, 


1.0, 0.0, 0.0, 


1.0, 1.0, 1.0, 
0.0, 1.0, 1.0 


1, numpy.float32) 


4 individual triangles 


indices = numpy.array([ 


1, numpy.int16) 


self.nIndices = indices.size 


# set up vertex array object (\/АО) 
self.vao = glGenVertexArrays(1) 


glBindVertexArray(self.vao) 


#vertex buffer 
self.vertexBuffer = glGenBuffers(1) 
glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(vertices), vertices, GL 
_ STATIC DRAW) 


# vertex buffer - cube vertex colors 
self.colorBuffer = glGenBuffers(1) 


glBindBuffer(GL_ARRAY BUFFER, self.colorBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(colors), colors, GL STA 
TIC DRAW); 


# index buffer 


self.indexBuffer = glGenBuffers(1) 
glBindBuffer(GL ELEMENT ARRAY BUFFER, self.indexBuffer) ; 


glBufferData(GL ELEMENT ARRAY BUFFER, 2*len(indices), indic 
es, 


GL STATIC DRAW) 


# enable attrs using the layout indices in shader 
aPosLoc = 1 


aColorLoc = 2 
# bind buffers: 
glEnableVertexAttribArray (1) 


glEnableVertexAttribArray (2) 


# vertex 


glBindBuffer(GL_ARRAY BUFFER, self.vertexBuffer) 


glVertexAttribPointer(aPosLoc, 3, GL FLOAT, GL FALSE, 0, No 
ne) 


# color 

glBindBuffer(GL_ARRAY BUFFER, self.colorBuffer) 
glVertexAttribPointer(aColorLoc, 3, GL FLOAT, GL FALSE, 0, 
None) 


# index 


glBindBuffer(GL ELEMENT ARRAY BUFFER, self.indexBuf fer) 


Я unbind МАО 


glBindVertexArray(0) 


# FBO 


self.initFBO() 


def renderBackFace(self, pMatrix, mvMatrix): 


"""renders back-face of ray- 
cube to a texture and returns it""" 


# render to FBO 


glBindFramebuffer(GL_FRAMEBUFFER, self.fboHandle) 


4 set active texture 
glActiveTexture(GL TEXTUREO) 
# bind to FBO texture 


glBindTexture(GL_TEXTURE 2D, self.texHandle) 


# render cube with face culling enabled 


self.renderCube(pMatrix, mvMatrix, self.program, True) 


# unbind texture 
glBindTexture(GL TEXTURE 2D, 0) 
glBindFramebuffer(GL FRAMEBUFFER, 0) 


glBindRenderbuffer(GL RENDERBUFFER, 0) 


4 return texture ID 


return self.texHandle 


def renderFrontFace(self, pMatrix, mvMatrix, program): 
"""render front face of ray-cube""" 
# no face culling 


self.renderCube(pMatrix, mvMatrix, program, False) 


def renderCube(self, pMatrix, mvMatrix, program, cullFace): 


"""render cube use face culling if flag set""" 


glClear(GL_COLOR BUFFER BIT | GL DEPTH BUFFER ВІТ) 
# set shader program 


glUseProgram(program) 


# set projection matrix 


glUniformMatrix4fv(glGetUniformLocation(program, b'uPMatrix 
'), 


1, GL_FALSE, pMatrix) 


# set modelview matrix 
glUniformMatrix4fv(glGetUniformLocation(program, b'uMVMatri 
x'), 


1, GL FALSE, mvMatrix) 


# enable face culling 

glDisable(GL CULL FACE) 

if cullFace: 
glFrontFace(GL CCW) 
glCullFace(GL FRONT) 


glEnable(GL CULL FACE) 


# bind VAO 


glBindVertexArray(self.vao) 


# animated slice 


glDrawElements(GL TRIANGLES, self.nIndices, GL UNSIGNED SHO 
RT, None) 


# unbind VAO 


glBindVertexArray (0) 


# reset cull face 
if cullFace: 
# disable face culling 


glDisable(GL CULL FACE) 


def reshape(self, width, height): 


self.width - width 


self.height height 
self.aspect = width/float(height) 
# re-create FBO 


self.clearFBO() 


self.initFBO() 


def initFBO(self): 
4 create frame buffer object 
self.fboHandle = glGenFramebuffers(1) 
# create texture 
self.texHandle - glGenTextures(1) 
4 create depth buffer 


self.depthHandle = glGenRenderbuffers(1) 


# bind 


glBindFramebuffer(GL_FRAMEBUFFER, self.fboHandle) 


glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 2D, self.texHandle) 


# set parameters to draw the image at different sizes 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE MIN FILTER, GL LI 
NEAR) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE MAG FILTER, GL LI 
NEAR) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE WRAP S, GL CLAMP ` 


TO EDGE) 


glTexParameteri(GL TEXTURE 2D, GL TEXTURE WRAP T, GL CLAMP ` 
TO EDGE) 


# set up texture 


glTexlImage2D(GL TEXTURE 2D, ©, GL ВОВА, self.width, self.he 
ight, 


0, GL RGBA, GL UNSIGNED BYTE, None) 


4 bind texture to FBO 


glFramebufferTexture2D(GL FRAMEBUFFER, GL COLOR ATTACHMENTO 


, 


GL TEXTURE 2D, self.texHandle, 0) 

# bind 
glBindRenderbuf fer(GL_RENDERBUFFER, self.depthHandle) 
glRenderbufferStorage(GL_RENDERBUFFER, GL DEPTH COMPONENT24 


, 


self.width, self.height) 


# bind depth buffer to FBO 


glFramebuf ferRenderbuffer(GL FRAMEBUFFER, GL DEPTH АТТАСНМЕ 
NT, 


GL RENDERBUFFER, self.depthHandle) 


# check status 
status = glCheckFramebufferStatus(GL FRAMEBUFFER) 
if status == GL FRAMEBUFFER COMPLETE: 
pass 
#print "Тро %d complete" % self. fboHandle 
elif status == GL FRAMEBUFFER UNSUPPORTED: 
print("fbo %d unsupported" % self.fboHandle) 
else: 


print("fbo за Error" 5 self.fboHandle) 


glBindTexture(GL TEXTURE 2D, 0) 
glBindFramebuffer(GL FRAMEBUFFER, 0) 
glBindRenderbuffer(GL_RENDERBUFFER, 0) 


return 


def clearFBO(self): 
"""Clears old FBO""" 
# delete FBO 


if glIsFramebuffer(self.fboHandle): 


glDeleteFramebuffers(int(self.fboHandle)) 


# delete texture 
if glIsTexture(self.texHandle): 


glDeleteTextures(int(self.texHandle) ) 


def close(self): 
"""Call this to free up OpenGL resources""" 
glBindTexture(GL TEXTURE 2D, 0) 
glBindFramebuffer(GL FRAMEBUFFER, 0) 
glBindRenderbuffer(GL RENDERBUFFER, 0) 
# delete FBO 
if glIsFramebuffer(self.fboHandle): 


glDeleteFramebuffers(int(self.fboHandle)) 


4 delete texture 
if glIsTexture(self.texHandle): 


glDeleteTextures(int(self.texHandle) ) 


# delete render buffer 


if glIsRenderbuffer(self.depthHandle): 


glDeleteRenderbuffers(1, int(self.depthHandle)) 
# delete buffers 


glDeleteBuffers(1, self. vertexBuffer) 
glDeleteBuffers(1, & indexBuffer) 


glDeleteBuffers(1, & colorBuffer) 
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def init (self, width, height, volume): 


"""RayCastRender construction""" 


4 create RayCube object 


e self.raycube = raycube.RayCube(width, height) 


# set dimensions 


self.width = width 


self.height = height 


self.aspect = width/float (height) 


# create shader 


e self.program = glutils.loadShaders(strVS, strFS) 
# texture 
e self.texVolume, self.Nx, self.Ny, self.Nz = volume 


# initialize camera 


ο self.camera = Camera() 
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def draw(self): 


# build projection matrix 

e 

pMatrix = glutils.perspective(45.0, self.aspect, 0.1, 100.0 
# modelview matrix 


ө 
mvMatrix = glutils.lookAt(self.camera.eye, self.camera.cent 


ег, 


self.camera.up) 


# render 


# generate ray-cube back-face texture 


© 
texture = self.raycube.renderBackFace(pMatrix, mvMatrix) 


# set shader program 


о glUseProgram(self.program) 


# set window dimensions 


glUniform2f(glGetUniformLocation(self.program, b"uWinDims") 


, 


float(self.width), float(self.height)) 


# bind to texture unit 0, which represents back- 
faces of cube 


e glActiveTexture(GL TEXTUREO) 
glBindTexture(GL TEXTURE 2D, texture) 


glUniformli(glGetUniformLocation(self.program, b"texBackFac 
es"), 0) 


# texture unit 1: 3D volume texture 
Q glActiveTexture(GL TEXTURE1) 


glBindTexture(GL TEXTURE 3D, self.texVolume) 


glUniformli(glGetUniformLocation(self.program, b"texVolume" 
), 1) 


4 draw front-face of cubes 


ο 
self.raycube.renderFrontFace(pMatrix, mvMatrix, self.progra 
m) 
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#version 330 core 


Ө layout(location = 1) in vec3 cubePos; 


layout(location = 2) in vec3 cubeCol; 


Ө uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 
© out vec4 vColor; 


void main() 


{ 


// set position 


о 
gl Position = uPMatrix * uMVMatrix * vec4(cubePos.xyz, 1. 


// set color 
e vColor = vec4(cubeCol.rgb, 1.0); 


} 
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#version 330 core 
in vec4 vColor; 


uniform sampler2D texBackFaces; 
uniform sampler3D texVolume; 


uniform vec2 uWinDims; 
out vec4 fragColor; 


void main() 
{ 
// start of гау 


€ vec3 start = vColor.rgb; 


// calculate texture coordinates at fragment, 


// which is a fraction of window coordinates 


Ө vec2 texc = gl FragCoord.xy/uWinDims. xy; 


// get end of ray by looking up back-face color 


Ө  vec3 end = texture(texBackFaces, texc).rgb; 


// calculate ray direction 


Ө  vec3 dir = end - start; 


// normalized ray direction 


vec3 norm dir - normalize(dir); 
// the length from front to back is calculated and 
// used to terminate the ray 


float len = Length(dir.xyz); 


// ray step size 


float stepSize - 0.01; 


// x-ray projection 


vec4 dst = vec4(0.0); 


// step through the ray 


for(float t = 0.0; t < len; t += stepSize) { 


// set position to end point of ray 


vec3 samplePos = start + t*norm dir; 


// get texture value at position 
float val = texture(texVolume, samplePos).r; 


vec4 src = vec4(val); 


// set opacity 
src.a *= 0.1; 


src.rgb *= src.a; 


// blend with previous value 


dst = (1.0 - dst.a)*src + dst; 


// exit loop when alpha exceeds threshold 
if(dst.a >= 0.95) 


break; 


// set fragment color 


fragColor = dst; 


} 
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er/volrender/[|[]raycast.py[ IL! 
import OpenGL 
from OpenGL.GL import * 


from OpenGL.GL.shaders import * 


import numpy as np 


import math, sys 
import raycube, glutils, volreader 


st rVS = ипи 


#version 330 соге 


layout(location = 1) in vec3 cubePos; 


layout(location = 2) in vec3 cubeCol; 


uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 
out vec4 vColor; 


void main() 


{ 


// set position 


gl Position = uPMatrix * uMVMatrix κ vec4(cubePos.xyz, 1.0) 


, 


// set color 
vColor = vec4(cubeCol.rgb, 1.0); 
} 


strFS = nu 


Zversion 330 core 


in vec4 vColor; 


uniform sampler2D texBackFaces; 


uniform sampler3D texVolume; 


uniform vec2 uWinDims; 


out vec4 fragColor; 


void main() 


{ 


// start of ray 


vec3 start = vColor.rgb; 


// calculate texture coords at fragment, 
// which is a fraction of window coords 


vec2 texc = gl FragCoord.xy/uWinDims. xy; 


// get end of ray by looking up back-face color 


vec3 end = texture(texBackFaces, texc).rgb; 


// calculate ray direction 


vec3 dir = end - start; 


// normalized ray direction 


vec3 norm dir = normalize(dir) ; 


// the length from front to back is calculated and 


// used to terminate the ray 


float len = length(dir.xyz); 


// ray step size 


float stepSize = 0.01; 


// X-ray projection 


vec4 dst = уес4(0.0); 


// step through the ray 


for(float t = 0.0; t < len; t += stepSize) { 


// set position to end point of ray 


vec3 samplePos = start + t*norm dir; 


// get texture value at position 
float val = texture(texVolume, samplePos).r; 


vec4 src = vec4(val); 


// set opacity 
src.a *= 0.1; 


src.rgb *- src.a; 


// blend with previous value 


dst = (1.0 - dst.a)*src + dst; 


// exit loop when alpha exceeds threshold 
if(dst.a >= 0.95) 


break; 


// set fragment color 
fragColor = dst; 
} 


class Camera: 
"""helper class for viewing""" 
def init (self): 
self.r - 1.5 
self.theta - 0 
self.center - [0.5, 0.5, 0.5] 
self.eye = [0.5 + self.r, 0.5, 0.5] 


self.up = [0.0, 0.0, 1.0] 


def rotate(self, clockWise): 
"""rotate eye by one step""" 
if clockWise: 
self.theta = (self.theta + 5) % 360 
else: 
self.theta = (self.theta - 5) % 360 


4 recalculate eye 


self.eye = [0.5 + self.r*math.cos(math. radians(self.theta) ) 
0.5 + self.r*math.sin(math.radians(self.theta)), 


0.5] 


class RayCastRender: 


"""Class that does Ray Casting""" 


def init (self, width, height, volume): 


"""RayCastRender constr""" 


4 create RayCube object 
self.raycube = raycube.RayCube(width, height) 
4 set dimensions 


self.width - width 


self.height height 


self.aspect = width/float(height) 


4 create shader 
self.program - glutils.loadShaders(strVS, strFS) 
# texture 


self.texVolume, self.Nx, self.Ny, self.Nz = volume 


4 initialize camera 


self.camera = Camera() 


def draw(self): 


# build projection matrix 


pMatrix = glutils.perspective(45.0, self.aspect, 0.1, 100.0 
) 


# modelview matrix 
mvMatrix = glutils.lookAt(self.camera.eye, self.camera.cent 
er, 
self.camera.up) 
# render 
# generate ray-cube back-face texture 
texture = self.raycube. renderBackFace(pMatrix, mvMatrix) 


# set shader program 


glUseProgram(self.program) 


4 set window dimensions 


glUniform2f(glGetUniformLocation(self.program, b"uWinDims") 


, 


float(self.width), float(self.height)) 


4 texture unit 0, which represents back-faces of cube 
glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 2D, texture) 


glUniformli(glGetUniformLocation(self.program, b"texBackFac 
es"), 0) 


# texture unit 1: 3D volume texture 
glActiveTexture(GL TEXTURE1) 


glBindTexture(GL TEXTURE 3D, self.texVolume) 


glUniformli(glGetUniformLocation(self.program, b"texVolume" 
), 1) 


draw front face of cubes 


self.raycube.renderFrontFace(pMatrix, mvMatrix, self.progra 


m) 


#self.render(pMatrix, mvMatrix) 


def keyPressed(self, key): 
if key == "1": 
self.camera.rotate(True) 
elif key == 'r': 


self.camera.rotate(False) 


def reshape(self, width, height): 
self.width - width 
self.height - height 
self.aspect - width/float(height) 


self.raycube.reshape(width, height) 


def close(self): 


self.raycube.close() 
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# set up vertex array object (МАО) 


self.vao = glGenVertexArrays(1) 


glBindVertexArray(self.vao) 


# define quad vertices 
vertexData = numpy.array([ 0.0, 1.0, 0.0, 
0.0, 0.0, 0.0, 
1.0, 1.0, 0.0, 
1.0, 0.0, 0.61, numpy.float32) 
# vertex buffer 
self.vertexBuffer = glGenBuffers(1) 


glBindBuffer(GL_ ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL ARRAY BUFFER, 4*len(vertexData), vertexData 


GL STATIC DRAW) 
4 enable arrays 
glEnableVertexAttribArray(self.vertIndex) 
# set buffers 


glBindBuf fer (GL ARRAY BUFFER, self.vertexBuffer) 


glVertexAttribPointer(self.vertIndex, 3, GL FLOAT, GL FALSE 
, 0, None) 


Я unbind МАО 


glBindVertexArray (0) 
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def draw(self): 
# clear buffers 
glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


# build projection matrix 


ЕБ = glutils.ortho(-0.6, 0.6, -θ.6, 0.6, 0.1, 100.0) 
# modelview matrix 
Ө mvMatrix = numpy.array([1.0, 0.0, 0.0, 0.0, 
0.0, 1.0, 0.0, 0.0, 
0.0, 0.0, 1.0, 0.0, 


-0.5, -0.5, -1.0, 1.0], numpy.flo 
at32) 


# use shader 
glUseProgram(self.program) 


# set projection matrix 


glUniformMatrix4fv(self.pMatrixUniform, 1, GL FALSE, pMatri 
x) 


4 set modelview matrix 


glUniformMatrix4fv(self.mvMatrixUniform, 1, GL FALSE, mvMat 
rix) 


4 set current slice fraction 
e 
glUniformlf(glGetUniformLocation(self.program, b"uSliceFrac 
"i 
float(self.currSlicelndex)/float(self.currSliceMax)) 
# set current slice mode 
ο 
glUniformli(glGetUniformLocation(self.program, b"uSliceMode 


г 


self.mode) 


# enable texture 
glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 3D, self.texture) 


glUniformli(glGetUniformLocation(self.program, b"tex"), 0) 


# bind VAO 


glBindVertexArray(self.vao) 


# draw 
glDrawArrays(GL TRIANGLE STRIP, 0, 4) 
Я unbind МАО 


glBindVertexArray(0) 
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# version 330 core 
in vec3 aVert; 


uniform mat4 uMVMatrix; 


uniform mat4 uPMatrix; 


uniform float uSliceFrac; 


uniform int uSliceMode; 


out vec3 texcoord; 


void main() { 


// x slice 
if (uSliceMode == 0) { 
e texcoord = vec3(uSliceFrac, aVert.x, 1.0-aVert.y); 
} 
// y slice 
else if (uSliceMode == 1) { 
ө texcoord = vec3(aVert.x, uSliceFrac, 1.0-aVert.y); 
} 
// z slice 
else { 


© texcoord = vec3(aVert.x, 1.0-aVert.y, uSliceFrac) ; 


// calculate transformed vertex 


gl Position = uPMatrix * uMVMatrix * vec4(aVert, 1.0); 
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# version 330 core 

Ө in vec3 texcoord; 

Ө uniform sampler3D texture; 
out vec4 fragColor; 


void main() { 


// look up color in texture 


© vec4 col = texture(tex, texcoord); 
о fragColor = col.rrra; 


} 
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def keyPressed(self, key): 
"""keypress handler""" 
if key == 'x': 
9 self.mode = SliceRender.XSLICE 
# reset slice index 
self.currSliceIndex = int(self.Nx/2) 
self.currSliceMax = self.Nx 
elif key == 'y': 
self.mode = SliceRender.YSLICE 
# reset slice index 
self.currSliceIndex = int(self.Ny/2) 


self.currSliceMax = self.Ny 


elif key == 'z': 
self.mode = SliceRender.ZSLICE 
# reset slice index 
self.currSliceIndex = int(self.Nz/2) 
self.currSliceMax = self.Nz 
elif key == 'L': 
self.currSliceIndex = (self.currSliceIndex + 1) % self.curr 
SliceMax 
elif key == 'r': 


self.currSliceIndex = (self.currSliceIndex - 1) 5 self.curr 
SliceMax 
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import OpenGL 


from OpenGL.GL import * 


from OpenGL.GL.shaders import * 


import numpy, math, sys 


import volreader, glutils 


strVS = 


# version 330 core 


in vec3 


uniform 


uniform 


uniform 


uniform 


aVert; 


mat4 uMVMatrix; 


mat4 uPMatrix; 


float uSliceFrac; 


int uSliceMode; 


out vec3 texcoord; 


void main() { 


// x slice 


if (uSliceMode == 


texcoord = vec3(uSliceFrac, aVert.x, 


) { 


1.0-aVert.y); 


} 
// у slice 
else if (uSliceMode == 1) { 


texcoord = vec3(aVert.x, uSliceFrac, 1.0-aVert.y); 


} 
// > slice 
else { 


texcoord = vec3(aVert.x, 1.0-aVert.y, uSliceFrac); 


// calculate transformed vertex 


gl Position = uPMatrix * uMVMatrix * vec4(aVert, 1.0); 


strFS = ипи 


4 version 330 core 


in vec3 texcoord; 


uniform sampler3D tex; 


out vec4 fragColor; 


void та1п() { 
// look up color in texture 
меса col = texture(tex, texcoord) ; 


fragColor = col.rrra; 


class SliceRender: 
# slice modes 


XSLICE, YSLICE, ZSLICE = 0, 1, 2 


def init (self, width, height, volume): 
"""SliceRender constructor""" 
self.width = width 


self.height 


height 


self.aspect - width/float(height) 


4 slice mode 


self.mode = SliceRender.ZSLICE 


4 create shader 


self.program = glutils.loadShaders(strVS, strFS) 


glUseProgram(self.program) 


self.pMatrixUniform = glGetUniformLocation(self.program, b 
'uPMatrix') 


self.mvMatrixUniform = glGetUniformLocation(self.program, 


b"uMVMatrix") 
# attributes 


self.vertIndex = glGetAttribLocation(self.program, b"aVert 
") 


4 set up vertex array object (VAO) 
self.vao = glGenVertexArrays(1) 


glBindVertexArray(self.vao) 


# define quad vertices 
vertexData = numpy.array([ 0.0, 1.0, 0.0, 
0.0, 0.0, 0.0, 


1.0, 1.0, 0.0, 


1.0, 0.0, 0.0], numpy. float32) 
# vertex buffer 


self.vertexBuffer = glGenBuffers(1) 


glBindBuffer(GL_ ARRAY BUFFER, self.vertexBuffer) 


glBufferData(GL_ ARRAY BUFFER, 4*len(vertexData), vertexDat 
a, 
GL STATIC DRAW) 
4 enable arrays 


glEnableVertexAttribArray(self.vertIndex) 


# set buffers 


giBindBuffer(GL ARRAY BUFFER, self.vertexBuffer) 


glVertexAttribPointer(self.vertIndex, 3, GL FLOAT, GL FALS 
E, 0, None) 


Я unbind МАО 


glBindVertexArray (0) 


# load texture 


self.texture, self.Nx, self.Ny, self.Nz = volume 


# current slice index 


self.currSliceIndex = int(self.Nz/2); 


self.currSliceMax = self.Nz; 


def reshape(self, width, height): 
self.width = width 


self.height 


height 


self.aspect = width/float (height) 


def draw(self): 


# clear buffers 


glClear(GL COLOR BUFFER BIT | GL DEPTH BUFFER BIT) 


4 build projection matrix 


pMatrix - glutils.ortho(-0.6, 0.6, -0.6, 0.6, 0.1, 100.0) 


4 modelview matrix 


mvMatrix - numpy.array([1.0, 0.0, 0.0, 0.0, 


0.0, 1.0, 0.0, 0.0, 


0.0, 0.0, 1.0, 0.0, 


-0.5, -0.5, 
at32) 


# use shader 


glUseProgram(self.program) 


-1.0, 1.0], numpy.flo 


# set projection matrix 


glUniformMatrix4fv(self.pMatrixUniform, 1, GL FALSE, pMatri 
x) 


# set modelview matrix 


glUniformMatrix4fv(self.mvMatrixUniform, 1, GL FALSE, mvMat 
rix) 


# set current slice fraction 


glUniformlf(glGetUniformLocation(self.program, b"uSliceFrac 
mj 


float(self.currSlicelndex)/float(self.currSliceMax)) 


# set current slice mode 


glUniformli(glGetUniformLocation(self.program, b"uSliceMode 
pr 


self.mode) 


4 enable texture 
glActiveTexture(GL TEXTUREO) 


glBindTexture(GL TEXTURE 3D, self.texture) 


glUniformli(glGetUniformLocation(self.program, b"tex"), 0) 


# bind VAO 
glBindVertexArray(self.vao) 

# draw 

glDrawArrays(GL TRIANGLE STRIP, 0, 4) 
à unbind МАО 


glBindVertexArray (0) 


def keyPressed(self, key): 

"""keypress handler""" 

if key == 'x': 
self.mode = SliceRender.XSLICE 
# reset slice index 
self.currSliceIndex = int(self.Nx/2) 
self.currSliceMax = self.Nx 

elif key == 'y': 
self.mode = SliceRender.YSLICE 
# reset slice index 
self.currSliceIndex = int(self.Ny/2) 
self.currSliceMax = self.Ny 

elif key == 'z': 
self.mode = SliceRender.ZSLICE 


# reset slice index 


self.currSliceIndex = int(self.Nz/2) 
self.currSliceMax = self.Nz 
elif key == 'l': 
self.currSliceIndex = (self.currSliceIndex + 1) % self.curr 
SliceMax 
elif key == 'r': 


self.currSliceIndex = (self.currSliceIndex - 1) 5 self.curr 
SliceMax 


def close(self): 


pass 
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# load volume data 

e self.volume = volreader.loadVolume (imageDir) 
# create renderer 

e 


self.renderer - RayCastRender(self.width, self.height, self 
. Volume) 
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def onKeyboard(self, win, key, scancode, action, mods): 


Я print ‘keyboard: ', win, key, scancode, action, mods 
# ESC to quit 
if key is glfw.GLFW KEY ESCAPE: 
self.renderer.close() 
self.exitNow = True 
else: 


о 
if action is glfw.GLFW PRESS or action is glfw.GLFW_REPEAT: 


if key == glfw.GLFW KEY V: 
# toggle render mode 
6 
if isinstance(self.renderer, RayCastRender) : 
self.renderer = SliceRender(self.width, self.height, 
self.volume) 


else: 


self.renderer = RayCastRender(self.width, self.height, 


self.volume) 


4 call reshape on renderer 


self.renderer.reshape(self.width, self.height) 
else: 


# send keypress to renderer 
© 
keyDict = {glfw.GLFW KEY X: 'x', glfw.GLFW KEY Y: 'y', 


glfw.GLFW KEY Z: 'z', glfw.GLFW KEY LEFT: 'l', 
glfw.GLFW KEY RIGHT: 'г') 
try: 
self.renderer.keyPressed(keyDict[key]) 
except: 
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er/volrender/[][]volrender.py[][][] 
import sys, argparse, os 
from slicerender import * 
from raycast import * 


import glfw 
class RenderWin: 
"""GLFW Rendering window class""" 


def init (self, imageDir): 


# save current working directory 


сыа = os.getcwd() 


# initialize glfw; this changes cwd 


glfw.glfwInit() 


# restore cwd 


os.chdir(cwd) 


# version hints 
gl fw. gl fwWindowHint (gl fw.GLFW CONTEXT VERSION MAJOR, 3) 
gl fw. gl fwWindowHint(glfw.GLFW CONTEXT VERSION MINOR, 3) 
манны рО GL TRU 


gl fw.gl fwWindowHint (glfw.GLFW OPENGL PROFILE, 


glfw.GLFW OPENGL CORE PROFILE) 


# make a window 

self.width, self.height = 512, 512 

self.aspect = self.width/float(self.height) 
self.win = glfw.glfwCreateWindow(self.width, self.height, b 
"volrender") 

# make context current 


glfw.glfwMakeContextCurrent(self.win) 


# initialize GL 
glViewport(0, 0, self.width, self.height) 
glEnable(GL DEPTH TEST) 


glClearColor(0.0, 0.0, 0.0, 0.0) 


# set window callbacks 
glfw.glfwSetMouseButtonCallback(self.win, self.onMouseButto 
n) 

glfw.glfwSetKeyCallback(self.win, self.onKeyboard) 


glfw.glfwSetWindowSizeCallback(self.win, self.onSize) 


# load volume data 
self.volume = volreader.loadVolume(imageDir) 


4 create renderer 


self.renderer - RayCastRender(self.width, self.height, self 
. Volume) 


4 exit flag 


self.exitNow = False 


def onMouseButton(self, win, button, action, mods): 


4 print 'mouse button: ', win, button, action, mods 


pass 


def onKeyboard(self, win, key, scancode, action, mods): 


4 print 'keyboard: ', win, key, scancode, action, mods 


# ESC to quit 

if key is glfw.GLFW KEY ESCAPE: 
self.renderer.close() 
self.exitNow = True 


else: 


if action is glfw.GLFW PRESS or action is glfw.GLFW REPEAT: 
if key == glfw.GLFW KEY V: 


# toggle render mode 


if isinstance(self.renderer, RayCastRender): 


self.renderer = SliceRender(self.width, self.height, 
self.volume) 


else: 


self.renderer - RayCastRender(self.width, self.height, 
self.volume) 


8 call reshape on renderer 


self.renderer.reshape(self.width, self.height) 
else: 


# send keypress to renderer 


keyDict = {glfw.GLFW KEY X: 'x', glfw.GLFW KEY Y: 'y', 


glfw.GLFW KEY Z: 'z', glfw.GLFW KEY LEFT: 'l', 
glfw.GLFW KEY RIGHT: 'г'} 
try: 
self.renderer.keyPressed(keyDict [key] ) 
except: 


pass 


def onSize(self, win, width, height): 
#print 'onsize: ', win, width, height 


self.width - width 


self.height height 
self.aspect - width/float(height) 
glViewport(0, 0, self.width, self.height) 


self.renderer.reshape(width, height) 


def run(self): 
# start Loop 
while not glfw.glfwWindowShouldClose(self.win) and not self 
.exitNow: 
# render 
self.renderer.draw() 


# swap buffers 


glfw.glfwSwapBuffers(self.win) 
# wait for events 
glfw.glfwwaitEvents() 

# end 


glfw.glfwTerminate() 


4 main() function 
def main(): 
print('starting volrender...') 
# create parser 
parser = argparse.ArgumentParser(description="Volume Render 
ing...") 
# add expected arguments 


parser.add argument ( ' - - 
dir', dest='imageDir', required=True) 


# parse args 


args = parser.parse args() 
# create render window 
rwin = RenderWin(args.imageDir) 


rwin.run() 


# call main 


if _ name == ' main ': 


main() 
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$ python volrender.py --dir mrbrain-8bit/ 
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#include "Arduino.h" 


гг LED pin numbers (digital) 
int pinRed = 4; 
int pinGreen = 2; 
void setup() 
{ 
ff initialize serial comms 
Serial .begin(9600); 


// set pins 
pinMode(pinRed, OUTPUT); 
pinMode(pinGreen, OUTPUT); 


`x 


void loop() 
£ 
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#include "Arduino.h" 


void setup() 


// initialize serial communications 
Ө Serial.begin(9600) ; 


} 


void loop() 

{ 
// read АӨ 

Ө int vall = analogRead(Q) ; 
// read Al 

Ө int val2 = analogRead(1); 
// print to serial 

Ө Serial.print(vall); 
Serial.print(" "); 
Serial.print(val2); 
Serial.print("\n"); 
// wait 

Ө delay(50); 


} 
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12.5 Python] 
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class AnalogPlot: 
# constructor 
def | init (self, strPort, maxLen): 
4 open serial port 


e self.ser = serial.Serial(strPort, 9600) 


о self.a0Vals = deque( [0.0] *maxLen) 
© self.alVals = deque( [0.0] Ἴπαχι еп) 


ο self.maxLen = maxLen 
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# add data 
def add(self, data): 
assert(len(data) == 2) 


self.addToDeq(self.aOVals, data[0]) 


° 


self.addToDeq(self.alVals, data[1]) 


© 


# add to deque; pop oldest value 
def addToDeq(self, buf, val): 
© buf .pop() 


ο buf .appendleft(val) 
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4 update plot 


def update(self, frameNum, a0, al): 


try: 
e line = self.ser.readline() 
o data = [float(val) for val in line.split()] 


4 print data 


if(len(data) == 2): 
e self.add(data) 
о a0.set data(range(self.maxLen), self.a0Vals) 
e al.set data(range(self.maxLen), self.alVals) 
except: 
о ра55 


return a0, al 
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4 clean up 

def close(self): 
# close serial 
self.ser.flush() 


self.ser.close() 


Umain(QQU000000matplotlib |[ II 


# set up animation 

fig = plt.figure() 

ax = plt.axes(xlim=(0, maxLen), ylim=(0, 1023)) 
a0, = ax.plot([], []) 


al, = ax.plot([], []) 


о о о о © 


anim = animation.FuncAnimation(fig, analogPlot.update, 


fargs-(a0, al), interval-20) 


# show plot 


Ө plt.show() 
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# create parser 


parser = argparse.ArgumentParser(description="LDR serial") 
# add expected arguments 


parser.add argument ( ' - - 
port', dest='port', required=True) 


parser.add argument ( ' - - 
N', dest='maxLen', required=False) 


# parse args 


args = parser.parse args() 


strPort = args.port 


4 plot parameters 
maxLen = 100 
if args.maxLen: 


maxLen = int(args.maxLen) 


--рогібОобрбоДОбОДОД0О0000000000Агачіпо 
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ΠΠΠΠΠΠΠΠΠΗΡΥΓΠΟΠΠΠΠΠΠΠΗ 
https://github.com/electronut/pp/tree/mast 


er/arduino-ldr/ldr.pyQQ000000 
import serial, argparse 
from collections import deque 
import matplotlib.pyplot as plt 


import matplotlib.animation as animation 


# plot class 
class AnalogPlot: 
# constructor 
def | init (self, strPort, maxLen): 


4 open serial port 


self.ser = serial.Serial(strPort, 9600) 


self.a0Vals deque( [0.0] *maxLen) 


self.alVals = deque([0.0] *maxLen) 


self.maxLen = maxLen 


# add data 

def add(self, data): 
assert(len(data) -- 2) 
self.addToDeq(self.aOVals, data[0]) 


self.addToDeq(self.alVals, data[1]) 


# add to deque; pop oldest value 
def addToDeq(self, buf, val): 
buf.pop() 


buf.appendleft(val) 


# update plot 
def update(self, frameNum, a0, al): 
try: 
line = self.ser.readline() 
data = [float(val) for val in line.split()] 


# print data 


if(len(data) == 2): 
self.add(data) 
a0.set data(range(self.maxLen), self.a0Vals) 


al.set data(range(self.maxLen), self.alVals) 


except: 


pass 


return a0, al 


4 clean up 

def close(self): 
# close serial 
self.ser.flush() 


self.ser.close() 


Я main() function 
def main(): 


# create parser 


parser = argparse.ArgumentParser(description="LDR serial") 
# add expected arguments 


parser.add argument (' - - 
port', dest='port', required=True) 


parser.add argument('-- 
N', dest='maxLen', required=False) 


# parse args 


args = parser.parse args() 


#strPort = '/dev/tty.usbserial-A7006Ygh' 


strPort = args.port 


print('reading from serial port %s...' % strPort) 


# plot parameters 
maxLen = 100 
if args.maxLen: 


maxLen = int(args.maxLen) 


# create plot object 


analogPlot = AnalogPlot(strPort, maxLen) 


print('plotting data...') 


# set up animation 


fig = plt.figure() 


ax = plt.axes(xlim=(0, maxLen), ylim=(0, 1023) ) 


ад, = ax.plot([], [1) 
al, = ax.plot([], П) 
anim = animation.FuncAnimation(fig, analogPlot.update, 


fargs-(a0, al), interval=20) 


# show plot 


plt.show() 


4 clean up 


analogPlot.close() 


print('exiting.') 


# call main 
if name == ' main ': 


main() 
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$ python3 --port /dev/tty.usbmodem411 ldr.py 
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// motor A connected to A01 and A02 


// motor B connected to В01 and B02 
€ int STBY - 10; //standby 


// Motor A 
int PWMA - 3; //speed control 
int AIN1 = 9; //direction 


int AIN2 = 8; //direction 
// Motor B 
int PWMB = 5; //speed control 
int BIN1 = 11; //direction 

Ө int BIN2 = 12; //direction 
void setup(){ 


©  pinMode(STBY, OUTPUT); 


pinMode(PWMA, OUTPUT); 


pinMode(AIN1, OUTPUT); 


pinMode(AIN2, OUTPUT); 


pinMode(PWMB, OUTPUT) ; 
pinMode(BIN1, OUTPUT); 


pinMode(BIN2, OUTPUT); 


// initialize serial communication 
O  Serial.begin(9600); 


) 
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// main loop that reads the motor data sent by laser.py 
void loop() 
{ 


// data sent is of the form 'H' (header), speedl, dirl, 
speed2, dir2 


Ө if (Serial.available() >= 5) { 

e if(Serial.read() == 'H') { 
// read the next 4 bytes 

e byte s1 - Serial.read(); 
byte 41 = Serial.read(); 
byte s2 - Serial.read(); 


byte 42 = Serial.read(); 


// stop the motor if both speeds are 0 
ο if(sl == 0 && 52 == 0) í 
stop(); 
} 
else { 
// set the motors' speed and direction 
e move(0, 51, 41); 
move(1, s2, d2); 
} 


// slight pause for 20 ms 


о детау( (20); 
) 
else { 
// if there is invalid data, stop the motors 
© stop(); 
} 
} 
else { 
// if there is no data, pause for 250 ms 


Ө delay(250) ; 


} 
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// set motor speed and direction 
// motor: A -> 1, В -> 0 
// direction: 1/0 
void move(int motor, int speed, int direction) 
{ 
// disable standby 


Ө digitalWrite(STBY, HIGH) ; 


LOW; 


6 boolean inPinl 


boolean inPin2 = HIGH; 


Ө if(direction == 1){ 
inPinl = HIGH; 


inPin2 = LOW; 


if(motor == 1){ 
ο digitalWrite(AIN1, inPinl); 
digitalWrite(AIN2, inPin2) ; 


analogWrite(PWMA, speed); 


else{ 


e digitalWrite(BIN1, inPinl); 
digitalWrite(BIN2, inPin2); 


analogWrite(PWMB, speed); 


} 
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digitalWrite ПДОДОДОДОДО0НІЄНОЇО000000 
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О00000000000РММрП0000апаіод\№гіќе()00 
00010, 2551П0000Агачіпо ОО0000000 
digitalWrite()QU00001 0000000000 
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void stop(){ 
//enable standby 


digitalWrite(STBY, LOW); 


} 


13.4 РуїһопПП 


ООДОДД000000ОРуспопООДОДООД00000000 
ОООО000РРТОООООООООАмаці по 00001 3.5000 
О00000000 
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р = pyaudio.PyAudio() 


ΠΠΠΠΠΠΗΡνϑ8α1ο0ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
де при  емсе()0000000000 


# get pyaudio input device 
def getInputDevice(p): 
Ө index = None 


Ө nDevices = p.get device count() 


print('Found %d devices. Select input device:' 5 nDevices) 
4 print all devices found 
for i in range(nDevices): 

e deviceInfo - p.get device info by index(i) 


ο devName = deviceInfo['name'] 


e print("%d: 555" % (i, devName)) 
# get user selection 
try: 
# convert to integer 
о index = int(input() ) 
except: 


pass 


# print the name of the chosen device 

if index is not None: 
devName = p.get device info by index(index) ["name" ] 
print("Input device chosen: %s" % devName) 


ο return index 


ПӨПППїпаехПППППМопе[П[ПїпдехП ӨПППППП 
О0000000МопердродО00О0О00000000Фп000 

де аеуісе соџпё()ПО000000000000000000 
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e[]|[get device info Бу іпаех()0000000000 
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ердробріари 00000000000000000000000000 
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# set FFT sample length 
Ө fftLen = 2**11 
# set sample rate 


Ө sampleRate = 44100 


print('opening stream...') 
Ө stream = p.open(format = pyaudio.palntl6, 
channels = 1, 
rate = sampleRate, 
input = True, 
frames per buffer = fftLen, 


input_device index = inputIndex) 


ΠΘΠΠΠΓΓΤΠΠΠΠΠΠΠΠΠΠΠΕΓΤΠΠΠΠΠΠΗΠΠΗ 
2048Π1211ΠΓΓΤΠΠΠΠΖΗΠΠΠΠΗΠΗΠΠΠΠΠΠΗ 
руаийаїоППППППП44100ПП44.1 kHz ΘΠΠΠΟΡ 
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. руацйдїо.раїпі16000000000260000 
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. ratellIUUUUUUUU44100 Hz[] 

. input[][]|]]True[] 

« frںames_per_bufferlUUFFTUUUUUDU‎ 

- input_device_indexQQU000 
getinputDevice()Q00000000 
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4 read a chunk of data 
e data = stream.read(fftLen) 


# convert the data to a numpy array 


e 
dataArray = numpy.frombuffer(data, dtype=numpy.int16) 


Пепоооооовообобе епррооопепооооооо 
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ООО000000РРТП 


# get FFT of data 
e fftVals = numpy.fft.rfft(dataArray)*2.0/fftLen 
# get absolute values of complex numbers 


6 fftVals = numpy.abs(fftVals) 
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# average 3 frequency bands: 0-100 Hz, 100- 
1000 Hz, and 1000-2500 Hz 


levels = [numpy.sum(fftVals[0:100])/100, 
numpy.sum(fftVals[100:1000])/900, 


numpy.sum(fftVals[1000:2500])/1500] 


О000000000000000003000000100 Н21100Г 
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4 'H' (header), speedl, dirl, speed2, dir2 


e vals = [ord('H'), 100, 1, 100, 1] 
# speed1 
e vals[1] = int(5*levels[0]) % 255 


# speed2 


e vals[3] = int(100 + levels[1]) % 255 
# dir 
d1 = 0 
о if levels[2] > 0.1: 
d1 = 1 
vals[2] = dl 


e vals[4] = 0 
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4 automatic test for sending motor speeds 
def autoTest(ser): 
print('starting automatic test...') 
try: 
while True: 
# for each direction combination 
о Тог dr in [(0, 0), (1, 0), (0, 1), (1, 1)]: 
# for а range of speeds 
for j in range(25, 180, 10): 
for i in range(25, 180, 10): 
vals = [ord('H'), i, dr[0], j, аг[1]] 
print(vals[1:]) 


data = struct.pack('BBBBB', *vals) 


o о o о © © 


ser.write(data) 
sleep(0.1) 
except KeyboardInterrupt: 
print('exiting...') 
# shut off motors 


Ө vals = [ord('H'), 0, 1, 0, 1] 


data = struct.pack('BBBBB', *vals) 
ser.write(data) 


ser.close() 
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13.4.7 ΠΠΠΠΠ 
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Я main method 
def main(): 


# parse arguments 


parser = argparse.ArgumentParser(description='Analyzes audi 
o input and 


sends motor control information via serial port') 
# add arguments 


parser.add argument ( ' - - 
port', dest-'serial port name', required=True) 


parser.add argument('-- 
mtest', action='store true', default-False) 


parser.add argument('-- 
atest', action='store true', default=False) 


args - parser.parse args() 


ООО000000000000000000000000000000000 
ОООО0000000000000000000000000 


О00таіп()00000000000000000 


4 open serial port 
strPort = args.serial port name 
print('opening ', strPort) 
Ө ser = serial.Serial(strPort, 9600) 
if args.mtest: 
manualTest(ser) 
elif args.atest: 
autoTest (ser) 
else: 


ө fftLive(ser) 


Le ri peyseriralt in i panim amm m n Oan p 
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4 manual test of motor direction and speeds 
def manualTest(ser): 

print('starting manual test...') 

try: 


while True: 


print('enter motor control info such as < 100 1 120 0 >') 


e strIn = raw input() 
e vals = [int(val) for val in strIn.split()[:4]] 
e vals.insert(0, ord('H')) 
ο data = struct.pack('BBBBB', *vals) 
Θ ser.write(data) 
except: 


print('exiting...') 
# shut off the motors 


9 vals = [ord('H'), 0, 1, 0, 1] 


data = struct.pack('BBBBB', *vals) 
ser.write(data) 


ser.close() 
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ΠΠΠΠΠΠΠΠΠΗΡΥΓΠΟΠΠΠΠΠΠΠΗ 
https://github.com/electronut/pp/tree/mast 


er/arduino-laser/laser.py[]L]LL] 
import sys, serial, struct 
import pyaudio 
import numpy 
import math 
from time import sleep 


import argparse 


4 manual test of motor direction speeds 


def manualTest(ser): 


print('staring manual test...') 
try: 


while True: 


print('enter motor control info: eg. < 100 1 120 0 >') 
strIn = raw input() 
vals = [int(val) for val in strIn.split()[:41] 
vals.insert(0, ord('H')) 
data = struct.pack('BBBBB', *vals) 
ser.write(data) 
except: 
print('exiting...') 


# shut off motors 


vals [ord('H'), 0, 1, 0, 1] 
data = struct.pack('BBBBB', *vals) 
ser.write(data) 
ser.close() 
# automatic test for sending motor speeds 
def autoTest(ser): 
print('staring automatic test...') 
try: 
while True: 


# for each direction combination 


for dr in [(0, 0), (1, 0), (0, 1), (1, 1)1: 
# for a range of speeds 
for j in range(25, 180, 10): 
for i in range(25, 180, 10): 
vals = [ord('H'), i, dr[0], j, аг[1]] 
print(vals[1:]) 
data = struct.pack('BBBBB', *vals) 
ser.write(data) 
sleep(0.1) 
except KeyboardInterrupt: 
print('exiting...') 
# shut off motors 
vals = [ord('H'), 0, 1, 0, 1] 
data = struct.pack('BBBBB', *vals) 
ser.write(data) 


ser.close() 


# get pyaudio input device 
def getInputDevice(p): 
index - None 


nDevices - p.get device count() 


print('Found %d devices. Select input device:' % nDevices) 


4 print all devices found 
for i in range(nDevices): 
deviceInfo = p.get device info by index(i) 
devName = deviceInfo['name'] 
print("%d: 555" % (1, devName) ) 
# get user selection 
try: 
# convert to integer 
index = int(input() ) 
except: 


pass 


# print the name of the chosen device 

if index is not None: 
devName = p.get device info by index(index) ["name" ] 
print("Input device chosen: %s" % devName) 


return index 


Я FFT of live audio 
def fftLive(ser): 
# initialize pyaudio 


р = pyaudio.PyAudio() 


# get pyAudio input device index 


inputIndex = getInputDevice(p) 


# set FFT sample length 
fftLen = 2**11 
# set sample rate 


sampleRate = 44100 


print('opening stream...') 
stream = p.open(format = pyaudio.paInt16, 
channels = 1, 
rate - sampleRate, 
input - True, 
frames per buffer - fftLen, 
input device index - inputIndex) 
try: 
while True: 
# read a chunk of data 
data = stream.read(fftLen) 


# convert to numpy array 


dataArray = numpy.frombuffer(data, dtype=numpy.int16) 


# get FFT of data 

fftVals = numpy.fft.rfft(dataArray)*2.0/fftLen 
# get absolute values of complex numbers 
fftVals = numpy.abs(fftVals) 


# average 3 frequency bands: 0-100 Hz, 100- 
1000 Hz and 1000-2500 Hz 


levels = [numpy.sum(fftVals[0:100])/100, 
numpy.sum(fftVals[100:1000])/900, 


numpy.sum(fftVals[1000:2500])/1500] 


£ the data sent is of the form: 
4 'H' (header), speedl, dirl, speed2, dir2 


vals = [ord('H'), 100, 1, 100, 1] 


# speedl 
vals[1] = int(5*levels[0]) % 255 
Я speed2 


vals[3] = int(100 + levels[1]) % 255 


if levels[2] > 0.1: 


dl = 1 


vals[2] = dl 
vals[4] = 0 
# pack data 


data = struct.pack('BBBBB', *vals) 
# write data to serial port 
ser.write(data) 
# a slight pause 
sleep(0.001) 
except KeyboardInterrupt: 
print('stopping...') 
finally: 
print('cleaning up') 
stream.close() 
p.terminate() 
# shut off motors 
vals = [ord('H'), 0, 1, 0, 1] 
data = struct.pack('BBBBB', *vals) 
ser.write(data) 
4 close serial 
ser.flush() 


ser.close() 


# main method 
def main(): 
# parse arguments 
parser = argparse.ArgumentParser(description='Analyzes audi 
o input and 
sends motor control information via serial port') 
# add arguments 


parser.add argument ( ' - - 
port', dest='serial_port_name', required=True) 


parser.add argument ( ' - - 
mtest', action='store true', default=False) 


parser.add argument ( ' - - 
atest', action='store true', default=False) 


args = parser.parse args() 


# open serial port 
strPort = args.serial port name 
print('opening ', strPort) 
ser - serial.Serial(strPort, 9600) 
if args.mtest: 

manualTest(ser) 
elif args.atest: 

autoTest(ser) 


else: 


fftLive(ser) 


4 call main function 
if name == ' main ': 


main() 
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5 python3 laser.py --port /dev/tty.usbmodem411 --atest 
('opening ', '/dev/tty.usbmodem1411') 

staring automatic test... 

[25, 0, 25, 0] 

[35, 0, 25, 0] 


[45, 0, 25, 9] 
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$ python3 laser.py --port /dev/tty.usbmodem411 
('opening ', '/dev/tty.usbmodem1411') 

Found 4 devices. Select input device: 

0: Built-in Microph 

1: Built-in Output 

2: BoomDevice 

3: AirParrot 

0 

Input device chosen: Built-in Microph 


opening stream... 
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rms = numpy.sqrt(numpy.mean(x**2)) 
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[1] “Relays and Optoisolators,” What- 
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0000600 


Pi-compatible power Pi case 


supply 


Pi-compatible USB 
wwe Wi-Fi adapter 


ADAPTER 


8GB or higher-capacity SD card 


ЛЕ 


000000000000 
GOD000000008 
000000000000 

{2000000006 ге 
500000000000 Composite video cable 
Qoseocaoonoo 





[14-3 ΠΠΠΠΠΠΠΠΠΠ 


О00000000000РІО0000000000001№р://е!іпихогд/ 
ЕРІ VerifiedPeripherals[] 


14.2 ПОООООП 


ООООО0О00000000Руєпопо000000000000000 
OUUUUUUUUUURaspberry Pi Foundation 
(“Getting Started with NOOBS” 00000000 
Π http://www.raspberrypi.org/help/noobs- 
setup/[][] 


14.2.1 ЦОГО 


ООДОООД00000000005вО000000000000000 
00000000ВазрібіапродроробО00000000000 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗΠπωΧ мікір 
П“АРІ Easy SD Card Setup” 000 
http://elinux.org/RPi Easy SD Card Setup[] 


ΠΠΠΙΠΠΠΠΠ“υ5ίπο МООВ5”ППППППППП 
14.2.2 ΠΠΠΗ 


ОДО0О000000000000000000520000000000 
ООООО00000000000000000000000000000000 
Д0Огавбрі-соплд ДООДОД00000000000000000 
Raspberry Pi Foundation] ]]]]]raspi-config|] 


ПП 





https://www.raspberrypi.org/documentatio 
n/ configuration/raspi-config.md 00000000 


1000ЕЄхрапа Ееѕуѕќет00005000 


2ПППЕпабіе Boot to Desktop/Scratch([]LHLL 
Desktop[] 


3000Change Time ΖΟΠΕΠΠΗ 
Internationalization Options II] 


A[]]ljAdvanced Options[][]|J Overscan[] 


500055НОПОООЕпавбте or Disable SSH 
Ѕегуег0ПАдамапсеа Оріїоп5ПД000000000 
LIE] 


ΠΠΠΠΓΙΠΙ5ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
14.2.3 ΜΜΙ-ΕΙΠΠ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΜΤΊ-ΕΠΠΠΠΗ 
ΠΠΠΠΙρ://6ε!1πΧ.ΟΓΩ/ΠΠΠΠΒΘ5ΡΟΙΘΠΠΠΠΠΠΠΠ 
ΠΠΠΠΠΗΠΠΠΠΠΠΗΠΗΠΗΠΠΠΠΠΠΠΠΗΙΡΉΠΠΗΠΗΠΠΗ 

Ν8ΠΟΠΠΠΠΠΗΠΠΠΠΠΠΠΗΝ8ΠΟΠΠΠΗΠΠΠΗΙΠΠΠΠΗ 
| μ.μ” 





ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΠΙ Χτειπιἰπα!ΠΠΠΠηΒαΣροίαη 
ООО00000000000 


$ sudo nano /etc/network/interfaces 


ООД000і пег'асев ДООООООДООО000000000 


auto lo 


iface lo inet loopback 


iface ethO inet dhcp 


allow-hotplug wlanO 

iface wlanO inet manual 

wpa-roam /etc/wpa supplicant/wpa supplicant.conf 
iface default inet static 

address 192. x.x.x _ 

netmask 255.255.255.0 


gateway 192. х.х.х _ 


ППППППППааагеввПпевтавКПдағемауППП 
О00000000000000000000255.255.255.0000 
ОДО00000ї їпихоДііпихДО0ієсоп 900 
Windows] Windows ДаДОДОДОЇ осопйа/ а ПО 
0005 ХППО5узіеги Preferences- 


> Мермогк ПОДООДООД000ТРООООДООООД00000 
ОДОІРОПО 


ΠΙΠΠΠΝΝΙΓΙ СопладоророрООМЛ-РІДОДОДОДО 


ООДОО0000000000000000 
https://learn.adafruit.com/[jJAdafruit 000 


ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΗΠΠΠΠΜΙἀΟΓΙΠΠΠΠΠΠΗ 


14.2.4 [I 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΒΡΊ.ΘΡΙΟΠΗ 
Bottle ΝΥΕΌΠΠΠΠΠΒΘ85ρΌΘΓΓν РІДОДОПРУ ОПП 


ПООООООООООоооооооо00000000000000 


$ sudo apt-get update 


ча 


sudo apt-get install python-setuptools 


fr 


sudo apt-get install python-dev 


ux 


sudo apt-get install python-rpi.gpio 


ча 


sudo easy_install bottle 


ПППП http://www.flotcharts.org/ ППППППППОвС 
JavaScriptQ 0000000000 о bp gp I ID 
ЕНЕЕЕЕЕЕЕ 


$ wget http://www. flotcharts.org/downloads/flot-x.zip 
$ unzip flot-x.zip 


$ mv flot myProjectDir/ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΑαα/Γμ Ργίποη ΡΗΤΗ 
[]https://github.com/ 
adafruit/Adafruit Python ὈΗΤ/ΠΠΠΠΠΠΠΠΠΙΠΠ 


ШШШОНТ11ШШШИШШШ 
$ 
git clone https://github.com/adafruit/Adafruit Python DHT.g 
it 


$ cd Adafruit Python DHT 


$ sudo python setup.py install 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 
14.2.5 ΠΠ55ΗΠΠ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΝ 
ΠΠΠΠΠΠΗΠΏΧΠΟ5 ХО00000005есиге Shell 
ы 


О00000000055Нр00 


€ moksha:- mahesh$ ssh pi@192.168.4.32 
Ө pi@192.168.4.32's password: 
Θ pi@raspberrypi - $ whoami 

pi 


pi@raspberrypi - $ 


ΠΠΠΠΠΠΠΠΘΙΙΠΙΠΠΞΞΠΙΠΠΠΠΙΠΠΠΠΠΡΙΠΠΠΙΡΠΠΠΠΗ 
ssh изеглате(атр_адагез ПИПООООООООООО 
ΠΠ55ΠΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΗΓΒ5ΡΌΘΠΥΠ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠνν πο8πη!ΠΠΘΠΠΠΠΠΠΡΙΠΠΠ 
ЕНЕЕЕЕЕЕЕШЕ 


НО ЫШЫ ЫШЫ 
ПСП 


14.2.6 Мер Во е 


ΠΠΠν/ΕΌΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΝΥΕΌΠΠΠΠΠΠΠ 
ΠΠΒΟΙΙΕΠΠΠΠΠΠΠΠΠΠΠΠΡγἴΠοη ΝΨΕΌΠΠΠΠΠΠΠ 
0000006оёе.руро000000000000Воє е0 
ООО00000000 


from bottle import route, run 
@route('/hello' ) 
def hello(): 


return "Hello Bottle World!" 


run(host='192.168.x.x', port=xxxx, debug=True) 


ΠΠΠΠΠΡΥϊποπΠΠΘτουτε[ΙΠΠΠΠΠΠΠΠΠΠΠυΡΙΗ 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ 000000 
run()000B ottleQQ000000000000000000000 
ПППППІРИППИППИППИГПППГППаеочо ІІ гиеП 
ΕΕ ΕΕ ΕΕ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕ 
http://192.168.4.4: 8080/hello[]ULIE]LILILIL] 
Bottle[]| LLULLU Hello Bottle Мога” 0o 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΝΥΘΕΌΠΠΠΗ 


ΠΠΠΠΠΠΠΠΙανα5εΠρΙΙΧΜΙ ΠΑΙΑΧΠΠΠΠΠΠΠΠΠΗ 
ΠΠΠΠΗΠΠΠΒοϊΠΕΠΠΠΠΠΠΠΠΠΑΙΑΧΠΠΠΠΠΠΠΠΠΠΠ 
ПШПШОчегуЦП 


[ΟΠΠ 


Python( 000000 @000000000 000000000000 III! 
О0000000000*00" 0000000000000 


@wrapper 


ОДОО000000 


ОДОРуєпопОДО000000000000 
14.2.7 ΠΠΟΓΠΠ 


ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΟΙΤΠΠΠΠΠΠΠΠΠΠΑΡΙΠΠΠΗ 
ОООДОО000000000000000000000000000 
ОНТМЕДОДОООООООО0ОД00000000РІС0Б00000 
ООООО000000000000005іпаріе-йоє. hem TE] 1L] 
ΕΠΗ 


<html> 
<head> 


«meta http-equiv="Content- 
Type" content="text/html; charset=utf-8"> 


<title>SimpleFlot</title> 
Ө <style> 
.demo-placeholder { 
width: 80%; 
height: 80%; 
} 
</style> 
Ө <script language="javascript" type="text/javascript" 
src="flot/jquery.js"></script> 
«script language-"javascript" type="text/javascript" 
src="flot/jquery.flot.js"></script> 
«script language-"javascript" type="text/javascript"> 
€ $(document).ready(function() ή 
// create plot 


Ө var data = ЇЇ; 


for(var i = 0; i < 500; i ++) { 


data.push([i, Math.exp(- 
1/100) *Math.sin(Math.PI*i/10) ] ); 


} 
Ө var plot = $.plot("#placeholder", [data]); 
}); 
</script> 


</head> 


<body> 
<h3>A Simple Flot Plot</h3> 
<div class="demo-container"> 
@ <div id="placeholder" class="demo-placeholder"></div> 
</div> 
</body> 


</html> 


ПФГПППІПС55ППвдето-ріасеһоаег(пПа ПП 
ПО000000000000000000000000000өП0000 
Онтмі00000000амаѕсгір0ачегу.ј50 
Поѓё.ј$8 0000 ОчегуПћо О000000000000000 
ΠΠΠΠΠΠΠΠΠΠΟΪΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


р0000амазст редород0О00009Ф000)Оцмегупо 
$(document).ready ОПООДО00000НТ7МООО00О 


обо000000000000000Ф0000000Шамазстірі 
ΠΠΠΠΠΠΠ5ΟΘΠΠΠΠΠΗΠΥ1ΠΠΠΠΠΠΠΠΘΠΠΠΠΠΠΠ 
ООО00000000000 x 0 y 000 


т(у) = e™ sin & Jn ΧΠΠΠΠΙΟ, 500] 


ПӨППППоЩППр!о О)ПППППӘӨППр!о()ППППППП 
ΠΗΤΜΙΠΠΠρΙαςεΠοἰαάειῇ1ΠΠΠΙαΓΙΠΠΓΠΠΠΠΠΠΠΗ 
ППИНТМИЦИПИПИПИП 14-4 ППІПІ 


A Simple Flot 


014-4 ΠΠΟΙΠΠΠΠΠΠΠ 


ΠΙΠΠΠΠΠΟοΙΠΙΗΠΠΗΠΠΠΗΠΠΗΠΠΠΗΠΠΗΠΠΠΠΠΠΠΠΗΠΗ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΟΙΠΠΠΠ 
14.4.2П00000000000000000000 


14.2.8 НИ 


ОООДОО000000000000000000000000000000 
ОДООО000000000000000000000000000000000 
055НОДО00000 


9 sudo shutdown -h пом 


ОДООООО0О00000000000000000000000000000000000000 


Linux] 


UUshutdownQQ00000000000000000002 endo 
ΕΕ ΕΕ ΕΕ 


14.3 [| I| 


ОООДО000000000000000000000000 


« DATLI 000 
. 4.7КОППЦ 

• 100ОППП 

• [J[JEED[] 

• 0000 


H14-SQUUUUUUUUUUUDAT119VDDQUUUU+ SV 
DODOHOHDU2UDDH T 110] b ATADOOUDUBOBOULU 
16ронтзІдемордорророросовмороповоб 
РОАТАПУРОДОООООА. 7 КОДООСЕОДДОООО00000 
1000ДрД00009еМмОороррооо0О00000000180 


pin #18 (board) 





DHT11 





1 2 3 4 
(VDD) (DATA) (unused) (GND) 







1000 


μων — ріп #6 (GND) 


—— pin #6 (GND) 


47КО рїп #16 (боага) 
014-5 ΠΠΠΠΡΗΤΙ1ΠΠΠΙ ЕО00000000 


ОДОО0О00000008Н7110-Еб90000000000014- 
еПОо000000000000000000000000 





014-6 ΠΠΠΠΡΗΤΙ1ΠΙ «ΡΠΠΠΠΠΗ 


14.4 ПП 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕΕ 
14.5000 


О00плаї 0000 


аеї та1п(): 
print ‘starting piweather...' 
# create parser 

e 

parser = argparse.ArgumentParser(description="PiWeather..." 
# add expected arguments 


parser.add argument('-- 
ip', dest-'ipAddr', required-True) 


parser.add argument('-- 
port', dest-'portNum', required-True) 


# parse args 


args - parser.parse args() 


4 GPIO setup 
O GPIO.setmode(GPIO.BOARD) 
Ө GPIO.setup(18, GPIO.OUT) 
O GPIO.output(18, False) 
# start server 


Ө run(host=args.ipAddr, port=args.portNum, debug=True) 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ--1ΡΠΠΠΠΠ 
ОДДІРОДО-рогеДОДОДО0ПО000900000РІОПОП 


О00000В0АКООООО00000000000000000000 
9ФПОО0Д0018000000000000000000000-2000 
еррорр0000РаїьерОСЕОПОООО00Д0000Д0ІРО 
ا‎ 


14.4.1 MU 
ОООО0000000000000000 


Ө @route('/getdata', method='GET' ) 

O def getdata(): 

Ө RH, T = Adafruit DHT.read retry(Adafruit DHT.DHT11, 23) 
4 return dictionary 


Ө return {"RH": RH, "T": Т} 


О0000000О0О /де ДатаЦПППӨ®ӨПППИПППППШШПП 
URL /getdata[][]getdata 0рр0000Фп000 
Adafruit рнтороооооооопепореооооооооо 
000000000амазегір р р0000)5О0МД"О0О0000 
П0050М00000—000000000000000000000 
000009екаава 00005О0М000000--00000000 
ОООАНОД0ОД00000100 


14.4.2 [OM 


Ρ!οϊ()ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗΤΜΓΙΠ 
<head> 00000" 000000CSS0 "00000000000 


JavaScriptUQ0U0000 


@route('/plot') 
def plot(): 

Ө return ''' 
<html> 

<head> 


«meta http-equiv="Content - 
Type" content="text/html; charset=utf-8"> 


<title>PiWeather</title> 
<style> 
.demo-placeholder { 
width: 90%; 
height: 50%; 
} 
</style> 
«script language-"javascript" type="text/javascript" 
src="jquery.js"></script> 
«script language-"javascript" type="text/javascript" 
src="jquery.flot.js"></script> 
«script language-"javascript" type="text/javascript" 


src="jquery.flot.time.js"></script> 


plot()QUU/plot URLDBottlen nini RLE 
ΠΠΠΡ!οϊ(})ΠΠΠΠΘΠΠΡ!ΟϊΟ)ΠΠΠΗΤΜΙ ΠΠΠΠΠΠΠΠΠΗ 


ΠΙΠΠΠΠΠΠΠΠΠΝΝΕΡΙΠΠΗΠΠΠΠΠΗΠΠΠΠΠΠΤΜΙΠΙΗΠΗ 
ООС55П0000000Йо8000000000000000014-40 
[rot III IILI III 


ΠΠΠ-«5ο4γ»ΠΠΠΠΠΗΤΜΕΠΠΠΠΠΠ 


<body> 
<div id="header"> 
e <h2>Temperature/Humidity</h2> 


</div> 


<div id="content"> 
<div class="demo-container"> 


о <div id="placeholder" class="demo- 
placeholder"></div> 


</div> 
© <div id="ajax-panel"> </div> 
</div> 
<div> 
о 
<input type="checkbox" id="ckLED" value="on">Enable Lighti 
ng. 
e «span id="data-values"> </span> 


</div> 


</body> 


</html> 


0е0000000000000000Ф0000 
<placeholder> 1 ПОПОООо ама5стр ТІП 


ΠΠΠΘΠΠΠΠΠΠΠΙῬΠ8]8Χ-ρ8ΠπΕΙΠΗΤΜΕΠΠΠΠΠΠΗ 
ΠΠΑΙΑΧΠΠΠΠΘΠΠΠΠΠΠΠΙΡΠΕΚΙΕΡΠΣΠΕςΚΌΟΧΠ 
ΠΠΠΠΠΠΠΠΠΠΠΠΕΕΡΠΠΠΠΠΠΠΠΠΠΙΌΠ4813- 
νθ!165ΠΗΤΜΕΠΠΘΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠΠΠ8ν85ςΓΙρ!ΠΠΠΠΠΠΠΙΌΠΠΠΠΠΠΠΠΠΠΠΠΗ 


О000000000000Јамаѕсгір 000000000000 
ΠΠΠΠΠΠΠΙ ΕΡΠΠΠΠΠΠ« script language 
-“|амуавсгірЕ”...>ЦПППНТМЦЦЦе<һеаа>ППЦ 


Ө $(document).ready(function() { 


// plot options 
Θ var options = { 
series: { 
lines: {show: true}, 
points: (show: true} 
}, 
© grid: {clickable: true}, 
о yaxes: [{min: 0, max: 10011, 


xaxes: [{min: 0, max: 100}], 


// create empty plot 


Θ var plot = $.plot("#placeholder", [[]], options); 


ΠΤΜΙ ΠΠΠΠΠΠΠΠΠΠΠΠΠΘΙ]Πγεααγ()ΠΠΠΗΘ/ΙΠΗ 
урООоріїоп5 ПОДОООООДОДОООДОДОДОДО0000 
er pmo d aad em p ap de 
ПППЗПППППр!о OD n amm PP 
DEBIDO ICD EIE IDLILIDICILILJe p tren ПІ 


ПППППИПППППИГПамазспір TTL] 


// initialize data arrays 
@ var RH = []; 
var T = []; 
var timeStamp = []; 
// get data from server 
O function getData() { 
// AJAX callback 
© function onDataReceived(jsonData) { 
о timeStamp.push(Date()); 
// add RH data 
e RH. push(jsonData.RH) ; 
// removed oldest 


о if (RH.length > 100) ( 


RH.splice(0, 1); 
} 
// add T data 
T.push(jsonData.T); 
// removed oldest 
if (T.length > 100) { 


T.splice(0, 1); 


} 
sl = [1]; 
s2 = []; 


for (var i = 0; i < RH.length; i++) í 
sl.push([i, АН[1]]); 
s2.push([i, T[ill); 

} 

// set to plot 

plot.setData([s1, 521); 

plot.draw(); 

} 


// AJAX error handler 


function onError(){ 


$('#ajax-panel').html('< p» Ajax error! < /р»'); 


// make the AJAX call 
10) $.ajax({ 
url: "getdata", 
type: "GET", 
dataType: "json", 
success: onDataReceived, 


error: onError 


} 


пе rrr mabipaatimesStam роб0О 
ППППИППИППЦППӨӘӨППППде Оаїа()ППППППППППП 
ППППППППӨПІІ Поп ОаїаКесе!уеа()ППШПЦППП 
ОДААХОДОДО2З3ДООО JavaScript IILI IILI L UI] 
ПППППППППППППППППППППОеФага (00000 
опраќѓаВесеімеа()00000000000000АЈАХО00 


onDataReceived()[| JavaScript Date 
ΠΠΠΠΠΠΠΠΠΘ/ΙΠΠΠΠΠΠΠΠΠΠΙ5οηαἴα[[[Π[Πῃ 
Οπαϊαβεςείνεαι()[Π[Πθ[ ([Π[ΠΙΠΙΠΙΠΙΠΙΠΙΠΙΠΙΠΙΠΙΠΙΠΠΙΠ 
ΠΠΠΘΠΠΠΗΠΠΗΗΠΠΗΠΠ1ΘΟΘΠΠΠΠΠΗΠΗΗΠΠΠΠΠΠΠΠΠΗ 
0000002 2000Агаміпо ДООО00000000000000 
ЗО00000000 


Lea E pd ad add eret n p rr 
О00000000000030000000 


Ш o, АН о], Lii, АН 11], ...1, П ѓо, Tob [11 
‚Т1]] 


Le rpg m b 


ΠΘΠΠΠΠΠΠΠΠΠΠΠΠΑΙΑΧΠΠΠΠΠΠΠΠΠΙΡΠΘ|8Χ- 
Ρ8ΠπΕε!ΠΗΤΜΙΕΠΠΠΠΠΠΠΠΠΘΦΘΠΠΠΠΠΠΠΑΙΑΧΠΠΠΠ 
ΠΠΟΗΙΠσ6[48ϊ8ΠΠΒΟΓΙΕΠΠΠΠΠΠΠΠΗΤΤΡΠ 
ΟΕΤΟΠΠΠΠΙ5ΟΠΠΠΠΠΠΠΠΠΠΠΠΑΙΑΧΠΠΠΠΠ 
ОпрааКесемеа о ПП ППППоОпПЕГГОГОП 
Ui μμ 


14.4.3 update()[][] 
[LIEIEIEIEILILL]getD ata O[Jupd ate OD 


// define an update function 
function update() { 
// get data 
e getData(); 
// set timeout 


6 setTimeout(update, 1006); 


// call update 


update(); 


update([nummeumgetbata (QD ern 
JavaScript[]setTimeout()[]L]L] 2O OO[TELILLIE]LILTLH 
Порае аса ООО000000000000000000 


14.4.4 ППрЕрЏама5бстркооо 


ΠΠΠΠΠΙΕΡΠΠΠΠΙανα5εΓΙρΙΠΗΠΠΠΠΗΝΝΕΡ/ΙΠΠΠΗ 
ППА)АХТПП 
// define the click handler for the LED control button 
€ $('#ckLED').click(function() { 
6 маг isChecked = $("#ckLED").is(":checked") 2 1:0; 
θ $.а]ах({ 
url: '/ledctrl', 
type: 'POST', 
data: ( strID:'ckLED', strState:isChecked j 
}); 
}); 


ΠΘΙ[ΙΠΠΠΠΗΤΜΙΤΙΗΠΠΠΙΒΠΕΚΙΕΡΙΠΙΠΗΠΠΠΠΠΗΠΗ 
ООО000000000000000000000000000000000 
ΠΠΠΠΠΘ/ΠΠΠΘ/ΙΠΠΠΠΠυΒΙ. /ledctriQHTTPOOOO 
РОЗ ТПОПААХООООДООООО000000 


DUAJAXUBUDBOBOUOHOHOUDUGPIOLDDUCLI 
00000 
Ө @route('/ledctrl', method='POST' ) 
def ledctrl(): 
O val = request.forms.get('strState' ) 
© on = bool(int(val)) 


Ө  GPIO.output(18, on) 


перпорпові /ledctriQBottleQQ000000000 

Іеасё()00000000ӨПОООВоШегеаиеѕ 0 
ПО00000000055аёеПоо00000000ө000000 
MEL PIO.output(L[LDm mp p 


14.4.5 (0000 


ООДОД00000000000 ЯоєО000000000000000 
ОООД0000ріосо000000Ооріїоп5000 
у" ble:true[]DO aman n p p CIC 


$("Zplaceholder").bind("plotclick", function (event, pos, i 
tem) { 


if (item) { 


° 


plot.highlight(item.series, item.datapoint); 


e var strData = ' [Clicked Data: ' + 


timeStamp[item.dataIndex] + ': Т = ' + 


T[item.dataIndex] + ', RH = ' + RH[item.dataIndex] 
+ ']'; 

e $('#data-values').html(strData); 

} 

7); 

y); 
«/script» 
</head> 


прорероопо моћну ВЕОООДОДОООДОД0000 
О00ӨПП000000000000000ааѓа-уаіџеП 
НТМіТОО00000000#о#0000000ето00000 
дО даса паех)р)0000000000000Огеаау 0000 


ООООО000000000000000000000000000 
НТМШПДӨП 


О00РуєһопО00000амаѕсгірЕО0000000 
Bottle UUU0WebQ 0000 avaScriptQ 000000 


@route('/<filename:re:.*\.js>') 
def javascripts(filename) : 


return static file(filename, root='flot') 


ΠΠΠΠΠΠΒοϊἑΠΠΠΠΠΠΠΠΟΘ/ΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΕΠΗ 


14.5 [ΠΠ] 


ПІ 
https://github.com/electronut/pp/tree/mast 


епрмеа ће  рмеа ег.руУППООООООООООО 
from bottle import route, run, request, response 
from bottle import static file 
import random, argparse 
import RPi.GPIO as GPIO 
from time import sleep 


import Adafruit DHT 


@route( '/hello') 
def hello(): 


return "Hello Bottle World!" 
@route('/<filename:re:.*\.js>') 
def javascripts(filename): 


return static file(filename, root='flot') 


@route('/plot') 


def plot(): 
return ''' 

«html» 

<head> 


«meta http-equiv="Content - 
Type" content="text/html; charset=utf-8"> 


<title>PiWeather</title> 

<style> 

.demo-placeholder { 

width: 90%; 

height: 50%; 

} 

</style> 

«script language-"javascript" type="text/javascript" 
src="jquery.js"></script> 

«script language-"javascript" type="text/javascript" 
src="jquery.flot.js"></script> 

«script language-"javascript" type="text/javascript" 
src="jquery.flot.time.js"></script> 


«script language-"javascript" type="text/javascript"> 


$(document).ready(function() ή 


// plot options 
var options = { 
series: { 
lines: {show: true}, 
points: {show: true} 
), 
grid: {clickable: true}, 
yaxes: [{min: 0, max: 100}], 


xaxes: [{min: 0, max: 100}], 


// create empty plot 


var plot = $.plot("#placeholder", [[]], options); 


// initialize data arrays 
var RH = []; 
var T = []; 
var timeStamp = []; 
// get data from server 
function getData() { 

// AJAX callback 


function onDataReceived(jsonData) { 


timeStamp.push(Date()); 
// add RH data 
RH.push(jsonData.RH) ; 
// removed oldest 
if (RH.length > 100) { 
RH.splice(0, 1); 
} 
// add T data 
T.push(jsonData.T); 
// removed oldest 
if (T.length > 100) { 


T.splice(0, 1); 


} 
51 = []; 
s2 = []; 


for (var i = 0; i < RH.length; i++) { 
sl.push([i, RH[i]]); 
s2.push([i, ТІ111); 

} 

// set to plot 

plot.setData([sl, 521); 


plot.draw(); 


// AJAX error handler 
function onError(){ 


$('#ajax-panel').html('<p><strong>Ajax error! 
</strong> «/р»'); 


} 


// make the AJAX call 
$.ajax({ 
url: "getdata", 
type: "GET", 
dataType: "json", 
success: onDataReceived, 
error: onError 


}); 


// define an update function 
function update() { 

// get data 

getData(); 

// set timeout 


setTimeout(update, 1000); 


// call update 


update(); 


// define click handler for LED control button 
$('ZckLED').click(function() { 
var isChecked = $("#ckLED").is(":checked") ? 1:0; 
$.ajax(1 
url: '/ledctrl', 
type: 'POST', 
data: ( strID:'ckLED', strState:isChecked } 
}); 
}); 


$("Zplaceholder").bind("plotclick", function (event, pos, i 
tem) { 


if (item) { 
plot.highlight(item.series, item.datapoint); 
var strData = ' [Clicked Data: ' + 


timeStamp[item.dataIndex] + ': Т = ' + 


T[item.dataIndex] + ', RH = ' + RH[item.dataIndex] 


к]; 


$('Zdata-values').html(strData); 


} 
}); 
}); 
</script> 
</head> 
<body> 
<div id="header"> 
<h2>Temperature/Humidity</h2> 
</div> 
«div id="content"> 
«div class="demo-container"> 
<div id="placeholder" class="demo-placeholder"> 
</div> 


</div> 
<div id="ajax-panel"> </div> 
</div> 


<div> 


<input type="checkbox" id="ckLED" value="on">Enable Lightin 
g. 


<span id="data-values"> </span> 


</div> 


</body> 


</html> 


@route('/getdata', method='GET' ) 

def getdata(): 
RH, Т = Adafruit DHT.read retry(Adafruit DHT.DHT11, 23) 
# return dictionary 


return {"RH": RH, "T": T} 


@route('/ledctrl', method='POST' ) 
def ledctrl(): 
val = request.forms.get('strState') 
on = bool(int(val)) 


GPIO.output (18, on) 


# main() function 


def main(): 


print ‘starting piweather...' 


# create parser 


parser = argparse.ArgumentParser(description="PiWeather... 


) 
# add expected arguments 


parser.add argument('-- 
ip', dest-'ipAddr', required-True) 


parser.add argument('-- 
port', dest-'portNum', required-True) 


# parse args 


args - parser.parse args() 


4 GPIO setup 
GPIO.setmode(GPIO.BOARD) 
GPIO.setup(18, GPIO.OUT) 
GPIO.output(18, False) 
# start server 


run(host=args.ipAddr, port=args.portNum, debug=True) 


# call main 


if name == ' main ': 


main() 


14.6 (OU 


DUDODOHDH T L1ECDUDUDUUS S НООООООО 
ОО0000000000000000001РО000000 


$ sudo python piweather.py --ip 192.168.x.x --port xxx 


ОДОО00000000000000000001РОДОО0000000 


http://192.168. x.x:port/ plot 


0000000 + 4- 7000000 


Temperature/Humidity 


100 
М) 
@ 
У Y Όρου ¥990.0000000550000 90065000009000000 0600000 











0 10 ) X 40 X 6 
@ Enable Lighting. (Clicked Data: Tue Jun 24 2014 15:55:30 GMT40530 (IST): T = 18, RH = 45] 


014-7 piweather.pyQU0000000000 


ООО000000000000000000000000000000000 
О00000100000000000000000000000000 


Enable LightinaUUUUUUUUULEDUL 
14.7 ||| 


ООООО0000000000000000УУеБОО00000000 
ООО0000000000000 


e 000000 

- ΠΠΠΠΠΠΘΡΙΘΠΠΠΗΠΠΠΗ 

• ДОНТ11000000000 

• (Python Web[j(]Bottle[]( | Web[]EHLL 
• ПЏамазспр По ООООО 

• ПО000--00000000 

• ΠΠΝΥΕΌΠΠΠΠΠΠΗ 


14.8 [| 


ООО0000000000000 


1 ОДО0Д00000000000000000000000000000 
ΠΤΠΒΗΠΠΠΠΠΠΠΠΠΠ/εχρο!Γ{ΠΠΠΒοΠΕΠΠΠΠΠΠ 
ΠΠΠΠΠΠΟΘ5νΠΠΠΠΠΠΠ/ρ!οϊΠΠΠΠΠΠΠΗΤΜΕΠΠΠΠ 
00000" Export ООООООО0О0О0000А/АХОДОООО 
[Dex portd ОО00000000000000С5У0000000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΕΕΕ 


20000000Н171100001000О000000000000000 
ООО000000000000000000000000000т08нр 


ОООООД000000000000000000Н 7МЕ0000000 
ОЦ/амазстірьЇооД000000000000000200000 
ООО0О0000000000000000000000000000000 


ОООО00000000000000000000000000000000 
SQLite[] 


ППА 0000 





ОДОДО000О00000РУспопОДОДООО0000000 
ОбО00240000000000000000000000000000 
ο ο μμ 2.7.8ЦРуїһоп 3.3.3000 


Α.Σ. UU 


[Hd] E]https://github.com/electronut/pp/[]L]LIL] 
Ор00000000000000Дом/піоаа ΖΙΡΠΠΠΠΠΠΠΗ 


ШПШШШШШШШШШШШШШШШШШШШШСОТШПТТтОПЦППЦПЦПШШШШШШ 
Црр-тасїег/соттопЦПЦЦЦРҮТНОМРАТНЦПП 


О000000000000000РуёһпопП00 


ОМіпаомѕПО0000000РҮТНОМРАТНОП0000 
О00000000000005 ХОО000000000000 


Π.ΡΓΟΠΙΕΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 


export PYTHONPATH=$PYTHONPATH: path to common folder 


Шпихо0000000005 ΧΠΠΠΠ 
П.БавһгсП.фавһ ргоћіе.сѕћгс/.Іодіп00000 


ППесһћо $5НЕШЦЦЦ000ѕ=һе1р 


ΠΠΠΠΠΠΠΠΓΜΝΙΠΩΟνν5ΠΟ5 ХОшпихцЦОРуЁћоп 
ΕΕ ΕΕ ΕΕΕΕΕΕΕ 


A.2 (МЛлпаом5 || 


0000 https://www.python.org/download/[][] 
ПООРућопр 


A.2.1 [|J]JGLFW 


О000000ПОрепсе (Да D0UUUUUUG ΓΝΝΠΗΠΗΗΗ 
http://www.glfw. org/download.htm!l ΠΠΠ 


Πνν!π4ονν5ΠΠΠΠΟΙΕΝΝΠΠΠΩΙ ΕΝ LIBRARYODO 
ППППППППЕЯЕ Environment Variables[ I] 
ΠΠα!ΐνν3.α!!ΠΠΠΠΠΠΠΠΟΙΕΝΜΠΡΥ΄ΠΟΠΠΠΠΠΠΠ 
ΠΠΠΠΠΠΠΠΠςΑ glfw-3.0.4.bin.WIN32\lib- 
msvc120\glfw3.dll[] 


ΠΡγ{ΠοΠΠΠΠΘΙΓΝΝΠΠΠΠΠΠΠΡΥ9!ἵννΠΠΠΠΠΠΠΠ 
ΡγζΠΟΠΠΠΠΠΠΠΠ9!ΐνν.ργΠΠΠΠΠΠΠΡΥ9!ΝνΠΠΠΠΠ 








ООООО000000Осоттопрооооооооооооо0о00 
О0000000№рѕ://9іћир.сот/ 
rougier/pyglfw/[] 


ООО0О0000000000000000000000000000000 
О000000000000000006РЧрП 


A.2.2 ΠΠΠΠΠΠΗΠΠΠΠΠΠΠΠΠΗ 


[IWindowsli|[|[I[IIIIIPythoniII IILI m pap 
000000000000000000000000000000032064 
ο μμ 


pyaudio 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#pyaudio 





pyserial 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#pyserial 


scipy 


http://www.Ifd.uci.edu/- gohlke/pythonlibs/ 
#scipy 





http://sourceforge.net/projects/scipy/files/s 
cipy/ 


numpy 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#numpy 


http://sourceforge.net/projects/numpy/files/ 
NumPy/ 


pygame 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#pygame 


Pillow 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#pillow 


https://pypi.python.org/pypi/Pillow/2.5.0#d 
ownloads 


pyopengl 


http://www.Ifd.uci.edu/- gohlke/pythonlibs/ 
#pyopengl 


matplotlib 


http://www.Ifd.uci.edu/- gohlke/pythonlibs/ 
4 matplotlib 


matplotlib[](]]l]dateutil[]pytz[]pyparsing[] 
51ΧΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


dateutil 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#python-dateutil 


pytz 


http://www.Ifd.uci.edu/- gohlke/pythonlibs/ 
#pytz 


pyparsing 


http://www. lfd.uci.edu/~gohlke/pythonlibs/ 
#pyparsing 


Six 


http://www.Ifd.uci.edu/- gohlke/pythonlibs/ 


#51Х 


A.2.3 |ІІІ! 


ΠΙΠΠΠΠΠΠΠΠΠΠΠΠΗΝΝΙΠάοννΞϱΙΓΙΠΠΠΠΠΠΠΠΠΠΠΗ 
ОСО0000000000 








https://docs.python.org/2/install/index.html 
#gnu-c-cygwin-mingw]000000 
http://www.scipy.org/install.html [][][]L]L] 
ΡΥΓΠΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


А.З [Πο ХЦ 


000005 ΧΠΠΗΡΥΕΠΟΠΠΠΠΠΠΠΠΠΠΠΗ 
A.3.1 [][|Xcode[|MacPorts 


ΠΠΠΠΠΗΧΣΟΩΕΠΠΠΠΠΠΑΡΡ 5ἵογα[[ΠΠΠΠΠΠΙΠΠΙΠ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΑΡΡΙΕΠΠΠΠΠΠΠΠΊρ5Η// 
демеіорег.арріе.сот/0000000Хсодецр 
ХсодеП000000000000000000МасрРогёѕ[0 
ΠΠΠΗΜΘΣΡΟΠ5ΠΠ 
Ihttp://guide.macports.org/#installing.xco 


αΕΓΙΠΙΠΠΠΠΠΠΠΠΗΠΠΗΠΠΗΠΠΠΠΠΗΠΗ 


МасРогі5 ДПОДОРУЄПОПООООООООООООДОООДО 
ο. аи 


A.3.2 ΠΠΠΗ 


ЦПЦМасРогї5$ЦЦЦЦЦИТегтїпа!ЦПЦПЦПЦПЦППППППрогї 
ΕΕ ΕΕΕΕΕΕΕΕΕΕΕ 


О00000000000000Руёћопо000 


$ port select --list python 


ОО0000РУуєћоп000000000000000000 
ΡγίΠΟΠΠΠΜΘΕΡΟΠ5ΠΠΠΠΠΠΠΠΠΠΠΡΥΙΠΟΠΠΗ 
2.700 


$ port select --set python python27 


ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕ 


sudo port install py27-numpy 

sudo port install py27-Pillow 
sudo port install py27-matplotlib 
sudo port install py27-opengl 
sudo port install glfw 

sudo port install py27-scipy 

sudo port install py27-pyaudio 
sudo port install py27-serial 


sudo port install py27-game 


МасРог5$ЦП[ЦПРуїһопЦППЦП/орїЛоса!/ПППППП 
Π.ΡΓΟΠΙΕΠΠΠΡΑΤΗΉΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
ΡγΥΓΠΟΠΠΠΠΠΠΠΠΠΠΗ 


PATH=/opt/local/Library/Frameworks/Python.framework/Version 
s/2.7/bin:$PATH 


export PATH 


ΠΠΠΠΠΠΠΠΠΡντποπΠΠΠΠΠΠΠΠΠΠΠΠ 
A.4 ΠΕίπυΧΠΠΗ 


ПіпихООДОРУЄПОПООООООООДОООДОО00000000 
Шпихо1О00000000ргроО000000000000ртр00 


ОО00000000 
http://pip.readthedocs.org/en/latest/installi 


па ЛЕ 
ΠΠΠΡΙΡΠΠΠΠΠΠΠΠΠΗ 


sudo pip install matplotlib 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ.9ΖΠ.Ζ!ΡΠΠΠΠΠΠ 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΡΕ 


sudo python setup.py install 


ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


HUB б000000 





ООО000000000000000000000000000000000 
ОООД0000000000000000000000 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ“ΠΠΠΠΠΠΠΡΙΥΠΠΠΠΠΠΠΠ 
ООО000000000000000000000000000000000 
ООО000000000000000000000000000000000 
ООООО0000000000000000000000000000 


В.1 (ІП 


ООО000000000000000000000000000000000 
ΕΕ ΕΕ ΕΕ 


ОДОО0О000000000000000000008-1000000 
ООО00000000000000 


BREAD BUARD 
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ПВ-1 ДООДО000000000 
В.1.1 (LL 


ООО000000000000000000000000000000000 
Tiri ο μμ μου 


B.1.2 ΠΠΠΠΠΙΡΒΠ 


СоКО0000000000000000000000000000000 
ΗΗΠΗΠΠΗ 


B.1.3 ΠΠΠΠΠΙΟΠ 


ІСОООДОДОО000001С00000000000000000000 
ОДІСООДОД0О0000000000000000000000000 
ОДООДОД00000000000000000000001С0555П 
ΕΕ ΕΕ ΕΕ 


B.1.4 ΠΠΠΠΠΠΡΟΒΗ 


ΠΠΙΠΗΠΗΠΗΠΠΠΠΗΠΗΠΠΗΠΠΗΠΗΡΟΒΙΗΗΠΗΗΜΗ 
РСВОООО0ОООООО0О000000000000000000000 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ 


OOOO“ DEP BOLD BUD ОПОДОРЄВОДОООДОДООДО 
ΗΗΠΗΠΗ 


B.1.5 ΠΠ 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΡΕ 


B.1.6 || 


ООО0О0000000000000000000000000000000 
000000000000000002. 7к0000002. 70000 

2 ТООПОДОДОО0О0000000000000000000000 
ΗΠΗΠΗΠΗΗΠΗΗ 


B.1.7 ΠΠΠΠΠΠΙΕΡΗ 


-ЕОП000000000000000000-ЕОПО00000000 
ООО0О000Д0000000000000000000000000000 
ОДООО0О0000000000С(Е000000007007000 


B.1.8 {ΠΠ 


ООО0О0000000000000000000000000000000 
ООО00000000000иғО0000000000000000 


B.1.9 [ПЦ 


ОООДОО000000000000000000000000000000 
АСПООООРЄПОООДОООДОО0О0000000000000 
ОООО00000000000000000000 


B.1.10 ПІП 


ООО0О0000000000000000000000000000000 
ООО0О0000000000000000000000000000000 
ОДД0008/ ООДОО000000000000МОЗРЕТОГО 


ПШШППШИПИПШИШППВ)ТИШШИИШШИИШИИШШИШИШИШИШ 
МОѕЅҒЕТОООООО000000000000000000008]ЈТ 
ОДОООДПОМОЗ5РЕТОДООООООООООО000000000 
ОДОООДО00000000000000С(Е2Е000000 


B.1.11 П/П 


О00000003090000000000000000000000000 
ООО000000000000000 


B.2 ПШ 


О00000000000000000000000000000000В-2 
ООООО0000000000000000000000000000000 





ПВ-2 ПООООООООО000000000000000000000000000000 
ОО00000000 


В.2.1 000 


ООО000000000000000000000000000000000 
ООО0О0000000000000000000000000000000 
00000000 


B.2.2 [IE 


ΠΠΠΙΠΙΠΠΗΠΠΠΠΠΠΗΗΗΠΗΠΠΠΠΗΠΗΗΡΟΒΗΠΗΠΗ 
ООООО0000000000000000000000000000000 
ОДОДОДОД0000000000000000000000ОРС8О0000 
ООО000000000000000000000000000000000 
ОООО00000000000000000000 


B.2.3 (LL 


ОДОО000000000000000000000000000000000 
ОДОО000000000000000000000000000000000 
ἱ μυ 


0В-2П000000000000000000000000000000 
ОО000000000РЄВОООО0000000000000000000 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΡΕ 


В.з OUOU 


ООО000000000000000000000000000000000 
ООО0О0000000000000000000000000000000 
ОДОДОДОДОРЄВОДОДООООООДОООООО000000000 
ΕΕ ΕΕ ΕΕ ΕΕ III 


О000000РСВООО000000000000000000РСВП 
ООО000000000000000000000000000000000 
ОДОДОО0Д0000000000008-З00000000000000 
ОДОООДОРСВОДО 


ППППППППППРСВППППППППППППППППППППЇПП 
ΠΠΠΠΠΠΠΠΠΡΟΒΠΠΠΠΠΠΠΠΠΠΠΕΑΩΙΕ!«!ῃ 
KiCad ІЗ'ПЕАС ΕΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΡΟΒΠΠΠΙΟΠΠχ8Π 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
ΠΠΕΑΟΙ ΕΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠ (0 























DB-3 ΠΠΠΠΠΠΠΠΠΠΠΠΠΡΟΒΠΠ 


EAG СЕПДООДОООООО0О00О00000000000000 
ОО0000000000000000000ЕАСГЕПОО0000000 
О00000000000000ЕАСГЕПОООООО00000000 
О000000ЕАСГСЕПОО000000РСВООО00000000 
О0000000000000000000000РСВО000000000 
О00ЕАСГЕЦ00В-4000000000000000000000 
ПЕАСГЕПОО00000УоитТирерВО000000 


























08-4 ΠΠΕΑαΙ ΕΠΠΗΠΠΠΠΗΠΠΠΠΡΟΒΠΗ 


ΠΠΠΡΕΒΠΠΠΠΠΠΠΠΠΠΠΠΠΡΕΒΠΠΠΠΠΠΠΠΠΠ 23 
О0000000000000000000000РСВО000000000 
ΠΠΘ6ΕΓΌΘΠΠΠΠΠΠΠΠΠΠΠΠΕΑΘΙΕΠΠΠΠΠΠΠΠΠΠΠΠ 


0 “ΠΠΗΠΗΠΠΡΟΒΠΠΠΗΠΠΗΠΗΠΠΠΠΟΣΗ Рагк 
Zn 


ΠΠΠΡΟΒΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
ΠΠΠΠ3ΡΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ2Ρ irj 
3D ΙΞΊΠΠΠΠΠΠΠΠΠΠΠΠΠΒΙΕΙ DecibelsQOOO000 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ «lm 


В.4 (LL 


ООО000000000000000000000000000000000 
ОДОО000000000000000000000000000000000 
О00000000000000Агаиіторо00000000000 
ООО0О0000000000000000000000000000000 
ООО000000000000000 


ОООО00000000000000000000000000000000 
ДОДОДОДОДООДіп5ігиссаріез) 
http://www.instructables.com/ 000000000 
ОООООДО00000000001У00000000 


[1] Ор000000000000Райі ScherzrSimon 


Monk[][]Practical Electronics for Inventors[] 
[]3H]]]McGraw-Hill[]20 13(][] 


[2] CadSoft EAGLE ΡΟΒΠΠΠΠΗ 
http://www.cadsoftusa.com/eagle-pcb- 


design-software/ П 


[3] KiCad ЕРАПЦООО http://www.kicad- 
pcb.org/display/KICAD/KiCad--EDA 4 Softwa 
re+Suite/[] 


[4] Jeremy Blumf]“ Tutorial 1 for Eagle[] 
Schematic Design"[]YouTube[]2012[]6[]14 
ΠΠΠΠ{ροῃ// www.youtube.com/watch? 
v=1AXwjZoyNno|] 


(З ΠΠΠΗΠΡΟΒΠΠΠΠΠΠΠΗ 


https://embeddedinn.wordpress.com/tutori 
als/home-made-sigle-sided-pcbs/[](L] 


[6] Akash Patel[]"Generating Gerber files 
from EAGLE"[]YouTube[]2010[]4[]7 [TL ]L] 
https://www.youtube.com/watch? 
v=B_SbQeF83XU[] 


[7] https://oshpark.com/[] 


[8] ΠΠΖΡΠΠΠΠΠΠΙηΚΞεαΡρε/[ 
http://www.inkscape.org/en/ [] 


[9] ППЗОППППППӘкекссһӘрП 
http://www.sketchup.comYj [] 





[10] Rich Decibels[]"Laser-Cut Project Box 
Tutorial[]" Ponoko[]2011[]8[]O[T[]L] 
http://support.ponoko.com/entries/203444 
37-Laser- cut-project-box-tutorial/[] 


ППС ППППППИПП 





Обр002400000000000000000000000000000 
ООО0О0000000000000000000000000000000 
ОООО000000000000000 


C.1 || Wi-Fi 


002. 400О00ДОДО0ДМ/1-РІОООООДОДОДДОДОУМ- 
РІООДООДОДОООВОДОДОДОООО000000000000000 
паподО0000000000 


$ sudo nano /etc/wpa_supplicant/wpa_supplicant. conf 


ОО000000000 


ctrl interface-DIR-/var/run/wpa supplicant GROUP=netdev 


update config-1 


network={ 
ssid- your-WiFi-network-name 
р5К= your-password 
proto-RSN 
key mgmt-WPA-PSK 
pairwise=TKIP 
auth_alg=OPEN 


} 


ΠΠΠΗΠΠΠΗΠΗΞΣΙάΓΠΙΡΞΚΙΠΠΙΠΠΠΝΝΙ-ΕΙΠΠΗΠΠΠΠΠΗ 
ο μμ.» 


C.2 ΙΠΠΠΠΠΠΠΠ 


ΠΠΠΠΠΠΠΗΠΠΡΙΠαΓΠΠΠΠΠΠΠΗΠΠΠΠΗΗΠΠΠΗΠΠΗ 
ріод 00000000 


ф ping 192.168.4.32 


PING 192.168.4.32 (192.168.4.32): 56 data bytes 


64 bytes from 192.168.4.32: icmp seq-0 ttl=64 time=13.677 m 
5 


64 bytes from 192.168.4.32: icmp 5ед-1 ttl=64 time=8.277 ms 
64 bytes from 192.168.4.32: icmp seq=2 ttl=64 time=9.313 ms 


-- Snip -- 


ΠΠΡΙΠΟΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ 
[| Request timeout...” ПОДООДОД0000000 


С.З ДОМА ОО0ОО00000 


ООО0000000ріпабОО0005НН00000058 Wi-Fi 
ООО0О0000000000000000000000000000000 
ОООО000000000000000 


5 sudo nano /etc/modprobe.d/8192cu.conf 


ООО000000000000 


я disable power management 


options 8192cu rtw_power_mgnt=0 


ΠΙΠΠΠΠΠΠΠΝΝΙ-ΕΙΠΠΠΠΠΠΠΠΠΠΠΗΠΗ 


C.4 ΙΠΠΠΠΠΠΠΠΠΠ 


О000000000000000000000000000000г5упс0 
ООООО0000000000000000000000000000000 
ΕΕ ΕΕ ΕΕ ΕΒΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ПВ 


ООО000000000000000000000000900000-п0 
ОД000"700700000геупс ПОООООДОДОДОДОДОДО 
ООО00000Д0000000000000000000000000000 
000000000005 XIII 


#!/bin/bash 


echo Backing up ЕРІ ΜΉ... 


4 set this to Raspberry Pi IP address 


РІ ADDR="192.168.4.31" 


4 set this to the Raspberry Pi code directory 
# note that the trailing slash (/) is important with rsync 


PI DIR="code/" 


4 set this to the local code (backup) directory 


BKUP DIR="/Users/mahesh/code/rpil/" 


# run rsync 
# use this first to test: 
я rsync -uvrn pi@$PI ADDR:$PI DIR $BKUP DIR 


rsync -uvr pi@$PI ADDR:$PI DIR $BKUP DIR 


echo ... 


echo done. 


# play sound (for OS X only) 


afplay /System/Library/Sounds/Basso. aiff 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΙπυχΠΟ5 ХО00гѕупср 
ΠΠΠννίπάονν5ΠΠΠΠΠΠΘγογης!1Π 


C.5 ΙΠΠΠΠΠΠΠΗΠΠ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ5ΡΗ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ5ΡΌΡΠΠΠΠΠΠΠΠΠΗ 
ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 
5[8ΟΚΕΧΟΠΘΠΩΘΠΠΠΠΠΠΠΠΠΠ!«!ΠΠΠΠΙἰπυχῃΠ 
OS ΧΠΠΠΩΩΠΠΠΠΠΜΠαονν5ΠΠΠΝνΙΠ32 Disk 
ΙΠΙΒΩΘΓΠΠΠ 


C.6 [IISSHI[| I IILI 


LO 40000000000000055НОБОООДОО0000000 
ΕΕ ΕΕ ΕΕΣ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ILI ΕΕΕΕΕΕ 
SSH(UUssh-keygen(QU000000000/0000000 
О00000000000000000000005 XI ILinuxI III! 
ОДДООО00000М// пом ОДОПРИТ ТУОДОДОДООО 


(З (OOO 0000000 0000000 0000000000!P O00 


$ ssh-keygen 
Generating public/private rsa key pair. 


Enter file in which to save the key (/Users/xxx/.ssh/id rsa 
): 


Enter passphrase (empty for no passphrase): 
Enter same passphrase again: 


Your identification has been saved in /Users/xxx/.ssh/id rs 
a. 


Your public key has been saved in /Users/xxx/.ssh/id rsa.pu 
b. 


The key fingerprint is: 


-- snip -- 


ООО0000000000000 


$ scp -/.ssh/id rsa.pub pi@192.168.4.32:.ssh/ 


The authenticity of host '192.168.4.32 (192.168.4.32)' can' 
t be established. 


RSA key fingerprint is fl:ab:07:e7:dc:2e:f1:37:1b:6f:9b:66: 
85:2а:33:а7. 


Are you sure you want to continue connecting (yes/no)? yes 


Warning: Permanently added '192.168.4.32' (RSA) to the list 
of known hosts. 


pi@192.168.4.32's password: 


id rsa.pub 100% 398 
0.4KB/s 00:00 


ΗΠΗΠΗΠΉΗΠΗΗ 


$ ssh р1@192.168.4.32 


pi@192.168.4.32's password: 


$ cd .ssh 
$ 15 


id rsa.pub known hosts 


$ cat id rsa.pub »» authorized keys 
$ 15 
authorized keys id rsa.pub known hosts 


$ Logout 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ55Π- 
keyg en 0000000000000 0000000000000 
ООДДДД0000000000000055НнО0000000000 
[“Working with SSH Key Passphrases” 00 


ΠΠΠΠΟίτΗυρΠΠ!1ῃ 
C.7 ΙΠΠΠΠΠΗ 


ΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠΠ5“ΠΠΠΠΠΠΠΠΠΗ 
О00000000000000005000000000001080000 
ЗОДУООООО0ОДО0О00000000000000000000 
raspistil ОООООДООО0О000000000000000000 


0000000000000000 ΠΠ 
raspistill gy | 
C.8 ПИШ 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ. 
ee Web Ῥ65!ϑΠΠΠΠΠΠΠΠΠΠΠΠΠΗ 


C.9 0000 
ПППППППППППППППППППППППППППППру5хПП 
ПООРУћоп0000000 {810000000000 


$ 
wget https://pypi. python. org/packages/source/p/pyttsx/pytts 
1.1.tar.gz 


$ gunzip pyttsx-1.1.tar.gz 
$ tar -xf pyttsx-1.1.tar 
$ cd pyttsx-1.1/ 


$ sudo python setup.py install 


ΠΠΗΠΠΠεΞρεακΠΗΠΠΠΗ 


$ sudo apt-get install espeak 


ОДООД0000000000000000000000 


import sy: 


import pyttsx 


# main() function 

def main(): 
# use sys.argv if needed 
print 'running speech-test.py...' 
engine - pyttsx.init() 
str - "I speak. Therefore. I am. " 
if len(sys.argv) » 1: 

str - sys.argv[1] 

engine.say(str) 


engine.runAndWait() 


# call main 
if name == ' main ': 


main() 


C.10 (HDMI 


О000ноМІП00000000000000000000000000 
ОО000000000000000005000000000000 
config.txtQ 0000 


hdmi_force_hotplug=1 


ОДОО0О000000000000Н0МІ00000 


C.11 [II 


ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕΕΕΕΕΕΕΝ 
ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ ΕΕ TEILTE TE TE] 
micro ОФВОДОДОДОДОДОДАпкКег Astro Mini 


30ΟΟΟΠΙΑΠ/ΠΠΗΠΠΗΠΠΠΠΠΗΠΠΠΗΠΠΗΠΠΗΠΠΖΟΠΠΗ 


C.12 [IB 
ИШТ 


$ cat /proc/cpuinfo 


ΕΕ ΕΕ ΕΕ 


processor : ϐ 

model name : ARMv6-compatible processor rev 7 (v6L) 
BogoMIPS : 2.00 

Features : swp half thumb fastmult vfp edsp java tls 
CPU implementer : 0x41 

CPU architecture: 7 

CPU variant : 0x0 

CPU part : 0xb76 


CPU revision : 7 


Hardware : BCM2708 


Revision : 000f 


Serial : 00000000364a6f1c 


О0р0000000гемізіоп ДОООП 
http://elinux.org/RPi HardwareHistory[][J[] 
лаа шаа 
OBOL 


[1] 0000http://grsync-win.sourceforge.net/0 
[]grsync - ПДМіпаом5Ддгзвупс ПДД 


(21 ППППобаск ЕхсһапдеЦППППППППП5ОППП 
ПП“Ноу/ do І backup my Raspberry Pi?”[] 
Stack Ехсһапде[ 
http://raspberrypi.stackexchange.com/ques 
tions/311/how-do-i-backup-my-raspberry- 
рі/П) 

[3] “Πονν to Create SSH Keys with PuTTY to 
Connect to а VPS"[lDigitalOcean[]2013[]7[] 
19000 https://www.digitalocean.com/ 


community/tutorials/how-to-create-ssh- 
keys-with -putty-to-connect-to-a-vps / [] 


[4] "Working with SSH Key Passphrases"[] 
GitHub JUL 





https://help.github.com/articles/working- 
with-ssh-key-passphrases/[] 


[5] 0000000000000000000 
http://www.raspberrypi.org/product/camera 
-module/ [] 


[6] The RaspberryPiGuy[]"Raspberry Pi- 
Camera Tutorial[]”YouTube[]2013[]5[]26[][][] 
https[]// www.youtube.com/watch? 
v=T8T6S5eFpqE[] 


[7] "Raspberry Pi—Getting Audio Working” | 
CAGE Web ДрезідпД2013П2П9ПДПРЕД// 
cagewebdev.com/index.php/ raspberry-pi- 
getting-audio-working/[] 


[8] ΠΠΠΠΠ{ρ5:// github.com/parente/pyttsx/ 
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12 Arduino[T] 
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12.5 Python] 
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